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1. — Introduction. 


The inquiry of this paper begins with the commutation rule for angular 
momentum and proceeds from it to analyse the algebra generated by the com- 
ponents of an angular momentum vector. (An algebra, over a field, is a 
system of elements and the operations of addition, multiplication, and scalar 
multiplication. Addition is commutative and associative. Multiplication is 
distributive with respect to addition and, for the angular momentum algebra, 
associative (12)). It is remarkable that the commutation rule suggests and 
permits a demarcation of the elements of the algebra into linearly independent 


(1) L. E. Dickson: Algebras and Their Arithmetics (New York, 1938). 
(2) D. E. LirtLEWwo0D: A University Algebra (London, 1950), chap. XIV. 
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sets of tensor operators. These tensors are related to the forms introduced 
by RACAH and others (**). The algebraic presentation directs attention to 
aspects of these tensors which reveal the structure of the angular momentum 
algebra, and it suggests manipulations. which inevitably lead to an operator 
calculus which effects the immediate computation of observable quantities, 
transition probabilities rather than probability amplitudes. 


2. — Tensor notation and a fundamental theorem. 


The manipulations to come are best handled by tensorial methods. The 
space is three dimensional and the summation convention is adopted: each 
dummy index appearing twice is to be summed over the values 1, 2, and 3. 
A very valuable entity makes it first appearance in the expression of the com- 
mutation rule. This is the completely antisymmetric unit pseudotensor of the 
third rank (°°), written as e,,,. A component of ¢;;, changes sign under the 
interchange of any two indices, and the non-vanishing components are phased 
by setting ¢,,.; equal unity. Two rules obeyed by e¢,;, are easily proved: 


(1) Esjg€oun = 912957 — Oin 95x 9 
19) pa} 
(2) E; uvEuvi > i 205; ? 


where, as usual, 


The commutation rule can be written as 


(3) Ss S;S; — (See S;) = 1€... S 


WO. C's 


Multiplication (and summation) by &,;; yields an alternative form 


(4) ASS) 
or 
(5) ISS, = 


Eq. (3) is similarly derivable from Eq. (5). 


0) A. R. Epmonps: Angular Momentum in Quantum Mechanics (Princeton, 1957). 
O M. E. Rose: Elementary Theory of Angular Momentum (New York, 1957). 

any L. Lanpav and E. Lirsnrrz: The Classical Theory of Fields (Cambridge, Mass. 
LJoO1L). 


6 A E = = 
(°) J. L. Synen and A. ScuILD: Tensor Calculus (Toronto, 1952). 


482 


ON THE ALGEBRA OF ANGULAR MOMENTUM 3 


In the construction of a linear algebra, it is necessary to find a set of basal 
units, a set of elements in which any element of the algebra can be expanded. 
The components of S are a good beginning, for they are linearly independent. 
It is a fact that if 


where a, and 6 are members of the field (*), then 


Op ae Uy 

i =. 
For, on the expression 

a,8, = b 


perform a commutation with S,, getting 


E 


O — 
nio lac =). 


Perform another commutation, with S,: 


O 
sea Sie 


pido“ ovk 


SO = a, 8, = 


py 


(6) 


Contract this last expression: 


Going back to Eq. (6): 


Division by the member of the field is permissable except in the case that 


The proof indicates more than the linear independence of the angular mo- 
mentum components. Actually, in the writing of the expression 


a, Dà =p 
(*) For purposes of this paper, the field need be no more than that of the complex 


numbers. 
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there is the tacit assumption of the existence of a modulus, a unit element, 
in the algebra.’ Henceforth, let the algebra have associated with it a modulus. 
The theorem just proved concerns the linear independence of the three com- 
ponents of S and the modulus. 


3. — Products of components. 


The search for basal units of higher powers of angular momentum com- 
ponents can be carried on in this way. The elementary quadratic form in 
the components of S is S,;S;. This product is not a «pure» or irreducible 
quadratic, inasmuch as 


S,S; = 


bo|H 


; DR TRO LE ate Ss. 1 : 
(SS; + SS) + = (S,S; —S8,S;) = DI (SiS; + S;8;) na Ergo + 


di 


The symmetrized form, S,S;+8,8;, seems to be the thing to consider. In 
fact, for a product of any degree, the fully symmetrized form, only, need be 
considered. The factors (effectively the indices) in a product of n components 
of S can be re-arranged into any of the n! permutations on » symbols by suc- 
cessive commutations. Each commutation gives a re-arranged product of 
degree n, plus a homogeneous polynominal of one degree less. By the addition 
of the n! re-arrangements it is possible to write 


ni S;S;...S, = (sum of all permutations of S,S;....8,)+ 


+ (a linear combination of products of lesser degree), 


The second degree product has been expressed as a sum of the symmetric 
forms of second and first degree. By induction the product of n components 


of S can be expressed as a sum of the symmetric tensors of degree n and lower. 
As notation I shall use 


(7) DI: = | SS; eee Sra ’ 
where the symbol | indicates that the form to the right of the symbol is to 


be brought to iull symmetry in all its non-dummy indices. 
From the definition, Eq. (7), it should be clear that 


(3) Dos = | Do FA | S; Sila ) 


with the convention that Det the modulus. As for the commutator of 
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>” with S, the symmetrization does not hinder the application of the rule 


(AL BO) = (AB) Bt, C) 
so that 


(9) (Sis (mn) )=4 | € (n) 


essre pio rana 
duo 


The indices to be symmetrized are here indicated beneath the double bar for 
further clarity. Eq. (9) can be looked upon as a generalization of the basic 
commutation rule, Eq. (3). The Y®, however, do not admit a generalization 
of the commutation rule expressed as 


CY CY E On CY 
Ei Si) S,) IZ 


for multiplication of both sides of Eq. (9) DY Equi gives 


Ù (n) SL (n) pe ae, (n) (n) 
Eani(Nu ’ ie) x VE vi | E io OF. UE) iE pio ae Eaui | Enio dra) dr 
j...n 


PELLI 


(10) Di i(2 ee Be (n a 1) oo 3A | Ò;a SN) 
Diveste 


Er(Suo De) = aL (n n 1) ae Sei | Osa SAI È 


Jie 


The contracted tensors disturb the harmony of the commutation rules. 
In fact, they prevent the linear independence of the 2” tensors, as an attempt 
at a proof will show. There is a compulsion to modify the + so that a con- 
traction yields a null tensor. To illustrate: 


Zi; = Sat ISPISPRE 


DE = 282. 
Let 

Di; = Zi; — $8055 
Then 

Dig 


D,; should be generalized. The task is to preserve the fully symmetric homo- 
geneous polynomial character, yet assure a vanishing contraction. An attempt 
at a construction is to begin with 
(n) a nel) O 
DE dl | Di; Sa 


7 


Del 
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following Eq. (8) for Yy,,, and then to subtract off 3n(n — 1) diagonal terms: 


I(n— 2) 
Dee JS = | De ns. = | Ò;; bi i 


where the B are to be fully symmetric and also of vanishing contraction. 
Given that the contraction of D‘-» vanishes, the condition that D\” behave 
likewise is 

0 = 2D%"8, —3Be,” — An — 2) BY.” , 


which can be solved to give 


DI 
(n— 2) = (n ù Y 
Br AGLA Sg Dar doh 
2n — 1 


Therefore, let a multipole tensor (*) be defined by 


9 


| (TDI 
Og je Pie 


je _|pevg Se 
ut) y | ey) 
with 
D©=1 
and 
Db ci: 


There is a direct way of modifying a symmetric tensor so that the modifi- 
cation has a vanishing contraction. It is a messy business and is demonstrated 
in Appendix I. The explicit reduction procedure applied to 2° ,, which obeys 
the iterative formula of Eq. (8), must lead to the DY? , which obeys the iter- 
ative formula of Eq. (11). The explicit procedure shows that D7, is a linear 


combination of 2°), and its contractions. A p-times contracted 2%? ,, can 
be man-handled into lesser degree symmetric tensors multiplied by functions 
of S*. Therefore, D” , is a linear combination of the symmetric tensors of 


degree n and lower, the coefficients being functions of S?. The roles can be 
reversed by and inductive argument based on 


ii = D+ 8 S?6,;. 


That is, a symmetric tensor, 2)", can be expressed as a linear combination 
of multipole tensors, but the coefficients are functions of S?. Any product 
of angular momentum components can be expressed as such a linear combi- 
nation of multipole tensors, but the coeffcients are not simply numbers. The 
multipole tensors are not quite basal units at the moment. 


(*) L. Lanpav and E. Lirsurrz (5), p. 104. 
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4. — Properties of the multipole tensors. 


DIRÀ ode: i 
D® is a homogeneous polynomial in the components of S. It is as true 


ij. 


5 es È ì È : È 
for Dj, a8 it is for v® that the operation of commutation with S, is 


n in 
(mn) 


distributive over the indices, each index in turn being replaced by de, Dio 
The only point to be checked is the commutation with factors like S?0;;. 


But this is all right, as one can write 


(So Std) 8° ExioSia) = 182(8,;; =e Ei) ‘ 


uio 


which is equal to zero, as it should be. Therefore, one form of the commu- 
tation rule is (*) 


(12) (RD 


19... 


)=ile n) 


pio Gin * 
dis 


The multipole tensors were introduced to permit the transformation of 
Eq. (10) into 


(13) EAST Da) = un sl (VEE 


the second form of the commutation rule. 
Take Eq. (12) and bring it to full symmetry in all the indices: 


| (8 


n? 


Dia) = iui 


pio Ginn si 


Symmetrization of the antisymmetric unit tensors will cause cancellation, 


so that 
(14) | (So Dee) = 0. 


From this equation, coupled with the condition of vanishing contraction, it 
follows that 


Cy (n eS. 
(15) CASO 


Eqs. (14) and (15) allow flexibility in the generating formula for the multi- 
pole tensors, Eq. (11), in that the multiplication can be performed from either 


side. 


(*) This equation could be derived by induction on Eq. (11). 


= 
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To get an identity of tantalizing appearance, replace DOS in the left- 
hand side of Eq. (13) by its formula given by Eq. (13): 


Esuy Evo (Syr (Sa; Di = M+ DEL po 

(De. BS) ) + (8, (So Did E DIO 
The last line results from the use of the identity 

(A, (B, 0)) + (B, (0, 4)) + (0, (A, B)) = 

and the application of Eq. (15). Going on: 

cid) DDD 
and finally 
(16) (Si (Sy De a(t ra) De 


LI Die ig) bait ® 


Because D’? , commutes with S?, the last equation can be written in two 


CICCO 


other ways: 


(17) SD 2S dan), 
and 
(18) 8,(8,,5 DI) = }n(n+1)D ,. 


A clue to the most important information contained in the multipole ten- 
sors is in the calculation: 


8, D,,= 8,(8,5, + 8,8,—876,,) = 


pui 


= 828,4 Si, SS, —$ 878, = £878, (S?—1)S, 


because of Eq. (17). In final form 


It can be expected that the tensor defined by 8 DI be proportional 
to Di”, because that tensor must be a fully 11 homogeneous poly- 
nomial in the components of S, and it must have a vanishing contraction. 
A proof by induction will confirm this expectation and will reveal the signi- 
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© 


ficant proportionality factor. It turns out that 


\ (n) —_ n° 2 1 2 n— 
(19) Da Dia = DR = i bi ce 7a (n° = 1) Daci E; 


Set n=2S +1 to emphasize the bracketed factor as 
[S°-S(S+1)]. 
The recursion relation for the multipole tensors can now be re-written as 


| 
(20) DIS È | SD 


2(n —1)? Weds 
PSA SE se yp LD 
(2n —1)(2n — 3) 


Assurance should be given that the identities derived here do permit a proof 
of the linear independence of the multipole tensor components. A lengthy 
proof appears in Appendix II The statement is that if 


(n) @) ff 
DI Op wifes =0, 


(n) 


where the scalar tensor a,” is fully symmetric in its indices and has vanishing 
contraction, then 
at, = 0. 

Of course, any tensor undergoing inner multiplication with a multipole tensor 
can be brought to full symmetry and modified so as to have vanishing con- 
traction, without adding anything to the product. This statement of the linear 
independence of the multipole tensor components means that all the linear 
dependence is recognized when account is taken of the symmetry of iD 
and of its vanishing contraction. It is therefore possible to count the number 
of independent elements of a multipole tensor, and this is done in Appendix ITT. 


The number of independent elements of TDS Dy SDS 


5. — Spherical tensors. 


The multipole tensors which were introduced here, in the search for an 
algebra, with no thought of the eigenvalue-eigenvector problem of angular 
momentum, assume much more significance when they are related to certain 
tensors which may be typified as spherical representations of the cartesian 
forms developed here (3), (4). These spherical tensors are gene ralizations of 


(er 
A 
_ 


10 A. MECKLER 


the well-known shift operators of angular momentum. These shift operators 


can be defined in the following invariant way: 
Let a be a real, three dimensional vector normalized so that 
i 


(21) @ a, = IR 
Let 

Do = a, Si, 
Define 

i 


and insist that 


DI ra) va | 7g 
(22) [DES 8 Gre 
The expansion of the last equation is 


14,07 8,84 = 4, 8, - 


VS ie U o 
Coefficients of S, can be equated. Therefore 
(23) E I MOL c 


Cp» | Lal Gg 


Introduce the antisymmetrical tensor dual to a : 


doi = E 4,3 — dg 
The converse relation 
=a 
a, REI Ena ve 


is also true. 
Eq. (23) can be re-written as 


(24) a. a 24a" 


This equation can be interpreted, properly, as an eigenvalue-eigenvector equa- 
tion, for the normalization condition, Eq. (21), implies that 


es! 3 21 
1= $8, 44 Erol, = 4(4,,4,, — d4,,) =—4Trd?. 


That is 
Tr d—=—2., 
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One eigenvalue of a real, antisymmetrical tensor of the third order must be 
zero; the other two must be conjugate imaginaries. The condition on the 
trace of d? means that the eigenvalues of d are 0, +i, and — 7. By Eq. (24), 
at is the eigenvector belonging to the eigenvalue +7. The eigenvector be- 
longing to the eigenvalue 0 is a,, for 


dig, = Ear ta = 0 


because of the antisymmetrical nature of ¢,,,,- 
If a, denotes the eigenvector belonging to the eigenvalues — 7, then 


satisfies the commutation rule 


(25) (You, eee 


(26) aia, =1. 


With this normalization it can be proved that 


(27) call ia, 
so that 
(28) (ee enor 


Now ata, + a, a; + 4,4, is a tensor which acts as the unit tensor on the 
ee pv pv 
j i ti si Nn (fo APP V > Pg TIE 
three linearly independent vectors 4,5 4, > and a, - These three vectors span 
the space, and so 
+ = sant SNA 

(29) a; 4, + a, 4, oe Oi da 7 
acter of the transformation.) By 


(This is a statement of the unitary char 


re 
oa 
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multiplication 


(30) S,= a,a*S+ ag a-S+ a, a'-S 
and 
(31) S°= (a:S)?+ at -Sa--S+a--Sat-S. 


Eqs. (28) and (31) can be combined to give 


(32) at-Sa--S = 3[S?—a:S(a-S —1)], 
(33) a--Sat-S=1[8° — a:-S(a-S+1)]. 


The spherical tensors are generated in this way: 
Define 


wm =(y{y 
pete Nk a) 
RD 


and so on down to Y,". (The normalization is terrible, but that’s not im- 
portant here.) A general commutation rule, easily proved, is 


(34) ue 


Eq. (34) leads to the interpretation of Y¥ as a shift operator which preserves 
S? but shifts the eigenvalue of a-S by M units. An angular momentum 
system whose highest possible multiplicity is 25-++1 can be characterized by 
the statement 
9 FPSS 
(35) Vrs I 0. 
Now, 
n A 
BRIO SER e A II 


tj... 


because 4, ... ay is a Symmetrical tensor of vanishing contraction. If 


then 
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The linear independence of the multipole tensor components would demand that 


put this is absurd. Y*=0 must mean that DI’ ,=0. The converse is ob- 
viously true. 


6. — Structure of the finite algebra. 


What are the consequences of constraining the angular momentum algebra? 
Let it be that 


(36) eq 
but that no lower tensor vanishes. Multiply Eq. (36) by Si: 
AU 
which by Eq. (19), becomes 
(Ch a ae 
pe #0, so it is permissable to perform another multiplication (by S;) and 


so on, arriving at 


(37) 


ASBI fee 
(38) tl |s*— =: |=0. 


Bq. (38) can be called the Hamilton-Cayley equation of S*. It is well 
known what to do with such an equation (°). Let 


N — k)? — 
(39) a ee 
and let 
ay, Pe SOND) 
Ehi i aie) 
ae VA 
zk 
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By Eq. (38) 
(41) S? E(k) = E(k)S? = a(k) E(k) . 


E(k) is a polynomial in S? of degree (N—1). > E(k) is a polynomial in S?, 
} 


defined over the field, of degree (N — 1). At the N points, (k=0,..., N—1), 
> E(k) has the value 1. This is a statement concerning the field coefficients 
k 


in the polynomial; it has nothing to do with the nature of angular momentum. 
Therefore, 


(42) Hey = 1. 


The summation of Eq. (41) yields 
(43) SS (kE) 
It follows from Eq. (40) and (41) that 
(44) E(k)E(k')= 0, kak, 
and that 
(45) Fi (ie) == Ek): 


Eqs. (44) and (45) characterize the £°s as orthogonal idempotents. 
Any component of S can be resolved in this fashion, for 


S:= > E(k)S, 


or 

(46) SI 
where 

(47) h;(k) = E(k)S; . 


The interesting thing is that 4,(k) is itself an angular momentum: 


(hi, h;) = (ES, BS,) = E(S,, 8,) = ie, ES, = de, h 


io a ° 


The above steps are possible because E commutes with SE 
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From Eq. (47) 
h? = ES;ES; == ES? 
(48) 
Moreover, 


(49) h;(k)h;(k')=0, kez ke 
h has a very simple Hamilton-Cayley equation: 
(50) heh; = ah, . 


The A is an angular momentum. Its components can generate a sub- 
algebra in the algebra of S. Let a multipole tensor formed with the com- 
ponents of A be denoted as 

(n) a 
Di La(k)] . 
Then 
TRIS 


i). 


[h(k)] = E(k)Dg.[S] = 0 


for k > N n, because E(k) would then contain at least all the factors ne- 
cessary to satisfy 


N-n=1 (N — p)? — 


ci (DETS]F:08 


Therefore, 


DI P[hk)]=0. 


To summarize, the important equations of this section can be written as 


(51) LS 

(92) S—Y Mh+1)E(h), 
(53) CASI TP 

(54) SIE 

n°) Die tele 9 5 


Uni 
D 
Su 
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where the range of h is 


in half-integral steps. Note also that no A exists for h=0. 


7. — Structure of the sub-algebra. 


As stated in the last section, the component of fh can generate a sub-algebra. 


In this case, because 
h?h; = h(h+1)h; 


the multipole tensors do form a basis for the sub-algebra. The modulus is 
E(h), as demonstrated by 


pre 
dii 


The sub-algebra is finite; the (2%-+1)st multipole tensor vanishes. The order 
of the sub-algebra is easily computed as 


Speed) = 441: 


In terms of the algebra of S, the equation 


> Bh) = 1 


indicates a decomposition of the modulus into a sum of orthogonal idem- 
potents. In the theory of algebras, the question arises as to whether an idem- 
potent is primitive. An idempotent, e, is not primitive if there exists another 
idempotent, w, such that 


Ch neu. 
A non-primitive idempotent can be written as the sum of two orthogonal idem- 
potents, for 


e=u+(e—u), 


where 


u) = e —Quetu —e—-u, 


ule—u)=(e—u)u=ue—u?= 0. 
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If E(h) is not primitive, the modulus of the S algebra can be further de- 
composed. In terms of the sub-algebra of A, the question of the primitiveness 
of E(h) can be phrased as a question as to the decomposition of the modulus. 

The sub-algebra is characterized by 


(56) TM 


Chose now an a,, a real unit vector such as was used in the discussion of the 
spherical tensors. Multiply Eq. (56) by 4,4,..: 


(57) ADESSO 
Let 
(53) A = a,4; oe @, Die, ) 


where n is any integer. Clearly, A™ is a product of E and a polynomial of 
degree n in a-h. For example, 


AV= a-h=Ea-h, 


A® = q;0;(h;h;+ hjh; —8Mh+1)Eòd), 
= 2(a-hP_8KA+1)F, 
= [2(a-h)—3hh+1)]E. 
Eq. (57) is equivalent to 
(59) AGrD =). 


Here again, as in the last section, we have a Hamilton-Cayley equation. It 
is possible to factor Ed. (59) so as to read 


where the 7, are the as yet unknown roots of the Hamilton-Cayley equation. 
Formally, it is permissible to introduce 


2h ahi 
e 


p=0 


Fm 
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and it would follow that 
Y on(m) = E(h) 


m 


a’ h = > Ain On(M) 


(60) 
on(m)en(m')= 0, m~Am 
0n(M) O,(m) = On(m) 
and 
a-ho,(m) = on(m)a-h = Aon On(m) 
(61) 


h?0,(m) = h(h-+ 1)0,(m) . 
To determine the /,,, bring in the shift operators: 
a-ha--ho,(m) = a~-ha-ho,(m) — a~-ho,(m) = (An—1)a~he,(m) - 
Multiply on the left by @,(n): 
(62) o,(n)a:-ha--ho,(m) = (An —1)o.(n)a~-he,(m) . 
But a-h can also act to the left, so that the last equation is equal to 
(63) Zn0(n)a-h 


On(M) . 


Comparison of Eqs. (62) and (63) leads to 


(64) on(n)a~-he,(m) = 0 


unless 
ig An 
Eq. (64) can be summed over all values of n. Because 
di: On(m) = E(h) 
it follows that 
a--ho,(m) = 0 


unless there exists a A=/,,—1. 
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There are only 2A+1 roots. Let A be the minimum root. Then 


a--hmo,(4)= 0. 
Multiply by a*-h: 
atha~-ho,(A) = 0. 


Eq. (32) can now be used to give 


[h? — a-h(a-h —1)]o,(4) = 0 
which requires 
[h(h +1) —A(A—1)Jon(A) = 0 
and therefore 


h(h+1)—A(A—1) =0. 


The two solutions of the equation are 


and, since Z is the minimum root (h is positive), 
(65) A=—h. 


It is impossible that 
a*+ho,(4) = 0 
for this would imply 
a~-ha*-ho,(a) = 0 
or, by Eq. (32) 
[h?—a-h(a-h + 1)]@,(4) = 0 
or 
h(h+ 1) —A(+41) = 0 


which is not possible unless A= h= 0, the case for which no h exists. 


fore 
a‘ -ho,(A) #0. 


It can be proved by the methods already applied that 


at-ho,(m) = 0 
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unless there exists a Z such that 


A=Amt1. 


In particular, there must exist a 4= A+1. In this way, all the roots are 
found to be the integers, or half-integers, going from 


(66) I) Wine 0s 


At this point, the modulus of the S algebra is expressible as a sum of idem- 
potents, H(h), which characterize S?, and these idempotents can be further 
resolved into idempotents, 0,(m), which characterize a-S. 


8. — The super algebra. 


The product of any two elements of an algebra is an element of the algebra. 
In particular, the product of any two basal units can be expressed as a linear 
combination of the basal units. The scalar coefficients can be arranged into 
a multiplication table. In this section, it will be indicated how such a mul- 
tiplication table could be constructed for an angular momentum algebra. 
However, it seems that for a large class of problems, it is not necessary to 
know each entry in the table. Only certain entries, or certain combinations 
of entries, are of physical significance. In any case, the constants of multi- 
plication (*) can be specified when it is recognized that to every angular mo- 
mentum algebra there is associated another algebra, a master or super algebra. 
The super algebra is formed in this way: denote the algebra in which the 
angular momentum is imbedded as A. Consider the product of any two ele- 
ments in A, wy. Represent the product as the result of an operation on y by « 
(from the left): 
LY > GY 
or, equally well: 
LY = LY op 
(from the right). This set of operations, abstracted from the basic algebra, 
is itself an algebra, the super algebra. It is not a trivial algebra, for consider 


the operation of commutation with S,, 


(67) Ga (Seana 


(elie. DICKSON (*),. Daal 
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g; is an angular momentum, for 


(Gi, 91)e = (Si, (S;, x) # (S;, (Si, w)) 
= (8, (85, @)) + (Se, (#, Ss) + (0, (Ss, 82) 
= (a, (S;, S,)) 


Lie, (85 2) 


or 

(68) (9; 93) = Eti0Io 
Note that 

(69) (PAN ie Ne 


Beyond its character as an angular momentum, g; is distinguished as a 
distributive operation (*): 


(70) gly) = (9:0)Y + (9,4) - 


Moreover, g; can be written as the sum of two commuting angular momenta. 
Define the operations, m, and n;, by 


(71) me = Sa 

(72) na =— td; . 

Each of these operations is an angular momentum: 
(73) (m;, Mj) = 1€;;,™, 

(74) (Ni, M5) = deo 

and they commute: 

(75) (m;, n) = 0. 

Egs. (73), (74), and (75) are easily verified. Clearly, 


(*) Sometimes called a derivation. 


a 


A. MECKLER 


bo 
bo 


Tf the constraint on S is 


(25 +1) Pr 
DI [S]=0, 
the constraints on m and n are 


(77) eWay) SAD 


mo S+ 

(78) Dn] = 9. 

In Appendix IV, it is proved that Eqs. (76), (77) and (78) lead to 
(79) Di "gl = 0. 


Here is a set of angular momenta g, m, n, to which the methods of this 
paper, and all other papers on angular momentum, can be applied. Even the 
vector space on which they operate is given; it is the algebra A. The modulus 
of the super algebra can be reduced into a sum of orthogonal, angular mo- 
mentum idempotents, each of which is a product of commuting idempotents: 
a (g?, a-g) idempotent, an m? idempotent, and an n? idempotent. Each ele- 
ment of A can then be written, identically, as a sum of projections onto simul- 
taneous eigenvectors of g?, a-g, m?, n?. Already certain eigenvectors have 
been identified. The spherical tensor 

VOR deo DOC A] 


L dee 


satisfies the eigenvalue equations: 
Gays e N(N+ na 
by Eq. (16), and 


5 E na i Y 
a:-gY," = 4, ...a,4,(8 


(N ag = (N) 
; yw Diy) = ia, 074, || DE 


— Na, ax. Oy GEG Dat co Maga è DOZZINA 
and 
mY <* Th] = hh+1)Y5"[h], 


n° Yo" [h]=Mh+1)Y;"[h]. 


Y,***, which belongs to the eigenvalue — N+K of a- g, is generated by K 
operations with a*-g, the shift-up operator for a-g. 

The multiplication table of A could now be constructed by application of 
the combinatorial methods of angular momentum (#4). To amplify this state- 
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ment, the case of an intrinsic angular momentum will be considered and only 
those constants of multiplication will be derived which are of physical signi- 
ficance. 


9. — Calculus of an intrinsic angular momentum. 


In this Section, S? is a constant, S(S+1) times the unit element. Without 
going into the details of how to assure the equivalence (+), let the algebra of S 
be taken as equivalent to a matrix algebra of order (25+1) (?). S; is hermitian. 

The primitive idempotents, 0g(m), are the projection operators for a-S. 
In a representation in which a-S is a diagonal matrix, og (m) has unity as its 
m-th diagonal element and zeros elsewhere. In a representation in which a-S 
is not diagonal, a diagonal element of og(m) is to be interpreted as the proba- 
bility that the system be in the a-S= m state. The expectation vaiue of any 
operators, X, in a state characterized by a definite value, m, of a:S is, (7): 


Tros(m)X . 


Therefore, only that part of the multiplication table which enters into the 
evaluation of traces of products need be investigated. 

The algebra of S is the vector space of its super algebra. A metric can 
be introduced: on the matrix algebra generated by S define the inner product 
of x and y as 

Traty = e], 


where «* is the hermitian conjugate of «. This trace operation satisfies all 

the conditions of an inner product (*). With this definition of the inner pro- 

duct, it is easily proved that the linear operator’s g, m, and n are all hermitian. 
In a previous Section, the operator 


(n) 
i). 


(n) — 0; x Un D 


was introduced. For the case of a spin of multiplicity 25 +1, there are 25+1 
such operators (0 <n < 25). Bach is a polynomial in a-S. The set OLA 
ig a set of commuting operators, a set of linearly independent commuting 
operators. The A can be simultaneously diagonalized and, in diagonal form, 
the array of eigenvalues of an A™, to be denoted as A™(m), can be considered 


) P, A. M. Dirac: The Principles of Quantum Mechanics (Oxford, 1947). 
(8) P. HaLmos: Finite Dimensional Vector Spaces (Princeton, 1948), p. 126. 
) M. Born: Atomic Physics (New York, 1946), p. 222. 
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a 25-+1 dimensional vector. The vectors, A™(m), span the space of 25+1 
dimensional vectors. In fact, the set of A(m) is an orthogonal set. That is, 


(80) DAMA (m) = 0, n#A N. 


m 


For, 
YAM(M)A" (n) — Tr AM AM = CAM) | Am?» 


m 


and A™ and A are each eigenvectors of g?, but each belongs to a different 
eigenvalue: 


GRA = n(n + 1)A™ ci 


Eigenvectors of a hermitian operator, belonging to different eigenvalues, are 
orthogonal. 

There is another way of proving this orthogonality relationship. The trace 
of a commutator is zero. Therefore, the trace of any multipole tensor, except 
the modulus, is zero: 


(81) Tr D'ORO ne 0: 
because 
n(n + 1)Di_ = (8,, (8 


nu) 


(n) 
Dj.) - 
What of the product of two multipole tensors? D7, is an eigenvector of g? 
belonging to the eigenvalue n(n+1); Di’, is an eigenvector of g? belonging 
to the eigenvalue n’(n’+-1). Remember that g;, an angular momentum, has 
a distributive property: 


gry) = (gir)y + @(Giy) . 


Because of this special property of g;, it is fruitful to think of the product 
of two multipole tensors as an outer product, an element of the direct pro- 
duct of the vector space with itself. The methods of angular momentum syn- 
thesis can be taken over. The product 


DO 


tj... 


(n') 
D5) 
presents a problem in the coupling of two angular momenta. If an angular 
momentum of g?= n(n + 1) is coupled to an angular momentum of CO fa 


=n'(n'+1), then the resultant angular momentum can take on the values 


7 PISTA ri Pr / 
nin’, meni 1s ch; ijn —n’ 
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Therefore, in the product, D' 


(n') : o 
lip there will appear only those multi: 
pole of orders 


i'j'...n' 9 


ntn', wen a wy In-n'|. 
a , 
Tf nn’, the modulus does not appear, and so 


(82) TD” D® .—0, eae 


in 8H N 


Consequently, 
Trp A AGIO NEN. 


The A™ are orthogonal. How are they normalized? The norm of A‘ is 


Tr [AP 2. VAM(M)A®(Mm) 2 


A™(m) is a polynomial in m. A recursion relation can be derived by mul- 
tiplying Eq. (20), the recursion relation for DI ,, DY 44; ... An: 


ij...n? 


2(n—1)? SIs ia 21) ao 
= ? 


Aw = na: SAS PR 
(2n —1)(2n — 3) 4 2 


or 


(83) Am) = nmA%(m) ie ian 33 CS +m) 25 — n+2) A0-D(m) . 


The equation indicates that the A(m) are polynomials of definite and alter- 
nating parity. 
Let 


(84) A(m) = JnPn(™) 5 


where p,(m) is to be an orthogonal polynomial normalized to unity: 


(85) Y prlM) Pulm) = Onn 


Then the norm of A™(m), 


(86) YAM(m)A(m) =. 


m 


There is a formula, called the Christoffel-Darboux formula, which holds for 
any orthonormalized polynomial (19) which, when applied to a polynomial 


(19) G. SZEGÒ: Orthogonal Polynomials (New York, 1939). 
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of alternating parity such as p,(m), reads as 


(87) Pa(M) cui F,MPna ae pies : 
where 


IG = LA RE È 


Substitution of Eq. (84) into Eq. (83) gives a recursion relation for p,. The 
Christoffel-Darboux formula is satisfied if 


wae 
(88) P= a 2 i 
and 
ng N\4 Da 
n “n [en iP Qn 4 È Za isi 
where 
n DS 
and 


(2n —1)!! = (2n —1)(2n — 3)(2n — 5)... 


For purposes of numerical computation, a generating function for p, can 
be derived. This is done in Appendix V. The result is 


(2n + 1) (p —n —1)!}8 q—p 
91 n == ,! — — — — —————— n 
(21) BAe Nie (p+ n)! a n n ; 
with 
m=—S+q, 0<¢<25 


and A is the difference operator: 


Af (a) = fla +1) — fla) 


È a RI, 
n} ni(g— n)! 


It is marvelous that these polynomials are already known and are associated 
with the name of the ubiquitous TCHEBICHEF (eds 


and 
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bo 
+1 


The initial definition of A°® was 


1 = (I 

A ad An Dip 

Tt is convenient to modify the scalar tensor, 4,4; ...A,, SO as to have a va- 

nishing contraction. A set of scalar multipole tensors can be constructed out 
of a; by the recursion relation 


(n) dl i Î 


3 a = — 
tj... n | 


9 | 
(n-1) ae || (n-1) 
Ad; 4j... “a Dn e 1. | MyWire... o 5 
a J | 


with a =1, a?=a,. The scalar multipole tensors are to be fully symmetric 
in the indices and are to have vanishing contraction. By the methods used 
in the case of the angular momentum tensors, it can be proved that 


(92) Qui. = 4 On 9 
so that the equation for 4;;., reduces to 
(93) pet a ina, 
ie ee Asch sealar multipole tensor formed from the unit vector d;, then 
(94) aS Dien = se P,(a-b), 


(on —1)!1 
where P,(x) is the ordinary Legendre polynomial normalized so that 
P.(1)=1. 


Eq. (94) is derived in Appendix VI. 
There will be a requirement, further on, to compute 


n ) 
Ty A BO ; 


where 
A = a DO 


ij.in 7 dn I 


B@ = bh DO 


iano Ae 


The computation begins with 


(95) Tr A9BM = a be _ Tr DO D3 n 


tim Ween 
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Now 


IDO 


‘I 
+ terms which yield zero when multiplied by multipole tensor coefficients. 


(n) (n) 


Multiply by 4@,, 4; : 
(n 1)? 


OF ES: 


Tr[Am}? = A,nla? ag. = 


so that 

(2% —1)!! , 
(96) An = (nije = Dis . 
Therefore, returning to Eq. (95), 
(97) Tr A®B® — g-P.(a-b). 


For the computation of an expectation value, the factored forms of 05 (M), 


S aS: —k 
ESE m — k | 


Fm 


os(m) = 


is not best manipulated. 05(m) is an operator which commutes with a-S and 
the A™. As a polynomial in a-S, 05(m) must be expandable into a sum 
of A™: 


(98) os(m) = > 0,A™. 


To evaluate the coefficients, multiply Eq. (98) by A™, and take the trace: 
TIA Oe (Gy) aC er Ate 
But TrA"os(m) is the diagonal element of A™: 


A™(m) — Tr A™ O05 (m) = (E Tr pace 
so that 


In Pn(m) = Cn dì 
or 


Pn(m) 
Gn 


The expression for the projection operator is 


(be — 


(99) MA pani AM, 


3.8 
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10. — An illustrative problem. 


Given a spin definitely quantized along a unit vector a with component m, 
what is the probability that it is quantized along a unit vector b with com- 
ponent m’? The answer is given, in effect, in the texts (?), (4). From the 
definitions given there, it follows that the probability is the square of the 
matrix element of the rotation matrix: 


lE az i do; BB) | b) 


where 8 is the angle between the vectors a and b. d®, (8) is computed as 


mm 


a5,,,(B) = [(S + m)!(S— m)!(S+m)!(5-m)!} 
DI (a1)! 
r (5 —m'—r)! (Stm—nir+m—m)irl 


b S4m-m'-2r B m-m+27 
‘| COS 5) — Sn 9 A 


The calculus of the preceeding Section effects a straight-forward answer 
to the problem and a more compact expression: 


aa == Lr o(m)o'(m') , 


with 
o(m) — Soe AM) 
In 
j n m' Jin) 
(mi) = SP! peo 
0) 
By the formulae of the preceding Section: 
(100) Pam = > Pa(M)Pa(m')P.(a' db) . 


41. The Tchebichef polynomials and Clebsch-Gordan coefficients. 


There is a startling identification of the Tchebichef polynomials, Pn (S; m) 
with certain Clebsch-Gordan coefficients. It is a fact that 


(101) DnlS;3 m) = (— 1)° "0(5, S, n; mi —M), 


(ec) 
(=) 
Usi 
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where the O-notation is that of ROSE (4). To establish this identity, examine 
the relationships of the eigenvectors of the super algebra: 
Because 


a-mos(m)= mos(m) 
an os(m)=—mos(m) 
mo5(m) = n°os5(m) = S(S+1)og(m) 
Tr os(mos(m)= <os(m)|os(m)) = 1 


the projection operator is, except for a phase factor, a normalized, simultaneous 
eigenvector of m?, n?, a-m, and a-n. 


Because 
AM 
ie fa ates 

Am AM 

9? ny n(n + ra 4 

AM AM Am 
Mm? = i? = S(S+ 1) ‘ 

In In n 


Am Alm | Am 
Tr |—| = = yee Da, 
In In | In 
A[g" is, except for a phase factor, a normalized eigenvector of g?, a-g, m?, 
and n?. The equations 
AM 


ae 
n 


eS(m) = > Pa(S; m) 


— = Y p,(S; m) 0S(m), 


are expressions of the transformation from the uncoupled representation to 
the coupled representation and inverse transformation. The Clebsch-Gordan 
coefficients are defined as the elements of the unitary transformation. Pn(S, m) 
is, therefore, proportional to O(S, S, n; m, — m). Comparison of formulae gi- 
ves the phase factor in Eq. (101). 
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APPENDIX I 


Reduction of a symmetric tensor. 


Given a tensor, t;;..,, Symmetric in all its indices, it is required to work 
upon it so as to assure a vanishing contraction. 

Consider a term of the following sort: select p pairs of indices and con- 
tract t,;... over each pair, ending up with a p-times contracted tensor such 
BS Luuvr...mn- (There are p dummy indices). Now distribute the selected 2p 
indices over p Kronecker deltas. Form the product of the Kronecker deltas 
Bnd ty... Repeat the process so as to bring the form to full symmetry. 
Call this symmetrized form, made of p-times contracted tensors and Kro- 
necker deltas, Ti}. .n- 

The number of terms in 7%. is 


n! il 
(n —2p)! p! 2°’ 


for, n!/(n — 2p)! is the number of ways of leaving behind (n — 2p) equivalent 
indices. 

1/p! is to take care of the permutation of the Kronecker symbols. 

1/2” accounts for the symmetry of each Kronecker delta (*). 

The modified tensor shall be 


(k) 
bai + alora Sa 200 + Cy iin SF cose 


Its contraction shall vanish. Under the contraction of i and j, let's say, Dee 
splits into two types of terms. There were the terms in which the pair, (è, 9), 
appeared in the t,,;;.. factor. These (n — 2) !/[(n —2—2p)!p! 2”) terms end 
up containing (p +1)-times contracted tensors. The second type of terms 
in the contracted T_,, are 


n! (n— 2)! (n £ 


DIRI 
(n—2p)tp!2” M—2—2p)!p!2? (n_2p)! 


in number, and these second type terms retain a p-times contracted tensor. 
The second type terms in the contracted 7” duplicate the first type terms 
in the contracted 7’?-». There is over duplication, since the ratio of the 


(*) n1/((m— 2p)! p! 2%) is an integer, since 


ni n! { (2p ! (2p — 2)! | 


) 
(n—=2p)!p!2? (n— 2p)! (2p)! (2p — 2)! 2! (2p—4)! 2! ag 


a product of integers which suggests another combinatorial synthesis. 


n 
dA 
Uci 


32 A. MECKLER 


number of 7 second type terms to the number of T°?” first type terms is 


(n — 2)! (Qn — - 2p —1) / (n — 2)! ee o rg 


(@—2p)!(p— )I29/n—2—2p +2)! —1) 12 


It must be remembered when i and j appeared in different Kronecker deltas 
in 7, the contraction gave only one contribution. When è and j appeared 
in the same Kronecker delta, the contraction gave three contributions. The 
over duplication is therefore 


(2n — 2p —1) +2 =2n—2p+4+1. 


To assure a vanishing contraction chose 


Cole = —[2(n —p) +1], 


or 


APPENDIX ell 


The linear independence of multipole tensor components. 


Given that 


(n) (n) 
» re e = 0 ’ 


(n) 


where af, is a field quantity, symmetric in all its indices and of vanishing 


contraction, to prove that 


(n) 
a ,=0. 


To begin with, it is no restriction to have a , fully symmetric and of 
vanishing contraction. Care must be taken in the application of this theorem 
to alter the coefficients so that the conditions are satisfied. 

Repeated application of Eq. (16), 


Ù 3 (n) 
(Sur (Su, DI.) = n(n +1)DP, 
will give rise to an unending set of simultaneous equations of the form 


SS n (n = 1)Pa CO, EDGE = (a 


(n) 
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There is no linear dependence of the equations since the coefficients, n(n +1)", 
in the k-th equation, being polynomials of 2% degree in n, are not dependent 
on the polynomial coefficients of the other equations. Therefore, it must be that 


(AII.1) OF Dim = 0 
The idea is to somehow get from (AII.1) to 
AII.2) of Di. = 95 
which is a simple product with no summation involved. To get from (AII.1) 
to (AII.2) one must, so to speak, release the indices. The inner product form 


in (AII.1) can be said to have zero indices released; the form in (AII.2) can 
be said to have n indices released. An intermediate form such as 


(n) (n) 
Apivte..n-D ijt... b) 


can be said to have two indices released. Given that the form with p indices 
released is null, does it follow that the form with p-+1 indices released is 
null? It does. 

Let 


) (n) 
(AII.3) OU mn... Dijon ee 0 , 


where p indices are released. Commute the equation with Sq: 


(n) (n) = 
Aup.. mV | E115Doi...m..m =. 
ij 


To the right of the double bar there will be a number of tensors in which 
indices are released. Therefore, by (AII.3), the equation can be reduced to 


(n) (n) a 
Wray... | cime Daan ME 


M...n 


Because of the symmetry of Osea) this last equation can be reduced further to 


(n) (n) È 
Cpe oun masa etiam SS 0 Q 


There’s no need to encumber the typography with the released indices. The 
last equation can be written as 


(n) (n) pe 
(AII.4) Emo Wini.n Da, min "e 0 b) 


where the single index to the left of the comma represents the set of released 


= 3 - Supplemento al Nuovo Cimento. 
Lie 
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indices. The necessary and sufficient condition that (AII.4) be true is 
(AII.5) nb aac 
Commute both sides of (AII.5) with S,, and multiply by &me: 
CILE) Di...) = Oe aie (Ss Dae b) 
(n "o ILA oi 055 GIL | tim Da; yin + als aerial Craig Os ace ca 


) (n) (n) (n) 
= = a, ei o) VIAN n On [n 1)] DAI Bs Ee 
) ( ) ) ) 
ag [n ae (p cin 1 Vasa Dona | Oni. na SCI —| A Dee Ò,z . 


x 


The second term on the right hand side of the equation can be dropped because 
of (AII.3). The third term can be dropped because a", has vanishing con- 
traction. The sixth term can be written, because of the D®,., symmetry, as 


(n) (n) 
| esa Oo dea Ove ’ 


which vanishes because of (AIT.3). All that is left is 


EI TIRATI, ea dei 


LC,v0.. 


(n) (n) 
+ | Art. Dea . 


The last term on the right hand side of this equation can be re-arranged by 
virtue of (AIT.5): 


(n) (n) <— (n) (n) (2) (n) ( 
een Dy nese age De = Osta Die, PARADE. 


Application of (AII.5) to the second term of (AII.6) allows a final combi- 
nation into 


[2n — (2p + 1)]agm.. Din. = 0. 
It is impossible that 2n — (2p +1)=0. Therefore, 
a™ DO PE 0 
VM TV aa ’ 


which is the statement that the form with P+1 released indices vanishes. 
There is a route from 


(n) (n) 
O DAI ; 


al, n) 
SISI (= 0 . 
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a. is a member of the field. Unless 

(n) 
x LIA = , 
it must be that 


ATDDIINIDI INNI 
The number of independent multipole tensor components. 


Let #,;., be any tensor symmetric in all its indices. There need be no 
connection with angular momentum. The number of indices is n; the range 
of each index is x, not necessarily equal to three. To count the number of 
non-redundant elements in t;;.,, list the elements in dictionary order. That 
is, order the elements t,;.., 80 that i<j<k<...<n. This is an exhaustive, 
non-redundant procedure. 

It is possible to partition the dictionary into disjoint classes whose sum 
is the dictionary. Two elements of the dictionary are termed equivalent if 
they have the same number of indices set equal to unity, the lowest value 
in the range. This equivalence relation is reflexive, symmetrical, and tran- 
sitive. The class (p) is composed of all those elements in which p and only 
p indices are equal to unity. Every element of the dictionary must be in one 
class and one class only. 

Denote the total number of elements in the dictionary by f(x, n). The 
class (p) has }(x — 1, n— p) elements, for it reads as a sub-dictionary based 
jon n— p indices each with a range of a — 1. Therefore, 

n 
f(a, n) = > fla—1, n—>p). 


p=9 


| By a change of summation index, the equation becomes 
(ATII.1) (an) > He 1p), 
p=0 


from which a recursion relation is deriveable as 
(ATIT.2) f(x, n) = f(w,n—1) + f(e@—1,%). 


It is necessary to set 

f(v, 0) =1. 

All scalars are linearly dependent. With this boundary condition (AIIT.1) 
is deriveable from (AIII.2). The solution is 


nm+a-1)! 
COR RENI 


This formula is the Bose-Einstein combinatorial formula. Let » be the 
number of particles in a particular one-particle energy state, # be the degen- 
eracy of that state. Then f(a, n) is the number of distinguishable arrangements 
of particles among the degenerate levels belonging to that energy state (°). 


Uci 
Uan) 
te] 


; | 
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To get the number of independent elements of D®.,, account must, be 
taken of the dependence brought in by the conditions of vanishing contraction. 
There are 3n(n — 1) such conditions. The number of independent elements 
of DS, is, therefore, 


f(3,n) —inl(n-1)=2n+41. 


INDI DIExe JOY 


Proof that D¢o +» [9] = 0: 


Given that 
DES Tt) a0, 
and 
Des +[n] =0, 
to prove that 
DES g] =0, 


where g=m-+n. The understanding is that no lesser multipole tensor of m 
or n vanishes, and that no lesser multipole tensor of g will vanish. 

The proof is carried out in the spherical representation. As g=m-+n, 
it follows that 


[ Vile] = Yt[m] +Yitnl, 


(AIV.1) we D p! VO) x 
| Vlg] ee nt Gn Yo-xLm]|Yx[n], 


because m and n commute. In particular, 


AS a So f È 
Yis[g] dl. +)! Vis im] Yi{n]oc Y »s[m] Yis[n], 
because Yi[n]—0 for K>25, and Woes | nn = 0rfor k= 28% 


Therefore, Yi [g] 0. However, 


eee ai e a | 


k=0 


Shi (4S se ee 
ei Sees Vis +1aCm] Yi[n] 


Di i (45 +1)! 2d tip 2S +1-k 
7A OS F1—ml asp el Fasi] 


UE 


which is equiva lent to 
D§s* "gl = 0: 
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APPENDIX V 


The Tchebichef polynomials. 


The Tchebichef polynomial, p,(m), is given by the diagonal matrix element 
of A™ in a representation in which a-S is diagonal: 


AM 
Pim) = <m|—|m>. 


In 


A™ must be proportional to the spherical tensor Y° which is generated by » 
commutations of at+:S with the spherical tensor Y;,” defined as 


vas sei (an -S)* a 
What -is, 
AME R(at-S, (at-S), ... (at-S), (a=: 5)", 


in which n commutations are performed. This can be written as 


A©= R(at-S, (a~-S)") . 
In general, 
n | 


In] NIE < | —* Xn VX 


so that 


DA x n! +. S)n-k Bo: \ 


The matrix product is reduced to a single term because of the shift nature 
of the operators 


Re 


Ink=0 


n | dre IO a 
Cr ee 


- (mn + k 


Pai(M) 
(a--S)"|m +-k) mt k|(at-S)|m). 


The matrix elements of the shift operators, given in many texts, are 
«m+ 1 |at-S|m)= Cm 5 
im —1|a~-S|m> = Crain 


where 
Jom |? = HLS(S +1) — mm + DI. 
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Therefore, 


n 


R 
palm) =— > (—1) k!( Il |Cmeo| Il | Omo]? « 
In k=0 
The products can be evaluated: 
n SIA (S+m)! 
IT eno — on-k (S-m)! (S+m_-n+k)!° 
k= ee (S—m)! (S+m+k)! 
IT lems! TG (S—m—k)! (S +m)! 


Going on: 
S( 1} n! (S—m+n—k)! (S+m+k)! 
2” Gn > k!(n—k)! (S—m—k)! (Stmt+tk—n)! 


Pr(m) = 


Let m=—S +4, so that 
R 


en n n! Giani Ri ro (gd+ k)! 
Pd = 1) k!(n—k)! (25—q—k)! (q+k—n)! 
Let 
QS = g-n) bg! 
No Gs qi eens 
so that 
( SE n ve n! I 
Prif 7 2" ga (n!) DI k! (n k) q Î na 
TORA i n n! 
Aa 2 ad k PS 
vg AT ira exp |i ta] 0 
sl (ARR, ù te 
= (—1) gigi te (exp 3-3) t(q) = (—1)” I g, (n!)2?Ant(q) 


where 


To simplify the form of t(q), note the identity 


PEt Pe OE a SEN 
Q5—g) © 


which permits 


with p= 25 +1. 
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To evaluate R, realize that the coefficient of (a-S)" in A™ is m!, and that, 
consequently, the coefficient of g” in p,(q) is n!/g,. Now, 


Anfo) = fle + On), 0<0<1, 


by the mean-value theorem quoted by SzEGÒ (19). Therefore, the coefficient 
of q" in (—1)"t(g) is 


2n(2n TIM 1) _ 201! 


(n De (n)? 


It must be that 


n! 
SL 


With the choice of g, as always positive 


{Reino x ‘¢—p 
A) 


APPENDIX VI 


Evaluation of al, Df? n+ 


The recursion relation for a, is 


a; ee = = in| ò; “gl gl 


Multiply by d,b; ... da 


2(n —1) n(n —1) 
(2n —1)(2n — 3) 2 


bb... ba? at Lia beta 
05 Cad n | 


Let ae 
n 
Zn= bb; 0 da (aa = by. nj. * 


ZZZ 
(AIV.1) Z,=6'b2Zn1— on) @n=3) 7°’ 


with Z,j= @a-b. 
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The Legendre polynomials, P,(), are completely defined by the recursion 
relation (1) 


n NA 


@Pya = 5 Pat ‘Paesi 


2n —1 z 


with the condition P,(x) = «. It follows from (AVI.1) that 


n! 
(n —1)!! 


(n) (n) 
Zn = Wij...n04i.,.n = 


(11) E. JAHNKE - F. EmpE: Tables of Higher Functions (Leipzig, 1952), p. 111. 


ate 


SUPPLEMENTO AL VOLUME XII, SERIE X We Up JOSE 
DEL NUOVO CIMENTO 20 Trimestre 


Ultrasonic Investigation of Nuclear Spin-Lattice Relaxation. 


G. S. VERMA 


Department of Physics, University of Allahabad - Allahabad, 
Ultrasonics Research Laboratory, Western Reserve University - Cleveland, Ohio 


(ricevuto il 17 Dicembre 1958) 


Conrents. — 1. Introduction. — 2. Theory. — 3. Experimental technique. 
— 4. Results and discussion. 


1. — Introduction. 


Ultrasonics has been utilised extensively in the past to study the proper- 
ties of gases and liquids which threw light on many molecular processes 
and has proved to be a very convenient tool in the study of relaxation 
phenomena. In solids it has been used very effectively in studying the elastic 
constants from velocity measurements, especially in single crystals. From 
damping and attenuation measurements there have been attempts to seek 
information regarding imperfections- dislocations and radiation damage and 
their kinetics. Recently Ultrasonics has proved to be very helpful in studying 
the properties of conduction electrons in metals, especially when the metal 
becomes superconducting (*). The number of thermally excited electrons can 
be determined from the measurements of ultrasonic absorption in supercon- 
ductors. As a matter of fact, MORSE and Boum (2) claim to have measured 
the energy gap in the electron energy spectrum in superconducting tin and 
indium, which, however, has been questioned by some workers (3). Ultra- 


(1) W. P. MASON and H. E. BomMEL: Journ. Acoust. Soc. Am., 28, 930 (1956). 
(2) R. W. MORSE and A. V. Boum: Phys. Rev., 108, 1094 (1958). 
(3) M. A: BIONDI, À. T. FORRESTER, M. P. GARFUNKEL and C. B. SATTERTHWAITE: 


Rev. Mod. Phys., 30, 1109 (1958). 
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sonic waves have been shown to induce transitions (45) in the spin energy levels 
of a particular nucleus through the nuclear spin-lattice relaxation, the possi- 
bility of which was suggested first of all by ALTSCHULER (°). The possibility 
of a strong interaction between ultrasonic waves and spin waves, especially 
when both the magnon and phonon frequencies and wavelengths are equal, 
has been discussed in ferromagnetic crystals and shown to provide an inde- 
pendent measure of the exchange energy constant (?). The potentialities of 
Ultrasonics in the Gigahertz region, where the acoustic frequencies have wave- 
iengths comparable to the wavelength of light and the phonon energy may 
approach thermal energy at liquid helium temperatures, are many but its appli- 
cation is, however, limited by the inherent experimental difficulties. The 
highest frequency which has been achieved as reported recently by BA- 
RANSKI (8), and BOMMEL and DRANSFELD (°) is only 2.5 GHz (*). 

In many solids nuclear spin lattice relaxation is known to involve the 
presence of impurities. It has been shown by ROLLIN and HATTON (!*!!) and 
by BLOEMBERGEN (??) that for nuclei with a spin J=4, which have no quadru- 
pole moments, the observed spin-lattice relaxation times can be satisfactorily 
accounted for by assuming that they are due to the presence of electronically 
‘paramagnetie impurities. However, for nuclei with a spin J > 4, POUND (18) 
has obtained strong evidence that in pure crystals at not too low tempera- 
tures the quadrupolar relaxation may be more important than the impurity 
relaxation. KRANENDONK (34) has developed a method for caleulating the re- 
laxation time due to the perturbation of this quadrupolar interaction by the 
lattice vibrations. The interaction between atomic nuclei and their sorroundings 
can be studied by nuclear paramagnetic resonance (1°) or nuclear induction (16). 
Such interactions affect relaxation times (1!°!8), line widths (171°) and frequen- 


W. G. Procror and W. TANTTILA: Phys. Rev., 98, 1854 (1955). 
W. G. Proctor and W. TANTTILA: Phys. Rev., 101, 1757 (1956). 
S. A. ALTSCHULER: Dokl. Akad. Nauk USSR, 85, 1235 (1952). 
C. KirtEL: Phys. Rev., 110, 836 (1958). 

K. N. BARANSKIJ: Dokl. Akad. Nauk SSSR, 114, 517 (1957). 


(*) Note added in proof. - Recently JAcoBSEN has extended the frequency to the 
vicinity of 10 GHz (Phys. Rev. Lett., 2, 249 (1959)). 
(°) J. Hatron and B. V. RoLLIN: Proc. Roy. Soc., A 199, 222 (1949). 


(1!) N. BLOEMBERGEN: Physica, 15, 386 (1946). 

(!?) B. V. RoLLin and J. Harron: Phys. Rev., 74, 346 (1948). 

(13) R. V. Pounp: Phys. Rev., 79, 685 (1950). 

(14) J. van KRANENDONK: Physica, 20, 781 (1954). 

(1°) E. M. PurceLy, H.C. Torrey and'R. V. Pound: Phys. Rev., 69, 37 (1946). 
(16) F. BLoca, W. W. Hansen and M. E. PACKARD: Phys. Kev., 69, 127 (1946). 
(1?) N. BLOEMBERGEN, E. M. PurcELL and R. V. PounD: Phys. Rev., 73, 679 (1948) 
(5) N. BLOEMBERGEN: Thesis (Leiden, 1948). i 
(19) E 


. M. PURCELL, N. BLOEMBERGEN and R. V. Pounp: Phys. Rev., 70, 988 (1946). 
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cies (2°21), Excitation of ultrasonic vibrations in a crystal offers another me- 
thod to study the strength and nature of the coupling between lattice vibra- 
tions and the nuclear quadrupole moment. Ultrasonic waves can induce nuclear 
spin transition in solids through the interaction of the quadrupole moment 
of the nucleus, with the electric field gradient at the site of the nucleus, which 
varies periodically with the frequency of sound waves. 


2. — Theory. 


The electrostatic interaction energy F of a nucleus with the charges consti- 
tuting its environment is given by 


(1) 3 [MT 


where the integration is taken over the nuclear volume, o(#) is the nu- 
clear charge density and V(z), the electrostatic potential arising from all 
charges other than those of the nucleus under consideration. Expanding the 
potential V(x) in a power series about the nuclear center of mass, we obtain 


0%; 2 ON; OL, 


a i 277 
<a azoto) {Ts x) DI ca )an+ | - 
j \ 0 j,k \ 


: OV 1 C2 We 
= FeV jo > te) (a, + 5 > di (at 


On; 


where 


fire) — Ze = nuclear charge , 
| d'xe(a)x; =P; = electric dipole moment, 
| da o(#) 0,0, = Vix = electric quadrupole moment tensor . 


The first term represents the electrostatic energy of a point nucleus, which 
is independent of nuclear size, shape or orientation and is of no interest to us. 
The second term which arises due to electric dipole moment vanishes by virtue 
of the fact that the nuclear ground state is non-degenerate and the time 


(29) W. D. KNIGHT: Phys. Rev., 76, 1259 (1949). , 
(21) W. G. PROCTOR and F. 0. Yu: Phys. Rev., TT, 717 (1950); W. C. DICKINSON: 


Phys. Rev., TT, 736 (1950). 
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averaged charge distribution satisfies the symmetry condition o(#)= o(— 2). 
Leaving all other non-vanishing terms in the series except the quadrupole 
terms, the interaction energy is represented by 


(3) a +> Qin Vir , 
Tk 


where V,= (0?V/0a;0x,) and Q,, and V;, are second rank symmetric tensors. 
Vin, however, can be expressed as V,=— 0£,/0%;, where E,== 0V/0%, 38 
the electric field at the nucleus. Thus V, represents the gradient of the 
electric field and interacts with the nuclear electric quadrupole moment. 
However, in order to have correct evaluation of the interaction energy, 
one has to consider the polarization of the charge cloud surrounding the nucleus 
by nuclear quadrupole moment, which has been discussed in detail by STERN- 
HEIMER and his associates (2224). The induced quadrupole moment can be 
much larger than the nuclear quadrupole moment and is such as to reinforce 
the nuclear quadrupole moment. This is called antishielding or polarization 
effect. FoLEYy, STERNHEIMER and TycKko (?5) find that in ions such as CI, 
Cu*, Rb* and Cs* the induced quadrupole moment is 10 to 100 times larger 
than the coupling obtained when there is no polarization effect. KRANENDONK 
has taken this effect into consideration by means of a multiplication factor 


(4) A’ = ast = «| 0(“) V (a) dea, 


where #’ is the interaction energy when polarization of the atomic core is taken 
into consideration. There is another effect which is also present. The direct 
field from the sorrounding charges polarizes the central ion and leads to an 
extra field in a direction opposite to the direct field. It is called the shielding 
effect and its effect on the quadrupolar spin-lattice relaxation is small. 
KRANENDONK also drew attention to the effects of covalency which have 
been discussed in detail by YostpA and Morty (2°). If the lattice is deformed 
by the lattice vibrations, mixing of p- or d-like orbitals into the s-like wave- 
functions of ions takes place and as a result the electrons in these non-spheri- 
cally symmetric orbitals give rise to an electric field gradient at the nucleus 
which interacts with nuclear quadrupole moments. The resulting quadrupole 


M. STERNHEIMER: Phys. Rev., 80, 102 (1950). 
M. STERNHEIMER: Phys. Rev., 84, 244 (1951). 
M. STERNHEIMER and H. M. FoLey: Phys. Rev., 92, 1460 (1953). 

. M. FoLey, R. M. STeERNHEIMER and D. Tycko: Phys. Rev., 93, 734 (1954). 
. Yosipa and T. Morrya: Journ. Phys. Soc. Japan, 11, 33 (1956) 
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spin-lattice coupling can be much larger than with the direct field alone. As 
a matter of fact YostpA and MoRrIvA have shown that the relaxation time 
due to covalency is shorter by a factor of 10-4 than in the pure ionic case. 

A detailed theory of quadrupolar nuclear spin-lattice relaxation has been 
formulated by KRANENDONK for crystals of the NaCl-type. He considered 
one parameter model where the nearest neighbours of a given nucleus are 
replaced with point charges of magnitude q= ye. The magnitude q of these 
charges is the only adjustable parameter of this model. The dimensionless 
quantity y is a measure of the strength of the spin-lattice coupling and e is 
the electronic charge. Values of y are obtained by comparing the experimental 
values of the relaxation time with the values calculated on the basis of the point 
charge model. For example, in ionic crystals of the NaCl-type spin lattice 
coupling is due to interaction of the nuclear quadrupole moment with the 
electric field arising from six equal point charges on the six neighbouring lat- 
tice sites and the rest of the charges are supposed to be absent from the crystal. 
According to his calculations, values of y are of the order 10?— 10°. 
KRANENDONK expressed the interaction Hamiltonian as 


(5) K' = > W,,2,,» 
n 


where Qi are the components of the nuclear quadrupole moment, the equi- 
valent operators of which in the subspace J = constant are given by 


On = 0 bi = 1) 
(6) Q,17 CVS) (L as iI) I.]/2, 
Qui OM +0, 


where 


eQ 


= .Q = CIR ID) - 
mo 3 IONE 


de 


The only non-vanishing matrix elements of @, are ‘m+u|0, |M}. Here mh 
and (m+)h represent the eigenvalues of J, if the direction of the external 
applied magnetic field is along the z-axis. The coefficients of W, are functions 
of the displacement of the nuclei relative to the central nucleus and hence W, 


can be developed in powers of 7;. 


(7) Me. si Aci » Aye TT > Alert ae 


The double dot indicates the complete inner product. 
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Substituting (7) into equation (5), 70' is expressed in powers of the r, 


I 


(8) Ha KH, +H, + Hye, 


where #,= > A,Q,.- 


The onan term 26, gives rise to the static quadrupolar interaction. 
This perturbs the originally equally spaced Zeeman levels of the nucleus in 
the presence of the external magnetic field and leads to splitting of the magnetic 
resonance line. In the absence of external field only pure quadrupole spectrum 
occurs. However, relaxation is caused by the higher order terms Fe) Hee 
The linear term #6, in the first order leads to the so-called direct processes in 
which the nuclear spin makes an upward (downward) transition and one of 
the lattice oscillators is de-excited (excited). In second order the linear term 
Hi gives rise to Raman processes where the nuclear transition is accompanied 
by the simultaneous excitation of one lattice oscillator and de-excitation of 
another one. The first order effect of 3, gives rise to Raman processes. The 
contribution due to first and higher order effects of 76, higher order effects 
of H,, and first order effects of Fe, #H,, ... are all negligible compared to the 
first order effect of the quadratic term in H,, which according to KRANENDOK, 
is the relevant relaxation process. Processes in which more than two Re 
are involved are important only at very high temperatures. 

The transition probability P(m, m+) that the nuclear spin makes a 
transition from the state m to the state (m+) as a result of the Raman 
processes is given by 


P(: 4 Ore È [pa | f(@) ie er 
(9) P(m,m + u) = On? dev? o (el) M,,(ka) do , 
0 


where Q m= <m+p|Q,|m>, d= mass density of the crystal, v= velocity of 
propagation of long wavelength sound waves independent of direction of pola- 
rization, k=1/m the wave number of lattice waves. o(w) is defined such 
that the number of lattice oscillators with frequency between w and w+do 
is given by 


(10) 0o(m) dw = A 2, (0) )k? daw , 


where V represent the volume of the crystal. 
Here 2=how/AT, 


(11) = I DIS pijlm Tie sim(ka) , 


ss’ ijlm 
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with 

2) Lia. tim a) = {B:(k, 8) B:(k, 8)}{ Bk’, 8!) B,(k', 8')}, pijim = Ayes > Agim - 
Here the B,(k, s) are defined by 
(13) B;= B;(k, s) = 6(s, 1)(cosk,a;,—1)+ 6(s, 2) sin k, a, , 


lwhere s is a number which is equal to 1 or 2 and a; = RR, R, and R, re- 
\ presenting the equilibrium positions of the ith and central nuclei, respect- 
ively; k and k’ represent the directions of wave propagation. The curly 
I brackets in the equation (12) represent averaging over all directions of k and k’. 

Das, Roy and GHosH Roy (27) extended Kranendonk’s method to body 
centered cubic structure. 

Yosina and Morrya have calculated nuclear relaxation times in ionic 
erystals on the basis of covalent effects alone. In the covalent state one 
p-electron is transferred from the closed shell of the halogen ion to one of the 
neighbouring metal ions and the hole created by the missing electron causes an 
asymmetric electron distribution which results in a very large field gradient 
at the site of the halogen nucleus. They introduced a parameter 2 which is 
a measure of the degree of covalency. KANDA (?§), as matter of fact, has caleu- 
lated the degree of covalency in each halide from his measurements of che- 
mical shift of halogen ions in several metal halides and shown that the reci- 
procal of the relaxation time is proportional to the degree of covalency. 
Taking as the zeroth approximation the purely ionic state YosmpA and 
i Mortva have calculated relaxation times on the basis of perturbation treat- 
ment. According to them the third order perturbation is given by 


(g | Fé,|n) (| Fq|n') (m' | Hi \9) 
a SENO 


) where g means the ground state and n, n' are the excited covalent states in 
| which one p-electron on the halogen ion is transferred to the neighbouring n 
| or n'-th metal ion. Here &; is the Hamiltonian which causes an electron of the 
halogen to jump to a neighbouring metal ion. Here terms such as (n|€,|') 
are usually very small unless n= n'. Neglecting such terms 


(15) AzE,=> alee (n| H_|m) = ¥ Anln| Helm) = Oud An (ml Wl)» 


n 


(VIVI: Das, Ds Ky Roy and S. K. Guosu Roy: Phys. Rev., 104, 1568 (1956). 
(28) T. Kanpa: Journ. Phys. Soc. Japan, 10, 85 (1955). 
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where 7=|(g|70,|n)}:/(E,—E,)? is a measure of the degree of covalency and 
is regarded as a function of the distance only between the halogen nucleus 
and the n-th metal nucleus. 

The second summation in equation (15), abbreviated as Was is expressed 
as follows 


W, =A’ » A, (n 


(16) We ZA DA (PMI, 
We 1 Ai » ha, (Ln F Vn)? 7) 


where /,, m, and n, represent the direction cosines of R,(X,, Y,, Zn), which 
denotes the co-ordinates of n-th metal nucleus and A'=s eC<1/r®7 

Expanding these expressions in powers of small deviations from equi- 
librium positions and retaining second order terms only 


7 


2A n 
Wi. = — S [wu {nF y,(6R,)? + 2(n-dR,)(yndRa + nTòZ,)} + 


db 
+ onTy,(n-6R,) + A6Z,6R7], 
(17) re 
Wass = gia Su {(n*)(OR,)? + An ORE (n° dR.)} + 


+ o(n*)?(n-dR,)? + 2(n7)? (ORF), 


where u=—4ad'— A, v=44—5aX'++a?7". 

n, Pn and y, are the direction cosines of the equilibrium position of the 
n-th metal ion measured from the halogen nucleus, 6R,(dX,, 0Y,, 6Z,) the 
displacement vector, a=|R,0|,} n is the unit vector whose components 
are on, Bn and yn, OR=>=6X, +1 6Y, and a,+if,=—n*. 

Following Kranendonk’s notation, the transition probability is given by 


(e 


A” On N Ta ex 
(18) P(m, m + ps) = Se i] Di Tenia Mala) de 
0 
where 
(19) Da (ka) = > DI {Bno Bag’ Be Bio} > jin (00) O py (o0') b) 
n,n’ s pp' 


where o stands for a set of k (wave factor), p (one of the three directions of 
polarization) and s, s taking values 1 and 2, 


B,(0) = 6(s, 1)(cosak-n—1) + d(s, 2) sinak-n 
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and 
2 01,(00')t0x(00*) 4 [15 (ttn + Ba Bae) Yar + Ant + BuBu)(nn) + 
+ 4(af + B3)y2 + dynyn(n-m')} + {on On + BnBn)Ynyw (me n)} + 
AAP uv {nn + PnBn)PnYn + (Ann + PaBa)ya (n) + 


+ 4(0È + B2)ynyn(me-me’)} + du {2(at + BR) + Ay} + 220 + BV 


(20) 


> (5,(T0") Man (00°) a ua ees | (0%, Bn) Xn ot Bo) +4 2( (Cn Xn» st a Pa) 


pp’ 


sin A (On Xn pù Da ae sine Bi) (oe, bi) i vt (xp, Ba + Bi) + 
SF 2 (Gn Xn: + Bn Bn) + A? 4 1 V {4 (Gn tpn)? (Oni + tBa)? + 
A. 2 (cn om» + Bn Bu)(a2 + B2)(me-10')} + AAU (x + Br) + 4 4005 + Bn)? 


which can be simplified further as (n-n’)=0 for n# + ‘andere ad 
for n= + n' crystals of the NaCl type. The transition probability, however, 
is given in the following convenient form 


4 
(21) Poms oi Ga LEAD): 
n=1 
where 
On = = 31 Q im P/aed? 0? a? (KRANENDONK) , 


= A? 08 |Q im Fl [2a adv (YOSIDA and MORIYA) 


Shi 0 Debye temperature , 


IT 
ue 


D,(T*) = T* i eas L,(eT*«) da , 


C = lim = (6n?)*, Ho = hvkm = KO , c= (Ww3 1/2)? 


for NaCl-type for CsCl-type 


D,(T*) has been expanded in descending powers of T* by KRANENDONK and 
is given as follows (for NaCl or simple cubic crystals) 


NE 0.20 — 0.0071/T*° +... 
DP = 0920. 055/T* + 

VADA Rat VASTA bp ee. 
D,(2*) = 0.28 — 0. 017/T* +... 


4 - Supplemento al Nuovo Cimento. 
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The coefficients NV, are given by 


Lo 
ND (EN ian Ie) Uap (n= 2, 3, 4), 
Nig ARIA 
A= 0 fori |A Ae Ae SAT, 
fee | CS Ed, CO 


Ass A (0 be 0159 AA 00 a A de 


where a is the equilibrium distance between the nearest neighbours. According 
to YosmpA and MORIYA, the coefficients N,_ are given by 


un 


* 
Nia A NOI sa SI > OTS )O jy (GO) o) 
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which after somewhat lengthy calculations are expressed as functions of 7, 
A’ and 4" 
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The degree of covalency 4 is a function of the distance between the ions and 
is estimated from the chemical shift of the halogen nucleus. Assuming 
Aa exp[— R/0.345], where R is the interionic distance in A units, such that 
M/A =—1/0.345 and 7"/A= (1/0.345)?, Yostpa and MorIvA have calculated 
the coefficient N,,. 

The relaxation time t can be defined by t= 1/P, where P is the average 
value of P, averaged over all values of m and mw as the relaxation time in the 
presence of quadrupolar interaction is not unique. However, Kranendonk’s 
values of relaxation times in the absence of covalent and polarization effects 
(y=1) are about 4 to 6 orders of magnitude larger than those observed. In 
order to account for the experimental relaxation times, he assumes mul- 
tiplication factor y of the order of 10? 10%. The values of nuclear relaxation 
| times in ionic crystals calculated by YosIipA and Morrya on the basis of 
covalent bonding effects alone are about 3 to 10 times larger than the ob- 
} served ones. This discrepancy is traced in the simple Debye treatment of the 
| lattice vibrations and ambiguities in the determinations of the Debye tem- 
perature 9 or sound velocity v. The interaction of the covalent electron with 
the quadrupole moment is assumed to be unaffected by shielding and anti- 
| shielding effects and the relaxation time is attributed to covalent effects alone. 
Kraus and TANTTILA (2°) have computed the probability per unit time 
that the sound waves will produce a transition of a nuclear spin between its 
Zeeman levels. The nuclear quadrupolar interaction Hamiltonian is given by 


2 
>) ca > Al > C:.Q: Be, , 
dea 


where g, is an external charge to a nucleus N, 7 is the distance between 
the nuclear center of mass and the external charge q, when both of them are 
displaced by sound waves, 0, are numerical constants, Q; are the components 
of the nuclear quadrupole moment given by @; =O, (r)A,dt,, (= 0, + I E2), 
The quadrupolar term of the electrostatic potential at a position 7, from the 
center of mass of the nucleus, where o, is the nuclear charge density due to the 
charge q, after displacement, is given by 


(26) Vise AR, 


where A and B are second rank symmetric tensors. 


(29) O. Kraus and W. H. TANTTILA: Phys. Rev., 109, 1052 (1958). 
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The matrix elements of the operator which are of interest are given by 


(Lie FC | Les) = (Lal 142) LABS, 


(In| |Im—a) = (Em|Q+2| Ina) > a 0,B®,. 


The displacement of a lattice point in the presence of unidirectional standing 
sound waves being propagated along the x-axis is given by 


S = Acos(ke—d)coswi, 


where k= 2z// with A the wavelength of the sound waves, A the amplitude of — 
the sound waves and 6 an arbitrary phase factor. The displacement of the 
nucleus is taken as Sî— A cos écoswt and for the charge gq, as SP = 
= A cos (ka — 6) cos ot. 

Expanding B*,/R, and B*,/R° in terms of the relative displacement of q, 
and utilizing the fact that 


3e0) 
(Im|Q+ ey) a Teer ee m 2) (1 Mm LA m 1)(I m)|*, 
(28) 360 
(Im |Q-2|Inta) = rr pl U+#+2)A-m-1)T+m+1)A-m)}, 


the matrix elements are given by 


arie 
(In| 2€|Imt) = sro aj 
T(E + m + 2)(L— m —1)(L +m +1) (Lm) == Ak sin 8 cos at 
(29) 
eye 
(In| Fe|Ima) = ror — 7 


‘[(L—m + 2)(L 4 m—A)(Z~-m +1)( +m) Ab sin 6 cos wt. 


where q,=ye. 

From the first order perturbation theory the probability that in time ¢ 
the nucleus will make a transition from the initial state m to one of the states k 
is given by 


SI | ato |? 1 


(30) Pai È OS 


i, 


he 
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where 76° is time independent and g(w) is the normalized shape function of 
the nuclear resonance line. Averaging over the frequency spread 6 of the 
resonance line the average transition probability per unit time P is given by 


(31) P =| Hin |?/40? dv . 
Thus the transition probability per unit time for Am= + 2 transitions is 
given by 


SOE Ey aes 
~ 64720) 12(21 —1)? 
Ake, 
-(I—m + 2)(I + m—1)(I—m + 1) (I pm) sin? d, 


P(m, m — 2) 


(32) 
81 ey? 
— 64%? dy I°(2I — 1) 


:(1+m+2)1-m—-1)(I1+m, 1)(L m) = sine d. 


P(m,m + 2) 


Kraus and TANTTILA have also derived expressions for the equilibrium 
nuclear magnetisation in the presence of ultrasonic and thermal lattice vibra- 
tions. The time rate-of-change of the e-component, m,, of the macroscopic 
magnetisation, which arises from the nuclear spin system is given by 


dm, dm, | dm, | dm, 
(33) na —— SEI a SS 
dt dt |im=+1 dt im= +2 dt |ultrasonic 
where 
dm, | 6(m., — Mz) ARSA m 
a de Om_19 


di ngi II + 18) (21 a5 1) m=1 
Il 
my = VM + 1n(ho]KT), 


— value of m, when the spin system is in 
thermal equilibrium with the lattice, 
T,= Temperature of the lattice. 


Here ©? , is the probability per unit time for a thermally induced transition 
of a spin from the state m to the state m — 1 and given by 


90° (2m —1)(E +m)(T-m+1)F0), 


fe mori 


where 71) is a funetion of the lattice co-ordinates only. 
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Similarly (dm,/dt)| m1. i8 given by 


35 dm, | 24 (Mz, — ™m,) Sx. 
( 5) di AO Ge + itl CE ci iB) cai m_-2 9 
where 


9 2()2 
OR: pr ap E+ MA + IY + 11m +240, 


and again G@(l) is a function of the co-ordinates only. 


If the frequency of the sound waves is chosen in such a way that it pro- 
duces Am =-— 2 transitions 


(36) dun —24 © pa 
dt near? p To + 1) (27 4 1) i m=I a 


where the ultrasonic probabilities are given by equations (32) 
Introducing relaxation times dm,/dt can be expressed as 


Es dm, m,—M, Mm, 
197) Gigi T,) 
where 
LI 1 
TR gi e IRE 
1 6 = 
TI) (E 1) mom) 
(38) 
1 24 -142 
POI) F1i) =," 
1 24 2 oe 
T, I +) Cl +1) Di - 


Here 7, is the observed spin-lattice relaxation time. 
time can be computed from equations (38) and (32). 

If the ultrasonic excitation which is applied to the specimen is constant 
and the time through which it is applied is longer than 7,, the magnetization 
is constant in a region where d = constant. Thus putting dm,/dt= 0, we obtain 


Ultrasonic relaxation 


Mas 
(39) My = = È = 


ENI Ren 
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or 


where 


(40) ee 


The total magnetisation is obtained by integrating equation (39) over the 
volume of the crystal and is given by 
Ada M, 


41 M,= ho — Mm. = = pel tI 
È i | mS ye gain] (pl 


where A is the cross-sectional area of the uniform cylinder, dò = 2rx/A, 
7= ultrasonie wavelength, M, = total magnetisation of nuclei in the absence 
of ultrasonic excitation. 

ABRAGAM and Proctor (3°) have also deduced the ratio (M,/M,) without 
going into the exact expressions for the contributions to the rates of change 
of populations from the spin-lattice relaxation and spin-spin coupling. The 
rates of change of population of the levels 7, =m of 2*Na irradiated by ultra- 
sonic waves at twice the Larmor frequency are given by 


"i = — W(ny — N_4) + Fs (Mm) + 3 (Mn) ; 
d a 
=: = — W(m— 0-3) + Fa(Mm) + Gal%m) è 
(42) 

du_i 

di 2 = — W(n_y— Mg) + F_4(%m) + Gal%m) 5 
dn_s 

Ti = = — W(n_2— 3) + F_4(%m) + G_4(%m) 5 


where 7, and G,, represent the contributions to the rates of change of the 
populations from the spin-lattice relaxation and spin-spin coupling, respectively. 
Writing (M>= (ny n_3) + H(%|— 24) the contribution of spin-lattice re- 
laxation to the rate of change of <M,>, 1.€., 


provided a single spin-lattice relaxation is assumed, we obtain 


d(M,) sr» ALE al —2W{(n- n) + (ny— N3)}; 


(43) dt via 


(39) A. ABRAGAM and V. G. Proctor: Phys. Rev., 109, 1441 (1958). 
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as the contributions due to spin-spin interaction is zero because of conser- 
vation of energy. As Boltzmann distribution is maintained throughout 


(Mg — n_3) = 3(m — n_;) and consequently <M.) = Lo (NM, —N_3) 
or 
d (M,»—M,, 8W i 
CAMS EEA ZZZ 
(44) dt ik Z/ ‘TA 5 


The limiting value of <M,) is given by putting <dM./dt) = 0 


(M,>|M,, = (1+ 8WT,[5). 


Since the signal amplitude £ is proportional to <M,>, we have 


8WT,\-! 
(45) Elfe= Mo], = (14S) 


3. — Experimental technique. 


The technique is essentially based upon the pulsed nuclear induction method 
of measuring nuclear magnetization as described by HAHN (81). The thing to 
be remembered is that if a short pulse of radio frequency magnetic field is 
applied at the Larmor frequency, the transient nuclear induction signal is 
proportional in magnitude to the static magnetization of the system, which 
is a function of the population of the various states. Thus the amplitude of 
the transient nuclear induction signal A, following the pulse of radio frequency 
flux is proportional to the population difference between the quadrupole states 
at thermal equilibrium. Once the equilibrium is disturbed the instantaneous 
population difference A will try to reach its equilibrium value A, at a rate 
which is determined by the thermal relaxation time 7, according to the 
relation 

A(t) 


(46) ay (1— exp [—4/T,)), 


provided the time interval between the two successive radio frequency pulses 
is greater than the thermal relaxation time 7, and A(0)=0 which is 
the case when the pulse length 7 is adjusted to give maximum transient in- 


(71) E. L. HAHN: Phys. Rev., 77, 297 (1949). 
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duction signal. If A, is the amplitude of the transient nuclear induction signal 
after a second pulse of radio frequency flux of the same duration is applied 
after a time ¢ (see Fig. 1), then 


A, 


(47) Tee (1— exp [—#/Z,]) . 


By measuring the ratio (A,/A,) as a Ae 
function of the time interval ¢ and CLI 
plotting log, (1— A,/A,) vs. t, the ther- Suo 
mal relaxation 7, is known from the Fig. 1. 
slope of the straight line. 

Next the ultrasonic excitation is applied at twice the Larmor frequency 
and the magnetization is again measured. The typical time sequence of events 
is shown in Fig. 2. If & is the amplitude of the nuclear induction signal in 


the absence of ultrasonic excitation and &, in the presence of ultrasonic 
excitation, the ratio can be ex- 


RF pulse pressed as 


Teas) (2) - 14+ 47%, 


Fig. 2. where V is the peak voltage ap- 
plied to the transducer. If we put 
M,=&, M,=¢ and dM./dt= 0 in equation (41), we obtain 
243 GQ yt Ark? 
4 DM eae SR ANNE ARIE OE 
po) eo; TP) (Ot tol) (Oh = le 
(I—m+2)\f+m—1U—m4 1)(1 +m), 


where dy is the ultrasonic line width, A is the amplitude of ultrasonic vibra- 
tions, k= 2a a/A, 4 is the ultrasonic wavelength, a is the lattice constant, e is 
the electronic charge and the nearest neighbours of a given nucleus are re- 
placed by point charges of magnitude ye. 

Thus from the experimentally measured values of X, (which is known from 
the slope of the curve (&/&) vs. v2), V, òv (which is known from the plot of 
(E/E) vs. ultrasonic frequency at constant power around 27, % being the 
nuclear magnetic resonance frequency), k and 7,, one can obtain the infor- 
mation about the parameter y provided the amplitude of ultrasonic vibrations 
in the crystal is known. Ultrasonic energy density in the crystal is calculated 
from the power transmitted to the crystal and the phonon relaxation times 
of the ultrasonic phonons from the relation 


(50) Vie = PT, 
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where V’ is the volume of the specimen, € is the ultrasonic energy density, P, 
is the power delivered to the specimen, 7’, the sound attenuation coefficient 
or phonon relaxation time. The power delivered to the sample is known in 
terms of the voltage across the quartz and is given by 


Ret 
5 DEALERS 
a VE A 


where R* is the effective load resistance and Rf is the effective internal 
resistance of the quartz and its circuits. The sound attenuation coefficients 
are known from measurements of HUNTINGTON (22) and GALT (3%). PROCTOR 
and RoBINSON (#4) have also 

measured it by studying the 

Plate attenuation of the nuclear in 
aR as duction signal as a function of 

the frequency of ultrasonic ex- 
citation both when one face of 
Ee pe the specimen is or is not in con- 
out of page tact with castor oil. In the for- 
mer case since some of the 
acoustic energy escapes into oil, 
the attenuation of the nuclear 
induction signal is less. From a 


Ultrasonic 
generator | 


Quartz 
NaCl sample 


AN/APS 


receiver 


knowledge of acoustic impe- 
dances of NaCl and castor oil 
and from the loss of attenua- 
tion of the nuclear induction 
signal they calculated the value 
CEI e: of 10 MHz phonon relaxation 
time which is in good agree- 
ment with the values obtained 

by HUNTINGTON and GALT. 
The block diagram of the apparatus which has been used by Proctor and 
his associates is given in Fig. 3. It consists of two main groups, one which 
is used in producing ultrasonic excitation and the other which measures the 
nuclear magnetization. A nuclear induction head of the type described by 
WEAVER (*5) is used. The transmitter and the receiving coils are at right 
angles to each other and a static magnetic field is applied perpendicular to 


| Gated R.F 
power 
oscillator 


Ina Bs 


(*2) H. B. HunTINGTON: Phys. Rev., 72, 311 (1947). 

(°°) J. K. Gatr: Phys. Rev., 73, 1460 (1948). 

(**) W. G. Proctor and W. A. Rospinson: Phys. Rev., 104, 1344 (1956). 
(°°) H. E. WEAVER: Phys. Rev., 89, 923 (1953). 
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the axes of both the coils. The quartz crystal is glued to one end of an NaCl rod 
of about 4 in. diameter and one in. length which serves the purpose of an ex- 
tension the of ultrasonie path length for convenience. The other end of this rod 
is cemented to the sample crystal of the same diameter. The sample crystal 
is well insulated from the ultrasonic circuits. PRooToR and ROBINSON used 
a steady magnetic field of 4220 oersteds and a short pulse of 50 us of radio 
frequency magnetic field at the Larmor frequency 4.75 MHz for studying the 
transient radio frequency nuclear signal whose magnitude is proportional to 


‘the nuclear magnetization just before the pulse. Ultrasonic excitation was 


applied nine seconds after the preceding nuclear magnetization measurement, 
at twice the Larmor frequency for 8 seconds and after a delay of 0.03 second 
the nuclear magnetization was again measured (see Fig. 2). 


4. — Results and discussions. 


The data on ultrasonic investigation of nuclear spin-lattice relaxation time 
are very limited. The only crystals in which ultrasenic excitation has been 
‘used to study its effects on nuclear spin-lattice relaxation are NaCl, Nal, 
Na010, and InSb. NaClO, was the first crystal to be studied by PROCTOR 
and TANTTILA. They observed that the population difference of the two pure 
quadrupole levels decreased when the sample was subjected to ultrasonic 
energy at the transitions corresponding to Am=-+ 1 and Am==2 2. 
However, it was realised that the magnetic fields due to currents in the 
ultrasonic system might affect the population difference. Consequently 
Procror and RoBINSON used ultrasonic excitation corresponding to Am=+2 
transitions in the case of NaCl and the specimen too was placed in a static 


magnetic field. 

Assuming that the quadrupolar relaxation is the most important mecha- 
nism in NaCl, Procror and ROBINSON calculated y on the basis of Kranen- 
donk’s six point model. The value of y comes out to be equal to 3 which is 
somewhat lower than the value of 50 estimated by KRANENDONK. However, 
as he points out, his value is a measure of the second derivative of the electric 
field gradient whereas the value of y which is calculated from ultrasonic measu- 
rements is a measure of the first derivative. The discrepancy can also be attri- 
puted to the lack of reliable measurements of ultrasonic attenuation as re- 
ported in the literature and in the ambiguity in equating phonon relaxation 
times with ultrasonic attenuation. In view of the fact that the knowledge 
of the amplitude of ultrasonic vibrations in the crystal helps a great deal in 
determining the strength of the coupling, it should be known more accurately 
than what is possible now. However, if one does not care for the absolute 
value of the quadrupolar coupling the problem of determining the correct 


ci 
Dci 
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value of ultrasonic energy or amplitude in the crystal can be avoided by com- 
paring the net magnetizations of two sets of nuclei in equilibrium with thermal 
and lattice vibrations caused by ultrasonic waves in the direction of a constant 
external magnetic field H,. This was done in the case of Nal crystal by 
JENNINGS, TANTTILA and Kraus (89). They determined experimentally the 
equilibrium magnetizations in the direction of the static external magnetic field 
of sodium and iodine nuclei as a function of the ultrasonic power delivered 
to the sodium iodide crystal which in turn is proportional to the square of 
the peak voltage applied to the quartz crystal. By plotting (5/8)? vs. V for 
both sodium and iodine nuclei, the slopes ky, and k, are known. However, 
the ratio of the two slopes is given by 


(52) la Tx Qiyt Ome 9.9.40) , 
Kya Twa dvi QNaVNa 


provided we assume that the amplitude of ultrasonic vibrations A 1s propor- 
tional to the applied voltage (A = CV, where C is a constant). JENNINGS, 
TANTTILA and KRAUS obtain the following values experimentally: 7,;= 0.0065 s, 
Tig = 58, (Kr/kx,)= 2.2, Ory, = 4.65 KHz/8, 6v,; = 4.49 kHz/s. Substituting these 
values together with values for quadrupole moments for sodium and iodine 
nuclei, they obtain 


(0'r/Yxwa) = 10.9, 


which is in qualitative agreement with the theoretical calculations of STERN- 
HEIMER (37) and DAs and BERSOHN (35). However, k is to be compared with 
1+y, Where y, is the polarization calculated by STERNHEIMER. According 
to Sternheimer’s calculations, y, for sodium ions is 4.2. Though he does not 
give any calculation for iodide ion, yet it can be guessed from the corresponding 
value in the cesium ion, in which case he obtains y_ = 143. Supposing that 
the y, for the iodide ions is not less than 143, (y,/y,,)=28 according to 
Sternheimer’s calculations. The discrepancy between the two results can, 
however, be traced in the simplified theoretical model which has been used 
by JENNINGS, TANTTILA and KRAUS. Moreover, imperfections in the crystal 
which are known very well to cause broadening effects on the satellite lines 
of iodine nuclei should also be taken into consideration. Imperfections in the 
crystal give rise to random but permanent electric field gradients which inter- 
act with nuclear quadrupole moments. This might contribute in the right 
direction towards the increase in the ratio of (y,/y,,). However, in the treat- 


(°°) D. A. JENNINGS, W. H. TantTILA and O. Kraus: Phys. Rev., 109, 1059 (1958). 
(*?) R. M. STERNHEIMER: Phys. Rev., 105, 158 (1957). 
(88) T. P. Das and R. BERSOHN: Phys. Rev., 102, 733 (1956). 
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ment given by KRAUS and TANTTILA, the spin system has been assumed to have 
a uniform temperature, which may not be valid in view of the fact that ultra- 
sonic excitation causes hot spots in the system and the spin conduction cannot 
be neglected. 

The above calculations are based on the concept of a spin temperature 
as distinct from lattice temperature whose validity has been discussed in great 
detail by ABRAGAM and PROCTOR (3°). The existence of spin temperature as 
a valid assumption leads to certain predictions which can be tested experi- 

mentally. For nuclei 7/>4 the equidistance of the levels /.= m of a spin in 
ca magnetic field H, makes it possible to define the spin temperature according 
to the relation 


P, > exp {— yhH,m/kT], 


where the populations P,, of these levels form a Boltzmann distribution. 
The Boltzmann distribution of the levels is maintained and even if it is not 
so initially at t= 0, spin-spin interactions bring it to that form. The test 
of such an effect is provided by ultrasonic saturation experiments. Equa- 
tion (45), which is derived on the assumption of Boltzmann distribution shows 
that M tends to zero as the ultrasonic transition probability becomes suffi- 
ciently strong. ABRAGAM and Proctor have been able to increase the den- 
sity of ultrasonic waves in the crystal to such an extent that the magnetic 
moment of 2*Na completely disappears. Ultrasonic experiments illustrate the 
role of the spin-spin interaction in establishing and maintaining a Boltzmann 
distribution between the populations of the spin system. 

Procror and RoBinson have also attempted to seek information regard- 
ing the number of imperfections in the crystal from ultrasonic line width. 
ROBINSON (2%) has shown that the average field eradient caused by dislocations 
is approximately given by q~ bey6*/a?, where ¢ is the number of dislocations 
per square centimeter. Equating the line width to twice the pure quadrupole 
interaction energy eqQ/2, the line width is given by 


be yQe? 


Ay~ th 


Procror and ROBINSON obtain the ultrasonic line width in NaCl equal to 

4kHz. Using this value ¢ comes out to be equal to ~ 10” dislocations/cm?, 

which is in fair agreement with the values obtained by other methods. 
MENES and Bower (4°) studied the absorption of acoustic energy at the 


(39) W. A. RoBINSON: Thesis (University of Washington, 1956), unpublished. 
(49) M. MenEs and D. I. BOLEF: Phys. Rev., 109, 218 (1958). 
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resonant frequency of 1°In nuclei in a single crystal of InSb at both the 
normal nuclear magnetic resonance frequency (Am= +1) and at the double 
frequency (Am=-+ 2). The nuclear resonance acoustic absorption was studied. 
by the change in the mechanical Q of the sample, which changed the electrical 
impedance of the transducer which controlled the Q of the tank circuit of a 
modified Pound-Watkins type of spectrometer. 

So far the ultrasonic excitation has been applied in only one fixed direction, 
say the 100 direction, but if it is applied in different directions (*), the picture 
which we may get with regard to anisotropy of y in the crystal may yield 
interesting information on electron distributions and their relative displace- 
ment under ultrasonic excitation. It will be further interesting if ultrasonic 
excitation is applied to solids where spin-lattice relaxation is due to magnetic 
dipole interaction. So far as the improvements in the technique are concerned, 
some better ways are yet to be devised which may give reliable values of 
amplitudes of ultrasonic waves in the crystal or displacements of lattice points. 
under ultrasonic excitation. 


Part of this work was carried out while the author was at Western Reserve: 
University, Cleveland, Ohio, U.S.A. The author gratefully acknowledges a 
grant from Office of Naval Research which made his stay in Cleveland pos- 
sible. The author is also grateful to Professor ERNEST YEAGER for his helpful 
comments and suggestions. 


(") Note added in proof. — Recently TAyLoR and BLOEMBERGEN (Phys. Rev., 113. 
431 (1959), have studied nuclear spin saturation by Ultrasonics in NaCl crystal. They 
found no angular dependence in the saturation curves. They have been successful 
in measuring the components of the fourth order tensor 8S, which connects the electric 
field gradient tensor at the nucleus with the strain deformation tensor e 


ViVi =F > Siyxt eer > (i, 9, k, l= a, Y, e) - 
kt 


Absence of angular dependence implies that Sj,4/S,;,=%, which differs from Cauchy 
relation S44/Sj;= —3. Also the experimental ratio (8,,)ce/($11)~a= 1.8, whereas the ionic 
point-charge model coupled with an isotropic Sternheimer antishielding factor gives a. 
factor 10. The results show inadequacy of ionic point-charge model and indicate the 
possibility of a considerable amount of covalent character and configurational inter- 
action on the Nat ion. 
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Die Zeit in der Quantenmechanik. 


F. ENGELMANN und E. Fick 


Institut fiir Theoretische Physik der Techmischen Hochschule - Munchen 


(ricevuto il 22 Dicembre 1958) 


InuHaLr. — 1. Finleitung. — 2. Bildunabhingige Bewegungsgleichungen 
im Hilbert-Raum. — 3. Spezielle Bezugssysteme. — 4. Der Zeitoperator. — 
5. Diskussion. 


1. — Einleitung. 


Die Beschreibung der zeitlichen Veranderung des Zustands eines quanten- 
mechanischen Systems erfolgt in der Schrédinger-Darstellung durch die Be- 
wegung des Zustandsvektors in einem Hilbert-Raum §, wàhrend in der Hei- 
senberg-Darstellung der Zustandsvektor zeitlich konstant gehalten wird und 
die Operatoren in § gerade eine solche Bewegung erfahren, daB in beiden 
Darstellungen die relative Bewegung des Zustandsvektors gegeniber den Ope- 
ratoren dieselbe ist. Durch eine unitàre Transformation ist es bekanntlich 
moglich, von einem Bild in das andere zu gelangen. 

In der vorliegenden Arbeit wollen wir darauf hinweisen, daB eine bildun- 
abhingige Formulierung der Bewegungsgleichung des Zustandsvektors mog- 
lich ist, die man unmittelbar aus der statistischen Deutung der Quanten- 
mechanik erhalt, wenn man die allgemeine Quantisierungsvorschrift, d.h. die 
Ubersetzung von Observablen in Operatoren hinzunimmt. Die zeitlich auf- 
einanderfolgenden Zustande werden in der Quantenmechanik bekanntlich durch 
einen Zeitparameter # beschrieben, der eine gewòhnliche Zahl (c-Zahl) und 
kein Operator ist. Die pildunabhingige Formulierung der Bewegungsglei- 
chungen fiihrt uns darauf, da® neben dieser e-Zahl-Zeit t noch ein Zeitoperator T 
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fiir die Observable « Zeit » (1) einzufiihren ist, dessen Erwartungswert ¢ ist. 
© ist explizit unabhangig von ¢, falls der Hamilton-Operator explizit von t 
nicht abhingt. Durch diesen Operator wird es méglich, die Unbestimmtheits- 
relation zwischen Energie und Zeit direkt zu beweisen, was in der ùblichen 
Formulierung der Quantenmechanik, in der die Zeit nur als c-Zahl vorkommt, 
auf Schwierigkeiten fuhrt. 


2. — Bildunabhangige Bewegungsgleichungen im Hilbert-Raum. 


Um die Bedeutung der Zeit in der Quantenmechanik verstehen zu kònnen, 
ist eine genaue Analyse der Quantisierung erforderlich. Zur quantenmecha- 
nischen Behandlung der physikalischen Erscheinungen gelangt man, wenn man 
von den folgenden Forderungen ausgeht: 


1.) Es besteht eine ein-eindeutige Zuordnung zwischen den Observablen 
F und den hermitischen Operatoren ¥ eines Hilbert-Raums . 


2.) Klassische Beziehungen zwischen den Observablen gelten fiir die zu- 
geordneten Operatoren, wenn man im Rahmen der Hamilton-Jacobischen 
Theorie die Poisson-Klammern {#G} durch die Kommutatoren (i/h)[FG]_ ersetzt 
und Produkte symmetrisiert. 


3.) Es soll eine statistische Interpretation méglich sein; d.h. es existiert 
ein statistischer Operator, der den Zusammenhang zwischen den Operatoren 
und den Erwartungswerten der Observablen vermittelt (2). 


Die Ubersetzung der Hamilton-Jacobischen Theorie liefert gemaB Forde- 
rung 1.) und 2.) fur jede Observable / eine Operatorfunktion 


(1) FFG Pes ths 


wobei g, und p; die Orts- und Impulsoperatoren darstellen, wàhrend die Zeit #, 
sofern sie explizit auftritt, lediglich als eine gewéhnliche Zahl (c-Zahl-Para- 
meter) eingeht. Es erweckt dies den Anschein, als ob der Observablen « Zeit » 
im Gegensatz zur Forderung 1.) kein Operator zugeordnet wird. Wir werden 
jedoch im Abschnitt 4 zeigen, daB ein Zeitoperator existiert, obwohl in (1) 
die Zeit nur als c-Zahl auftritt. 


(') Wir beschrànken uns in der vorliegenden Arbeit auf die nichtrelativistische 
Quantentheorie. In einer relativistischen Formulierung iibertragen sich die Uber- 
legungen auf die Eigenzeit. 

(2) J. von Neumann: Mathematische Grundlagen der Quantenmechanik (Berlin, 1932); 
G. Lupwic: Die Grundlagen der Quantenmechanik (Berlin-Gottingen- Heidelberg, 1954). 
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Auf Grund der Forderung 2.) gilt, wenn H die Hamiltonfunktion ist, die 
\Korrespondenz 


iDabei bezeichnen wir entsprechend der Zuordnung 


ni 


ar OF 
Ot 


‘mit F=D den Operator, der nach 1.) der Observablen F= D entspricht 
i(z.B. mg =p). F beschreibt also keine zeitliche Verànderung (« Bewegung ») 
(von F in Him Gegensatz zu (¢F/Ct),,, das die Differentiation von ¥ nach 
seiner expliziten Zeitabhingigkeit bei festgehaltenen Operatoren p und q be- 
deutet (*). Eine Bewegung der Operatoren in © kommt also lediglich durch 
ideren explizite Zeitabhangigkeit zustande. Die Bewegung von F ist nach 
|Gl. (2) durch diejenige von F bestimmt. 

: Die Giiltigkeit der Gl. (2) ist an die Existenz eines Hamilton-Operators 9 
igekniipft (« reversible Bewegungen »). Sie gilt deshalb bekanntlich nicht fir 
‘ein Teilsystem, das Wechselwirkungen mit weiteren quantenmechanischen Sy- 
istemen hat, und auch nicht fiir Messungen. 

| Ein reiner Zustand eines quantenmechanischen Systems wird durch einen 
|Zustandsvektor ¥ in 9 beschrieben, dem ein Projektionsoperator My zu- 
igeordnet ist, der vermoge Dyf= (7, f)\¥ auf ¥ projiziert und der den sta- 
\tistischen Operator des reinen Zustandes nach Forderung 3.) darstellt. Dem 
| Operator Dy entspricht als Observable die Frage, ob sich die Zuordnung vom 
\statistischen Operator Py (bzw. vom Zustandsvektor Y) zum Zustand des 
Systems zeitlich andert. Dies ist definitionsgema& nicht der Fall, so dab 


1(4) Py =O 
| gilt, weil unmégliche Eigenschaften dem Null-Operator zugeordnet sind. Aus 


l denselben Griinden gilt auch fiir den allgemeinen statistischen Operator 0 


l einer gemischten Gesamtheit 


(5) =e |, 


0 
von Ql. (2) lassen sich fir F = OF /Ot Rechenregeln 


(3) Vermége der rechten Seite 
rentiationsgetzen entsprechen. 


aufstellen, die den normalen Diffe 
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fn die Zuordnung des statistischen Operators zu dem Zustand der durch 

ihn beschriebenen Gesamtheit bleibt auch allgemein dauernd bestehen. Fur 
reversible Bewegungen erhilt man nun als Spezialfall der Gl. (2) die Be- 
wegungsgleichung des statistischen Operators in als 


i OW 

a Bg 
(6a) 7 [ACW |_ + | x È 
fiir ein Gemisch (4) und speziell 

4 ODy 

-T[KR, eee 
(6b) (AD + ( È E 


fiir einen reinen Zustand. 
Wendet man die Gl. (65) auf irgendeinen Vektor aus © an, so findet man 
die Bewegungsgleichung des Zustandsvektors 


(7) 


Es ergibt sich damit folgender Sachverhalt: Alle bisherigen Uberlegungen 
haben sich ohne ein bestimmtes « Bild » formulieren lassen, d.h. die angeschrie- 
benen Gleichungen sind unabhaingig davon, welches Bezugssystem im Hilbert- 
Raum man zugrunde legt; sie gelten also alle sowohl etwa im Heisenberg- als 
auch im Schrédinger-Bild (Abschnitt 3). Als einzige Zeitabhangigkeit tritt bei 
allen Operatoren nur die explizite, d.h. jene bei festem p und g, auf. Fiir die 
statistischen Operatoren ist diese explizite Zeitabhangigkeit wesentlich. Gerade 
durch sie wird die zeitliche Anderung des Zustandes einer Gesamtheit erfaBt. 

In Anwendung der obigen Gleichungen laBt sich etwa das Ehrenfestsche 
Theorem bildunabhangig als direkte Folge der Forderungen 1.) bis 3.) erkennen. 
Der Erwartungswert 7 einer Observablen F laBt sich bekanntlich mittels des 
statistischen Operators schreiben: 


(8) F = Spur (WF) = f(t). 


Diese GréBe ist eine vom Bezugssystem in $ unabhangige reine Zeitfunktion 
f(t). Fur die zeitliche Ableitung dieser Funktion ergibt sich (unabhangig von 


(4) Diese Bewegungsgleichung der Dichtematrix entspricht bekanntlich dem Liou- 
villeschen Satz der klassischen Mechanik. 
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einer etwaigen Bewegung des Bezugssystems) 
= 4 OW lo 
FG == Spur i =| Dif Eis w-(F) | 
° | Pq t pal 
td i ) 
= Spur Toi eae —¥ W [HF || + WF) = 


0 
= Spur (WF), 


ne 


(9) ao F | (Ehrenfest) (?). 


3. — Spezielle Bezugssysteme. 


Die bisher angegebenen bezugssystem-unabhangigen Gleichungen sind da- 
durch charakterisiert, daB8 in ihnen keine totalen Ableitungen d/dt von Ope- 
ratoren und Zustandsvektoren vorkommen. Diese Ableitung d/dt einer Grobe 
in $ muB sinngemàB die vollstàndige Anderung dieser GréBe bezogen auf ein 
speziell gewaihltes Bezugssystem in 9 beschreiben. Erst bei Angabe eines 
speziellen Bezugssystems in © ist also die Bedeutung von d/dt definiert und 
umgekehrt laBt sich durch Angabe der Bedeutung von d/dt ein Bezugssystem 
definieren. 

Zum Schrédinger-Bild gelangt man durch die Forderung, daB man fiir die 
totale Zeitableitung eines Operators verlangt 


a ca) 


Man legt also ein Bezugssystem zugrunde, gegeniber dem die Operatoren p 
und g ruhen: 
dsp dsq 
IO 
a dt dt 


Im Heisenberg-Bild definiert man die totale Zeitableitung durch. 


duF 9 
(12) sa 
so daB in diesem Bezugssystem der statistische Operator ruht 
dW 
——_ =O. 
(13) di 


(5) Der Erwartungswert von QW ist also zeitlich konstant (allgemeines « Integral 
der Bewegung »). 


t 
< 
lc) 
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Die bildunabhingigen Bewegungsgleichungen liefern damit die bekannten 


Gleichungen im 


Schrédinger-Bild 


Heisenberg-Bild 


pet: SI (E) a a14(3), 

a F=ias1+ 5, =. 

© - SE ar a) (a 
a i [2290]. =o, 

> (8), af rf), 
ee eee = 0). 


SchhieBlich laBt sich das 
d€ = F°+ FF ausgeht, 
stellung gewinnen: 


Wechselwirkungsbild, das von einer Aufspaltung von 
natùrlich ebenfalls aus der bildunabhangigen Dar- 


Wechselwirkungsbild 


Det.: ra I [eo] +(£ DI 
QF =s9eg) DI, 
(6) sa — [aw]. E 
=~ ix), 
=~ daey 
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| 4. — Der Zeitoperator. 


Zur Beschreibung der Bewegung von Operatoren und Zustandsvektoren 
| in $ wird die Zeit als c-Zahl-Parameter eingefiihrt. Dariiber hinaus muB jedoch 
| der Observablen « Zeit » des quantenmechanischen Systems nach Forderung 1.) 
auch ein hermitischer Operator zugeordnet sein. Wir wollen zeigen, wie man 
i einen solchen Operator © einfiihren kann und welcher Zusammenhang mit 
| der c-Zahl-Zeit t besteht. - 

Zunéichst wollen wir die vereinfachende Annahme machen, da8 der Hamilton- 
| Operator nicht explizit von t abhangt. Aus der klassischen Mechanik ist be- 
i kannt, daB man vermége der Lésung der klassischen Bewegungsgleichungen 
| die Zeit als Funktion von p und q in der Form 


i (14) t —t) = T(p, 9) — TPo; %) 
darstellen kann, wobei 


| G@lgl 
| (15) =, T=1 


| ist. Die Operatorfunktion, die nach den Forderungen 1.) und 2.) der klas- 
| sischen Funktion 7 zugeordnet ist, stellt den Operator der Observablen 
«Zeit » dar 

| (16) C= (p,q), 


| fiir welchen gilt 
| 0 
i (17) Ca 


| Der Zeitoperator T ist explizit von der c-Zahl-Zeit t unabhàngig, wenn dies 
| auch fiir 70 gilt. Die Bedeutung der « Ring-Ableitung » kommt durch Gl. (17) 
besonders deutlich zum Ausdruck. 

Aus Gl. (17) folgt unmittelbar, daB #0 und © zueinander kanonisch konju- 
giert sind, und somit die Unbestimmtheitsrelation zwischen Energie und Zeit 


bo| Ss 


| (18) AE-AT > 


€ 


(6) Das heiBt im Schrodinger-Bild dg /dt=0 und im Heisenberg-Bild da C/dt=—1- 
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Wenden wir uns jetzt dem allgemeinen Fall zu, daB die Hamilton-Funktion 
H auch explizit von t abhangt. Die Funktion 7, die nunmehr auch explizit 
von t abhangig sein muB, gentigt der Gleichung 


aT or 
(19) ar = (7, +(5) =I, 


Pd 


Die Observable « Zeit » wird durch einen Operator 
(20) C= Tp, q, t) 
dargestellt, fiir den jetzt gilt 


0 i OT 
(21) C = [RC] + ea DI] 


Man sieht, daB 70 und © i.a. nicht mehr zueinander kanonisch konjugiert sind. 
Nach dem Ehrenfestschen Theorem [Gl. (9)] 


(22) coon 
erhalt man 


(23) C(t — TC (ty) = + — to. |, 


d.h. der Erwartungswert des Zeitoperators ist also stets die c-Zahl-Zeit t. Dabei 
bedeutet TC (to) = Spur {W(p, q; to) Tp, q; to)}- 
SchlieBlich werde noch die Streuung des Zeitoperators 


(24) AT =V(C 1G) =VE_— eo 
untersucht. Aus 
(25) T= CCLCET HIE 
folgt mit Gl. (23) 

clacio 


so daf fiir die Streuung von © gilt 


(26) AT(t) = AT(t) . 
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Die Streuung des Zeitoperators ist also bei reversiblen Bewegungen stets zeit- 
lich konstant. 
Fiir den zu 3(p, g, t) konjugierten Operator A(p, q, t) hat man allgemein 


ldie Bedingung 


7 h 
(27) [RL = sl. 
{Es gilt also 
0 
1(28) K=1+(S) 
Ct /p.q 
jund damit fiir den Erwartungswert 
i 
: ponte (AR 
(29) R(t) — R(t) = t—th+ Di dt’. 
ot Bed. 


In den Erwartungswert von A geht also auver ¢ auch noch seine explizite 
i Zeitabhingigkeit ein. Im zeitunabhargigen Fall fallt diese fort und es ist 


| (30) R(p.G) = T(P, 9). 


5, — Diskussion. 


Die angestellte Analyse zeigt, daB® die Zeit in der Quantenmechanik eine 
Doppelrolle spielt. Insoweit sie einen reinen Ordnungsparameter darstellt, der 
|a priori zur Beschreibung des topologischen Nacheinander vorgegeben ist, kann 
| sie keine (streuende) Observable sein und mu8 daher als e-Zahl t in die Theorie 
eingehen. Andererseits ist die Zeit eines quantenmechanischen Systems, die 
inach (16) durch das dynamische Verhalten des Systems selbst gegeben ist, auch 
jeine beobachtbare physikalische Grobe, also eine Observable, die den Quan- 
| tisierungsforderungen gemà8 durch einen hermitischen Operator © beschrieben 
| wird. Wéahrend diese zweite Figenschaft — von der man aber meist nicht 
etisch unproblematisch ist, scheint die Existenz eines 
hst im Rahmen der Quantenmechanik unver- 
aB auch makroskopische 


| spricht — quantentheor 
scharfen Zeitparameters zunac 
| stindlich zu sein (7). Hierauf ist jedoch zu sagen, d 
n denen makroskopische Zeitmessungen méglich sind, einer quanten- 
Die Messung der Observablen 
Festlegung des Zeit- 


) Systeme, a 
) mechanischen Beschreibung zuginglich sind. 
«Zeit » an solchen « Uhren » liefert eine physikalische 


(7) Vgl. E. SCHRODINGER: Berl. Ber., 238 (1931). 
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parameters ¢, wodurch seine Existenz auch im Rahmen der Quantenmechanik 
legitim wird. Die Forderung nach dem Vorhandensein des Begrifts « Nach- 
einander » in der Quantenmechanik fihrt also iiber die Observable « Zeit » zur 
Moglichkeit einer Bestimmung der Zeitordnung. 

Die Tatsache, daB der Erwartungswert des Zeitoperators irgendeines Sy- 
stems stets dem systemunhabhiangigen Zeitordnungsparameter ¢ gleich ist, 
kénnte man dazu ausnitzen, die c-Zahl t allgemein durch den Erwartungs- 
wert des Zeitoperators sogar eines beliebigen Systems zu ersetzen. 

Explizite Zeitabhingigkeiten sind nur durch die Betrachtung von Ver- 
anderungen in einer vorgegebenen Zeitordnung definierbar. Demnach kann 
hier nur der c-Zahl-Parameter t eingehen. Hier liegt der Ansatzpunkt fir die 
Ableitung der quantenmechanischen Bewegungsgleichungen, die allein aus der 
Grundforderung der korrespondenzmàBigen Ubersetzung bildunabhingig még- 
lich ist. Nur die expliziten Zeitabhangigkeiten haben hierbei demnach eine 
tiefergehende Bedeutung. Da die statistische Deutung der Quantenmechanik 
eine andauernd eindeutige Beschreibbarkeit der physikalischen Zustànde er- 
fordert (vgl. (4) und (5)), laBt sich rein deduktiv auch die Bewegungsgleichung (7) 
des Zustandsvektors bildunabhangig gewinnen. Eine weitere direkte Folgerung 
ist die ganz allgemeine Gultigkeit des Ehrenfestschen Theorems (9). — Totale 
Zeitabhangigkeiten von der c-Zahl-Zeit sind nur im Zusammenhang mit einem 
willkirlich wahlbaren Bezugssystem im Hilbert-Raum sinnvoll (z.B. Schré- 
dinger- oder Heisenberg-Bild). 

Allein mit Hilfe der Observablen « Zeit » làBt sich fiir den Fall von Syste- 
men, deren Hamilton-Operator und damit auch deren Zeitoperator explizit 
von der Zeit unabhangig ist, aus der Eigenschaft (17) (8) die Unbestimmt- 
heitsrelation zwischen Energie und Zeit begriinden. Fir Systeme mit explizit 
zeitabhangigem Hamilton-Operator ist der Zeitoperator Tim allgemeinen selbst 
explizit zeitabhàngig. Da hier eine direkte Berechung von © analog zu (14) 
und (16) nicht méglich ist und die Eigenschaft (21) allein © nicht vollstandig 
bestimmt, ist der Zeitoperator in diesem Fall unendlich vieldeutig. Der Grund 
dafùr liegt in der Tatsache, daB ein explizit zeitabhingiger Hamilton-Operator | 
ein System beschreibt, das entweder in Wechselwirkung mit einem makrosko- 
pischen System steht oder auf ein bewegtes Koordinatensystem bezogen ist, 
so daf in den Zeitoperator ein in weiten Grenzen willkiirlicher Anteil der Zeit- 
skala des makroskopischen Systems bzw. der vorgeschriebenen Bewegung des 
Koordinatensystems eingeht. 

(8) Die von W. PauLI (Handb. d. Phys. (S. Fliigge), Bd. V/1; S. 60) gegen die 
Beziehung (17) vorgebrachten Bedenken sind nach unserer Ansicht nicht stichhaltig, 
weil ein in der Argumentation beniitzter Satz von Dirac (Quantenmechanik, S. 56) iber 
die Eigenwerte konjugierter Operatoren nicht zu gelten braucht, wenn die beim Beweis 


verwendete unitàre Transformation keinen Automorphismus zwischen den ursprung- 
lichen und den transformierten Eigenfunktionen der Operatoren liefert. 
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sive interpretation. 2. Formal definition of symmetry. — 3. Consequences 
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stricted distribution law. Linearly independent elements. 6. The ope- 
rator function. 7. Superselection subspaces. 8. Continuity. 9. Re- 
duction to unitary and antiunitary transformations 0. 10. Dependence 
on the group element g. Representation up to a factor. 


Remark. 


This note does not contain any new result. But the proof that symmetries 
of quantum mechanical systems give rise to either unitary or antiunitary 
representations takes a new and simple form. In order to state the problem, 
it seemed, however, necessary to discuss the implications of symmetries in 


some detail. 


4. — Introduction. 


11. Complete set of observables. — We assume a definite physical system S 
to be given. The system will be described in a Hilbert space #. Take any 
observable (self-adjoint operator) A’ and construct a set % in the following 


A'EA, 
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way Ae A if and only if A commutes with all other elements of A. Now 


if AeA; A'eA then 
AA'eA 


aA + bA'EA, (a, b real). 
Hence from the distributive and associative law it follows: 
2 is a maximal abelian ring with real coefficients. 


Let BEM. This B defines in the same way an abelian ring 8. In this way 
we may construct a whole set of such rings A, B, €, D, ... (*). The set of 
all observables is not a ring since the product of two self-adjoint operators 
need not be self-adjoint. All operators of any such ring can be brought simul- 
taneously into diagonal form (spectral representation) and hence can be 
measured without interference. There is a proof by v. NEUMANN (!) that in 
such a ring X there always exists an element A, say, such that all elements 
of the ring are functions of A,, namely A,;=f,(A,). In principle the measu- 
rement of A, alone would determine completely the state of the system. This 
might be technically impossible. In fact, normally one takes one convenient 
operator A,, say. If now all other elements of % cannot be represented as 
functions of A,, one takes another one A, and so on, till one has as many ope- 
rators as are necessary in order to generate the whole ring A. The set of ele- 
ments A, ... A, is called a complete system of commuting observables. In 
actual cases, one frequently seems to be able to construct such a system. For 
every ring XU, 8, ©, D... we assume such a set of generating elements: A,, 
B;... and call the set of rings {X, 8, € ...} complete, if every observable belongs 
to at least one of the rings. 


12. Superselection rules. — Obviously the multiples of the unit operator 
commute with all our observables. If there exists, apart from that, another 
operator X’ commuting with all observables, then the Hilbert space can be 
decomposed into superselection subspaces. No observable will have matrix 
elements between different superselection subspaces (that is the reason for 
the name) (2). 


Proof: Together with X', X'* (the hermitian conjugate of X) also com- 
mutes, hence X= 3(X'+X'*) is hermitian and commutes also. We assume 
now that XY has a discrete spectral representation: 


(*) These rings overlap of course. 
(!) J. v. NEUMANN: Math. Grundlagen der Quantenmechanik (Berlin 1932). 
(?) C. G. Wick, A. 8. Wicurman and E. P. WiGNER: Phys. Rev., 88, 101 (1952). 
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with 
xX Sa; DCR) (Ak) and X}=—X,, 
X, projects into that subspace 3, of 0, where X has the eigenvalue 4,. 
(Always 4;4A,, if 2; occurs n; times, Fe; has dimension n; < oo). 
From 
[XA] = YA:[X;4}=0 


follows by multiplication with X, to left and right 


Let pic k;, Pre Hy, t#AK. Then at least one of A, A, is #0, Say 
Un #0. Thus 


Pi, Ax) = X iPir AXxPx) = CX iPi, X Ax» = XxX is Apr) = 0). 


Since A can be any observable, the proof is complete. Any observable 
can then be brought to the following form: 


and the dimensions (in general infinite) of these boxes on the diagonal will 
be the same for all observables. Examples of such superselection rules are 
the charge conservation and baryon conservation. 

One sees that in the presence of superselection rules arbitrary unitary 
transformations cannot be admitted, since they would in general mix up the 
superselection subspaces. Those unitary transformations, however, which trans- 
form from a basis {y,} belonging to the ring A to a basis {p;} belonging to the 
ring B, leave the subspace invariant and are admitted. We may even admit 
those transformations which furthermore induce a permutation of the super- 


selection subspaces. 


2. — Definition of symmetry. 


A transformation group @ may be defined in many Ways, @.9.: 
— Interchanging of particles, 

— Charge conjugation (better perhaps: matter conjugation), 
— Lorentz transformations, 


— Gauge transformations, ete. 
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21. Active and passive interpretation. — There are two ways of interpre- 


tation: 
i) passive, 
ii) active. 


Example: The three dimensional rotation may be interpreted passively 
as a rotation of the co-ordinate system or actively as a bodily rotation of the. 
system S under observation. 

In what follows we shall adopt throughout the active interpretation, which 
is always feasible, whereas the passive one looks very strange, if not impos- 
sible in some cases. 


22. Formal definition of symmetry. — Now we define what we mean by 
saying that our physical system S is invariant under a group G, which we 
shall then call a symmetry group of S. 

A material system $ is said to have a symmetry group G, if for every 
geG there exists 


i) another material system S’= gS (symbolically), 


ii) a definite function F, for all observables 
such that 


a) A'= F(A) is again an observable, and 


B) the possible results of measuring A’ in S’ are the same as for measur- 
ing A in S and both have the same probability distribution. 


In order to put this physical definition in a more mathematical form, we 
introduce a special notation: 

A physical state does not determine an element of #, since together with 
g also mp with |w|=1 represents the same physical state. 

For reasons of clarity, we shall denote by: 


f, 9, h ... elements of 2, 


®, Y, P state vectors, 7.e. elements with norm 1, 
(d1) 

d, È, @ unit-rays, representing physical states, 

Î, 9, h rays. 


A ray f is the set of all elements Af with 0 <|A|< co (f fixed). A unit-ray O) 
is the set of all state vectors mp with |@|= 1 (p fixed). Rays are not vectors. 
Addition, etc., are not defined. 
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We now give the formal definition of s try: 
g symmetry : 


If we denote the state vectors of S by y, y,... and those of S’ 
by 9’, y', ..., then there is a one-to-one correspondence of the phy- 
sical states (rerpesented by unit-rays) 


) poker 
such that the probabilities are conserved: 
<P, POI =1@, w) | 
if g'e@, y'e d' and PED, pep. 


We shall quote this definition henceforward as (D). 


3. — Consequences of the existence of a symmetry group. 


3:1. Ewistence of a unit-ray transformation. — It follows from (D) that the 
existence of a symmetry group @ implies only the existence of a group T of 
unit-ray transformations §, which maps the set of all unit-rays onto itself 


P—7'= 09. 


That this transformation group is uniquely defined by @ and is isomorphic 
to @, follows from physical experience. In this form, however, it is useless 
for quantum mechanics, since the superposition principle is true for state 
vectors rather than for unit-rays. The problem is therefore to find a group 
of transformations 0 defined for elements of the whole Hilbert space, which 
is homomorphic, or even possibly isomorphic to the group 6, and as simple 
as possible. 


312. Class of transformations 0 leading to the same 9. — Let T be a group of 


unit-ray transformations. If 6¢7, we consider all transformations 0 of # 
which correspond to 6. Correspondence is defined by 


6~6 it 

9pe dg for all ped. 

Consider the set of all non-singular 0, which correspond to all Ò e È. These 0 
form a group 7, which is homomorphic to the group T. This. is a very large 


homomorphism, since many 0 correspond to the same 0; for instance 0, and 0, 
may be defined completely different ly for all elements with norm #1 and yet 


~ 
Uni 
Usi 
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correspond to the same 0. Consider the unit element of 7, 6= E. The set of 
transformations 0 ~ E forms an invariant sub-group 7, of 7 (the kernel of the 
homomorphism). The factor group 7/7) is then isomorphic to IS Go 0, then 
the coset 67, = 7,0 contains all transformations which correspond to 6; we 
write 07,~6. The transformations 0 e 7, have the property of multiplying 
state vectors g only by unimodular complex numbers (leaving their «length » 
and «direction » constant). On other elements with norm #1, they may be 
defined arbitrarily. 

Any coset of 7, may be represented by one of its elements, but if we select: 
one element of each coset, then in general these elements will not form a group. 

It will be shown that elements 0 can be selected from the cosets such that. 
they are either unitary or antiunitary and form a representation up to a 
factor. 


33. Conservation of the norm. — First we postulate 


iù 


By this 0 is defined in the whole of # once it is defined for state vectors. These 
0 still form a group. The invariant subgroup 7, now contains all those 6 
which merely multiply each element of A with a unimodular complex factor 
such that 0f=©,,'f: The factor depends only on {/|f|. We shall hence- 
forward call this invariant sub-group 2. This is in fact a very strong re- 
striction. Its effect is that two transformations 0 and 6’, which both cor- 
respond to one single 0, differ at most by a 00€ 2. 

Obviously the restriction is not too strong, since it does not affect the 
transformation of state vectors. 


3°4. Complete orthonormal systems. 


Theorem 1: Complete orthonormal systems are transformed into com- 
plete orthonormal systems. 


Proof: Let {y,} be a complete orthonormal system, and 


then 00'~E. y;—>y,=6y; and from (D) 


(3.2) |<Oyi, Ove» | =1<wi, va) | = da - 
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It remains to show that {yt is complete. Assume it is not; then there exists 
I 

a y, orthogonal to all y,. Then d'w=%, would be orthogonal to all y; 

against the assumption. 
3°5. The restricted distribution law. Linearly independent elements. 
Theorem 2: 


If f, g, h... are linearly independent then 0f, 6g, 0h... are also. 


i) 
ii) O(f+g9) = off, g)Of + (9, f)Og with |w(f,g)|=1 for any pair of ele- 
ments f, g. 


Proof: 
i) By induction. Let it be true for n vectors Pec 
ifn Span a subspace M,, 
Of, ... Of, span a subspace Mî. 


Let f,.,=}f be linearly independent of f, --»fns Then there must exist 
a h such that 


Ch, f= 9; i=1... N, 
<h, {>} #0. 

It follows that 
<Oh, 0/9) = 0 ; A 
<0f, Of) #0. 


Therefore Of is linearly independent of Of, ... 0f,. Since the Theorem is 
true for n=1, it holds for all n. 
ii) For f=Ag it is trivial. Let f, g be linearly independent. They span 
a sub-space M. For any he M~ 


|<oh, 0f+g9>|=|Mhf+@®|= 9. 
Hence 0(f-+g) belongs to the sub-space spanned by 0f and 0g: 


O(f +9) = Af, Of + nf, 9)09 - 


Define 
G) PD 


Si —— g; <k gr a es 
lg} °' i 


k(f, 9) = f 
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Then 
|<6k, Of +9)>|=|<k, f+ 9\/=1% DI] 
but also 


Hence |A(f,g)|=1. From interchanging f and g it follows that A(f,g) = 
= Wg, f) = olf, 9), q-0.d. 
In particular: 


(3.3) 0Xow:=Ye0v, with  |e,|=|e,|. 


3°6. The operator function. - We can now write down explicitly the ope- 
rator function 
(A T= RGA 


g 


Let 6~ @(g). If 


A= > a; | pi) <p: | , 
then 


(3.4) F(A) = A’= > a,\6y><6y.| = > a,pp><yi| - 


We have to show that this is unique. Thus any 0~ 6(g) must give the same 
result. Take 0/'— 0,0; 062: 


A" = > a;|6'p><O'yi| = > €: OP ><oy,| => 10,0; VAICA =A’. 
It follows that A’ has the same eigenvalues as A and the transformed state 
vectors y, are the new eigenvectors. 
Theorem 3: The expectation values are conserved. 
Proof: 


= Leys  0p=Yc0,, 
then 


<Op, A'bp> = Y eFa,0,<Op;|Op><Oy. |Oy.> = Y a,|e,[, 
(p, Ap) = Dale |. 
These are equal on account of (3.3). This means: Corresponding measurements 


in corresponding states will always yield the same results. Hence the whole 
physics is unchanged. That is the meaning of our definition of symmetry. 


3'7. Superselection subspaces. 


Theorem 4: All 6€7 will, apart from a mapping of superselection 
subspaces onto themselves, at most induce a permutation of these spaces. 
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Root 


i) From Theorem 1 it follows that orthogonal subspaces are transformed 
into orthogonal subspaces. 


2 


ii) From Theorem 3 it follows that the transformed observables have 
no matrix elements between the transformed superselection subspaces. 


iii) The transformation of the observables is a mapping of the set of 
of all observables onto itself. This set defines the decomposition dt = 
- FK, K,© #,@... up to the order in which the #; are counted, thus a 
a certain 0e 7 induces a well-defined permutation of the #; (which might 
be the unit-permutation). 


3°8. Continuity. — There are two types of continuity : 


i) By symmetry operations on the physical system neighbouring physical 
states are transformed into neighbouring physical states. This is a physical 
experience. It is contained in our definition of symmetry (D): 


|<y’; PY |= |<y, »? le 


Let g = p(t) follow a continuous path, then both sides have to be continuous 
functions of t. This does not imply that g'(t) is also continuous, since it may 
contain any discontinuous phase factor a(t). 
ii) The symmetry group G may be continuous or partly continuous: 
There may exist continuous paths g(t). 
We shall presently discuss the first type of continuity. For this we define 
distance and neighbourhood of unit-rays: 


(43) The distance of two unit-rays @ and @ is d= min, |p — yp|| where 
ped and ped is an otherwise arbitrary pair of state vectors. 


(d4) A e-neighbourhood of @ is the set of all unit-rays @ whose distance 
from @ is < e. 


With these definitions the first type of continuity can be stated as. 


Theorem 5: The transformations 0, which correspond to elements g e @ 
are continuous. 

This is obvious since neighbouring rays are transformed in neighbouring 
rays. 

Theorem 6: In the set of all 6 ~ È there is at least one which is con- 
tinuous in the whole considered superselection subspace. 


6 - Supplemento al Nuovo Cimento. 
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Sketch of a proof: By construction. Take an arbitrary state vector 
vo. It belongs to $. Take now an arbitrary state vector p, € 9’, where 0G = P'. 
Define 0g) = Ag, for all 2 (0 <|A|< 00). Hereby 0 is defined for multiples 
of go. We now define it for any element g: 

Let p be any state vector and # its unit-ray. Find © such that: 


|g.—oy| = min; (y= oy is nearest to go). 
Find further © such that, with y'e0$ 
lo — 0‘9'|= min; (y,=o'y' is nearest to 9). 


Define 0uy)= wy, for any 4 (0 <|u|< 00). We have to show that this 0 is 
continuous. Given f, we must find a neighbourhood N, of f such that for 
all ge N, |0f — 0gjl< e. 

Put f=w, g=uD, where w and ©, are nearest to m in the above 
sense. Then 


||@f — og =| Av. — abi <IAl- vo — Poll + 1A — #1]. 
To make this smaller than e, choose 


|\A |< 6/2, 
7 I E 

Do — ca 

To achieve the last inequality, take the set of all D, fulfilling it. To each of 
them belongs a unit-ray D'. Call the set of these unit-rays {D'};. Transform 
each of these unit-rays with 07!. This gives a set {D};: Find the nearest 
element ®, to gy of each De {D};: These D, form an open set containing yp, 
and defining a neighbourhood of wi: Any g=u®D, for which ©, lies in this 
neighbourhood and for which |2—w|< e/2, fulfils the condition. 


3°9. Reduction to unitary and antiunitary transformations 0. 


Theorem 7: Given any 0, we can choose such a 0-0 that 0 is either 
unitary or antiunitary in the whole considered superselection subspace. 


This theorem was proved by E. P. WIGNER (*) without taking into account 
the antiunitary operators. This has been done explicitly also by E. P. WIGNER 


(*) E. P. WIGNER: Gruppentheorie und ihre Anwendungen auf die Quantenmechanik 
(Brunswick, 1931), p. 251. 
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in his Lorentz-chair lectures, University of Leyden 1957 (not published). His 
proof is rather complicated and the following simple proof is the only essen- 
tially new thing in this paper. 


Proof: The theorem means that given a 0 and a 0’~ 6, we can find such 
a (€ 2, that 0.0’= 0 is either unitary or antiunitary. We construct 6). 


a) Construct a continuous 0’~ È. following the lines of the proof of 
Theorem 6. From Theorem 2 it follows that for linearly independent ele- 
ments f, g, h 


(ft h+g) = olf +h, g)[alf, h)0'f + o(h, POR] + og, f+ h)0'g - 
Interchanging g and h and comparing the coefficients leads to 


olf +g, half, 9) = of +h, 9)0(f, h) coefficient of 0'f, 
o(f+g, h)w(g, f) = og, f+h) coefficient of 0’g, 


o(f +h, g)a(h, f) = w(h, f+) coefficient of 0/h. 


This gives, by dividing the first two and eliminating o(f--h, g) with the help 
of the third equation 


(3.5) olf, g) = 


obs f+9) og HON 
oth) alg, f+) 


From 6/(f-+0) = w(f, 0)0'f = 6'f follows w(f,0)=1. We put (0, f) = 1/u(f) (*). 
Now letting h—0 one obtains 


u(f) 


olf, 9) = uf +9) . 


u(f) is unimodular by definition. That it is continuous and independent of 
the norm |||, follows from the corresponding properties of 0'. 
We have therefore a unimodular function «(f) such that 


uf +9)0"(f+9) = uNOF+ (909. 


u(f) is an element 0,62. Now 6,0'=6 is distributive: 0(f+9) = 6f+0g, and 
0 is again continuous. 


BR) (A) <Of +9), OF +9)> = +9, +9 = Ifll+|g| +2 Re 49, 
LO(f +9), Ff + 9)> =|lf+||g]| + 2 Re <7, 09) . 


(4) The following argument is due to Dr. J. M. JAUCH (private communication), 
(*) This is not as obvious as it first seems. See the discussion in the appendix, p. 85. 
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Hence 
Re <0f, 0g) = Re <j, g 


and from (D) 


|K0f, 09 |=|<F @|. 
Only two solutions exist: Im <0f, 09) = + Im <fg> thus: 
Either CG. Ug> rg 0 is unitary, 
or «Of, 09% =<9, f 0 is antiunitary 


Since 0 is continuous, it is either unitary in the whole subspace or antiunitary 
in the whole subspace. Note that the same cannot be concluded for the 
whole of dC. 


3°10. Dependence on the group element g. Representation up to a factor. — 
We now have unitary or antiunitary 6’s. For this we have paid heavily. 
Beginning with a given 0’, we found the unimodular «(f) uniquely defined up 
to a unimodular factor @(0), which no longer depends on f. 

Assume we had started with a 0”, we would have found a «"(f) represented 
by the group element 07€ 2. Now 6,6’ and 00” can differ only by another 
90€. But this 07 can only be a constant unimodular ©, since otherwise 6, 
and 6, were not uniquely defined (apart from a constant factor). 

Therefore we can postulate that 0 shall be unitary or antiunitary (we cannot 
choose between these two) but we pay for this by loosing all freedom but for 
the freedom of multiplying 6 by a unimodular (6). Apart from that, 0 is 
then a uniquely defined function of È and therefore of the group element g. 

With a choice of the remaining free function ©(0) everything is fixed. We 
have then chosen exactly one transformation 6(g) from every coset. We cannot 
expect these transformations to form a group, but we have always 


(91) O(g2) = ©(9192) 9 (Gre) - 


This is called a «representation up to a factor ». 

It should be clear that if we admit arbitrary phase factors, then the 0 form 
a group. Its invariant sub-group corresponding to = is now the group @ 
of complex unimodular numbers, and the factor group 7'/@ is still isomorphic 
to the group G. It is by fixing these phase factors that the group property 
is lost. It may happen, however, that these phase factors can be chosen such 
that the 6’s still form a group. Under which circumstances this is possible, 
is discussed in full generality by V. BARGMANN (Ann. o. Math., 59, 1 (1954)). 
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APPENDIX 


= 


In the proof of Theorem 7 we put (0, f)=1/u(f). 
One may wonder whether (0, f) makes sense, since from 


AFL) = off, M}04+ olh, {)0" 


no information is obtained regarding (0, f) and even the lim (h, f) cannot 
be found from it. se 
That such a limit exists if h — 0 in a well defined way, is shown as follows: 
From the proof of Theorem 6 it follows that 0'Af=76'f for any A. From 
Theorem 2 ii) it follows that 


O' (f+ Af) = (1+ 4)0'F = [o(f, Af) + 20047, NIG, 


hence 


1— off, Af) 
1— ©(4f, f) ” 


A= 


Since |w|=1, the numerator and denominator of the ratio lie both on 
a unit circle with centre at +1 in the complex plane. In general the ratio 
and with it 2 will be complex. Real / are possible only for @(f, Af) = olAf, f) 
for ~1 this gives A——1. Other real 4 are possible only for ©(f, Af) = 
= (Af, f)=1. This is also a trivial solution for complex 2. Any solution 
for complex A must coincide with this solution when A becomes real. 

From (3.5) follows that the limit 


oh, f+9) 


ie oh, f) = Of, 9) 


exists independently of how h—0. Since g is arbitrary, f+g is arbitrary. 
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Therefore with f/+g=k 


o(h, k 
m Ma, exists independentlu of how h>Q0. 
h>0 (Ah, Î) 


Choose now h=A-ky); where k,=%/|K| and take 7 real. Then 
(Aki, ky) = 1 
for all 4> 0 and it follows that 


o(h, k) 


1 
=o o(h, ‘anes REA (Ao, i) =U, (f) 


exists. Since k, is arbitrary, we conclude that 


We have then 


w(}:g) = lim no) = iat ae. ee i Uno (7) : 
di h_->0 oh, Î) We. ah, f) Un (f+ 9) 


Since this is independent of h,, it follows by comparing it with another 
limit using kj, that 

Uy, (fu, (k) 

Un: (7) E Upi (k) 


=... = F(h, hh). 


Therefore uw, (f) and u,(f) differ only by a unimodular factor which does 


not depend on f. Since this is irrelevant for the further arguments, we may 
simply write w(0, f)=1/u(f) as was done in the proof of Theorem 7. 

One may wonder how 6'Aj=A06'f for arbitrary complex / could be true, 
since we finally arrive at a 0 which may be anti-linear. But obviously this 
property has gone then: 


OAf = w(Af)O'Af =A-u(Af)O'f and it may be that wu(Af)~u/(f). 
In a way the property of 0’: 


OAL = ABT 
and the property of 0: 


O(f+ 9g) = Of + 0g 


are opposite extremes, which coincide only for unitary 0. 
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Two Theorems on Scattering (*). 
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University of Oklahoma - Norman, Okla. 


(ricevuto il 19 Gennaio 1959) 


CoxtENTS. — Introduction. — 1. The evanescence of wave packets. 
2. Uniqueness of the free hamiltonians of a scattering system. 


Introduction. 


In this paper we shall prove two theorems which belong to the still growing 
body of literature on the formal theory of scattering (1). Any wave packet 
representing a free particle develops in time in such a manner that the proba- 
bility of finding the particle in any finite region of space vanishes in the limit 
as t>-+oo (Theorem 1). 2) In a multi-channel scattering process the set 
of (free) Hamiltonians defining the asymptotic motion of the fragments is 
unique (Theorem 2). Although Theorem 1, called for obvious reasons the 
theorem on the evanescence of wave packets, is well known, a rigorous proof (1) 
does not appear to be in the literature. In Theorem 2 we assume the appro- 
priateness of Jauch’s definition (?) of a multichannel scattering process, a 
definition which is an extension of his previous definition (?) of a single scat- 


(*) Supported by the National Science Foundation. 

(1) L. I. Scurr: Quantum Mechanics (New York, 1949). A particular case of the 
Theorem I referring to a one-dimensional minimum wave packet is considered on 
pp. 54-59. 

(2) J. M. JaucH: Helv. Phys. Acta (in print). 

(3) J. M. Javon: Helv. Phys. Acta, 31, 127 (1958). 
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tering process. Although these two theorems appear to be unrelated, it is 
convenient to consider them together, since the proofs are based on very 
similar mathematical results. This similarity is due to the fact that both 
theorems are formulated in terms of essentially the same kinetic energy ope- 
rator, cP?, where c?>0. 

The proofs of the theorems are based on the theory of linear operators 
defined on a Hilbert space # of (square integrable or L?) state functions. 
The delta function and non-normalizable functions are thereby completely 
avoided (4). The principal symbols used are: || ||, norm of an element in # 
and of a bounded operator (defined everywhere in #); —, strong convergence; 
— weak convergence; and (5) “+, convergence in the usual sense of point-wise 
convergence. 


1. — The evanescence of wave packets. 


It is perhaps physically reasonable to expect that a wave packet which 
represents a free particle will develop in such a way that the part of the packet 
which differs significantly from zero at #=0 will in the course of time fill 
all space. Because the normalization of the packet is conserved, it is clear 
that the probability of finding the associated particle in any finite portion 
of space should converge to zero as t > co. Before stating Theorem 1 which 
justifies this expectation, we introduce the following notation. 

Let # be an abstract space; #0(X), a realization of FH formed by the set 
of L?-functions, f(X), defined over the infinite three-dimensional Euclidean 
space E, of vectors X= (x, y, 2); and A(K), a realization of # formed by 
the Fourier transforms (also a set of Z?-functions) of the elements of F(X). 
If the abstract element f e # is represented by j(X) e A(X) and by f(X) e #(K), 
then the momentum operator P (with 7#= 1) has a simple representation (5) 
in 90(K): 


(1) (P?f)(K) = K? {(K) 


(*) An outline of the mathematical theory will be found in reference (3). More 
detailed expositions will be found in B. von Sz.-Nagy: Ergebnisse der Mathematik 
und Ihrer Grenzgebiete, 5, 5 (1942); F. Riesz and B. von Sz.-Naqy: Functional Analysis 
(New York, 1955). 

(°) This symbol is not standard and is used in this paper for the sake of clarity. 
Normally the symbol for strong convergence is also used for point-wise convergence 
because its meaning is rarely ambiguous in context. 

(5) J. von Neumann: Math. Ann., 104, 570 (1931). It may not be superfluous 


to point out here that the terms wave packet and element of a Hilbert space are 
synonymous. 
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provided [ K®{(K)} is integrable, è.e., provided K?j(K) e A(K). Also, if t is 
a real variable and ©(t) is the limit of a sequence of polynomials, then 


(2) (@(P?)f) (K) = o(K°)f(K), 


provided again @(K?)f(K)e #(K). 

Let A be an arbitrary Borel set of points of E, with finite Lebesgue mea- 
sure u(A). Also, let E, be a projection operator, which when applied to 
_f(X) e #(X), has the following property: 


OA 
(3) (2, f)(X) = 


On otherwise. 


Therefore, if a particle is represented by the normalized wave packet 
f;= vv at time i, the probability of finding it in A at time t is 
|E,wde/|vi|, or, if U; is a unitary group and y= U,y, then this probability 
is given by |£,U.y|/|¥| 
in the following way. 


2, The theorem to be proved may now be stated 


THEOREM I. Let y be an arbitrary element of 36; U,=exp[— tH yt]; 
H,=cP?, where c?>0. Then |£,U,y|| > 0 as t > = co. Although it is as- 
sumed that the non-relativistic Hamiltonian describes the free motion of the 
particle whose wave packet at t= 0 is y(X), it will be clear that the proof 
can readily be extended so that the theorem refers to a wider class of Hamil- 
tonians, one which includes the relativistic Hamiltonian V/m? + P*, for example. 

The theorem will be proved with the aid of six lemmas. The first four 
establish the validity of the theorem for a set of elements dense in 30. With 
the aid of the last two it is shown that the theorem holds for all elements 
of FC, i.e., for every wave packet. Since the proof of the first lemma is readily 
available, it is stated without proof. 


Lemma 1. (Wiener’s Displacement Theorem) (7). If g(X) belongs to L° 
‘and if almost everywhere its Fourier transform differs from zero, then if p(X) 
belongs to L? and e > 0, there exists an integer N together with a set of real 
vectors A, and complex numbers A, (n=1, 2, ..., N) such that 


III PX) FA, p(X — Bs) PRO 


n=1 


where dX= de dy de. In this form, Lemma 1 is a slight generalization of the 
displacement theorem as stated in reference (7), where the theorem refers to 


(7) N. WIENER: The Fourier Integral (Cambridge, 1933), p. 100. 
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functions defined on the line. However, it is clear that by considering the 
closed linear manifold formed by products of three (or more) functions, each 
defined on the line, Wiener’s theorem must hold also for L?-functions defined 
on a space of three (or more) dimensions. 

For the generator y of the displacement theorem we may use g(K)= 
= exp [— a?K?] with a? > 0 and fixed, since (1) its Fourier transform vanishes 
nowhere in E; and (2) 


lp|P =[[Jexr [— 2a K°] AK A (a), 


1.6., ped. It we. define 
(4) Pr (K)= ~(K — Ky) = exp[—aK—K)?], 


we readily deduce that gx ¢d and that the Fourier ‘transform of U,yg )(K) is 


(5) (U @x.)(X) = (ba)! exp [—a?(K — K,)? —ictK? + iK-X]aK = 
= [ 2(a? ole ict) |-?exp {— iK,| 5 (ctK, — X)] exp | (ctK, "75 2 "| 
a? + act 


and that 


12) | (Uipx,) (X) |? =[4(a4+ c2t2)]-® exp 


2a? (4X — ctK,)? 
at+ ce?t? : 


This result leads to 


Lemma 2 (8). |(Uipg,)(X) P< [4 (a4 + ct]. 
This lemma together with (3) implies that 


(7) |E1Ugr, 


= [(T,ge) IATA 0a | dX= p(4)[4at+ 2 8)}-t. 


«/ 


A 


Therefore we obtain 


Lemma 3. |H,U.p,|| 0 as t + + co (even uniformly in K,). 
The result is exactly Theorem 1 when the latter is applied to the set of 


(°) This result was obtained in another context by J. M. Cook: Journ. Math. 
Phys., 36, 82 (1957). 
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elements {px} in which all vectors K, are represented. In addition, since the 
operators {E,U,} are linear and the manifold D formed from the set {gx} is 

0 
dense in # according to Lemma 1, we have verified 


Lemma 4. A linear manifold D dense in & exists such that for all 
feD, | E,U.f|>0 as t > + 00. 

If an arbitrary sequence {tn} is selected from the interval 0 <t< co so 
that t, > oo as n+ oo, or if a similar sequence is selected from the interval 
— co<t<0 so that t, > — co as n > co, then Lemma 4 implies the 


Corollary. For elements in the domain of Lemma 4, |£,U, f+ 0 if 


t, > +1 co aS nN > 00. 


To complete the proof of the main theorem we must show that Lemma 4 
is valid for all elements in #, i.e., for each element in D, the closure of D. 
This extension of Lemma 4 can be carried out with the aid of the following 
two lemmas. The first of these is the Banach-Steinhaus Theorem (°). Since 
this theorem is one of the basic theorems in the theory of linear operators, 
it will not only be stated here in a form useful for the present demonstration, 
but a relatively transparent proof will also be given. 


Lemma 5. Let {A,} be a sequence of linear bounded operators such that 
1) ||A,|| <M for al» 
2) A,f converges strongly as n + oo for each f in a domain D, which 
is dense in €. 
Then A4,-> A on all #, where A is a bounded operator and |A|<M. 
Proof. Since the Lemma is trivial if M= 0, we shall assume that M > 0. 


1) To prove that {A,} converges strongly in #, it must be shown that 
given any e > 0 there exists an N such that for any ge dé 


|Amg —Amgl| < €; m, m'>N. 


Since D is dense in FH, corresponding to any g E #, there exists a sequence 
{9,} in D such that g, > g as n oo. Therefore 


| Ang —Aw 9l| =|(An Am) Ga—9 + 9a) < 


= | (An —Aw)(Gn— 9)| gu | LA A ou) Grill < 


< 2M||gn = g| da | (Am Am)In 


(®) E. HILLE: Functional Analysis and Semi-Groups, vol. 31, in Amer. Math. Soc. 
Collog. Pub., Amer. Math. Soc. (New York, 1948). 
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Since g,—>g, N, exists such that |g,—g|<«/4M for n>WN,. Also, since 
gnED and {A,,g,} converges strongly as m + co by hypothesis, N. exists 


such that (An Am)gx|| </2 for m, m'>N,. Therefore 
5 € 
[Ang Am gl < 2M ta I 


for m, m'>Max (N, N.), and {A,g} is a Cauchy sequence for arbitrary ge FH. 
Since # is complete, {A,g} converges to an element in #0 as n— oo, an ele- 
ment which may be defined as Ag. This proves the first part of the lemma: 


A, => A om all # as n co. 


2) To show that A is bounded with bound M, we notice that, since . 


A, >A aS n-+ co, corresponding to any £>0 and any geF, N exists 
such that 


[Ag] =|(4—4, + Agg <|A- Ag] + |A,9] <e+ Ml, >a 


Therefore, ||A||< M. 
Finally we have 


Lemma 6. If for all fe D, where D is a linear manifold dense in 
Fe, | Af <elfl, |A|<-M, and c is independent of f, then ||A|< MinfM, ch. 


Proof. It clearly suffices to assume that M > 0. By hypothesis corres- 
ponding to any gEé2, {g,} exists with g,e 2 such that In >g. Therefore for 
for any e> 0, N exists such that |g —g,|<e/(M-+c), n> N. For such n we 
have the relations 


|-Ag | Ss A (da In) = | Aga <M- i re [gx] < 
MSI 


<5 38 + els — aal+ olo] <e+ ell 
Mag a! Ri 


which shows that |A|<c. This result together with the assumption that 
|A|| <M proves the lemma. 

On combining Lemmas 4, 5 and 6 we obtain Theorem 1. To see this we 
notice that {#,U,} satisfies the requirements of Lemma 5 if t,> + co as 
n—> oo, since 1), ||E,U.|<|Z;|:|U,|<%; independently of n, and 2) EU, 
satisfies the second requirement on account of the Corollary to Lemma 4. 
Therefore, {H,U,,} converges on all # to an operator A such that | Al) <1. 
But on the dense domain of Lemma 4, A satisfies the requirements of Lemma 6 
with c=0. Therefore E, U, +0 if t, > +t co asn+ co. This convergence 
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to zero is however independent of the particular sequence {t,}. Therefore 
Hy U,=-0 on all of df as t>+oo. (qe.d.). 


2. — Uniqueness of the free Hamiltonians of a scattering system. 


The wave operator of a single-channel scattering system has been studied 
by many authors (1°) and has been recently defined most succinctly and pre- 
cisely by JAUCH (*). In a more recent paper (?) JAucH has extended this for- 
mulation to the case of a multi-channel scattering system. Since in the present 
paper we are interested in such multi-channel scattering systems, we shall 
outline the definition of such systems and summarize the properties of the 
associated wave operators needed for the proof of Theorem 2. 

A scattering system is characterized by a countable set of non-interacting 
fragments at t=— co and by another sueh set at t= + co. During finite 
times the incoming fragments are in a finite region of space where they interact 
and yield a long time later (t ++ co) what is recognized experimentally as 
the second set. Each fragment at t= + co is assumed to consist of a bound 
set of fundamental particles, and each fragment is characterized by a mo- 
mentum operator P,, which is the sum of the momentum operators of the 
particles constituting it. The free or asymptotic motion of all the fragments 
is described by the self-adjoint Hamiltonian 


1 il 
(8) H,= PIA E 


= | 
2m, 2M, 


Li 
oy all 


where x is a label specifying the channel energy, and the set {m,} consists of 

the masses of the fragments appearing in this channel. The self-adjoint Hamil- 

tonian of the entire system is H and the interaction is defined by H— H,. 
Let 


(9) UP = exp [— Ht] } V, = exp [— iHt], 

where {U®} and {V,} are unitary groups which clearly satisfy the relation 

(10) IS Ue eae 0% forall tae pe ti, to. 

It is further assumed or it can be proved from the assumptions (2) that if € 
(19) C. Motter: Kgl. Dan. Videnskab. Selskab, Mat.-fys. Medd., 23, 1 (1945); 


B. Lippmann and J. SCHWINGER: Phys. Rev., 79, 469 (1950); M. GELL-MANN and 
M. L. GOLDBERGER: Phys. Rev., 91, 398 (1953). 
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is the underlying abstract Hilbert space, then (A) a set of bounded wave 
operators {Q®} eists such that lim V*U® exists on a closed domain 2D, 


to 
and defines Q\ as a partial isometry fan isometry on D,; 0 on the orthogonal 
complement, DI, of D,); (B) the ranges of QY are FR (and of course the. 
R® contain elements different from and in addition to the trivial null ele- 
ment) and the closed linear manifolds R, formed from the sets {R®} and 
{R®}, respectively, are such that R,= R_- N, where N is the sub-space of 
ge associated with the continuous spectrum of H; (C) the number of channels. 
is at most denumerably infinite. With these remarks in mind we may 


now state 


THROREM 2. Let {U®} and {U®}, «=1,2,..., be two sets of unitary 
eroups such that 


CLI) UY? = exp|—iH,t]; U® — exp[— Hi], 
where 
(12) Da -> dl Se Memes be SC oe 

VIZIO! 


and {a®}, {af} are sets of real numbers. Also let {V,} be a unitary group: 
such that for each « lim V*UY and lim VÎU® exist on closed domains. 


t+ o t>Fo 


D, and D,, respectively, and define the partial isometries QF and Q®, 
with ranges Rand R™, respectively. Finally assume that R, —R,=R_=R_=N. 
Then the sets {H } and VDR } are identical, except possibly for the order of their 
elements. 

What this theorem states is that for any multi-channel scattering system 
as defined in (A4)-(C) there is only one set of non-relativistic free Hamiltonians 
available to describe the free or asymptotic motion of the fragments. Ex- 
pressed otherwise, the theorem states that if the total Hamiltonian can be 
split into two parts of which one part, H,, is suitable for the definition of a 
scattering process, this splitting can be done in only one way. One should 
add to this formulation the provision that the H, have the form of (12). 
However, it will be clear that the theorem can also be proved when relativistic 
expressions replace those in (12). Still it remains a moot point whether 
Theorem 2 is true if the form of H, is unspecified. Relative to this more 
general problem Theorem 2 may be regarded as a lemma. 

The proof of the theorem consists in showing (in Theorem 3) that if H, 
and H, are of the form (12) such that H,#H,, and if in conjunction with a. 
given unitary group {V,} they define wave operators QY and Di as partial 
isometries, then the ranges RQ and R” of these wave a are ortho- 
gonal to each other. por to the sae {H,} and {H,} of Theorem 2, this. 
result implies that if one «’ exists such that H, 4H, for all «, RY is ortho- 


H 
~ 
wD 
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gonal to all R®, that is, to V, and consequently also to R, and finally to 
itself, since by assumption R®°CR,. Therefore R® contains no non-null ele- 
ment, contrary to the assumptions. Consequently {H,} c {H,}. Similarly, one 
must conclude that {H,} c{H,}. Therefore, {H,} = {H,}. (q.e.d.). The proof 
of Theorem 2 will therefore be completed with the proof of 


THEOREM 3. Let {U,} and {U,} be two unitary groups defined for 
—oco<i< co such that 


(13) U, = exp [—1H,t] ; U, = exp[— i,t], 

where 

(14) 13 ee eae Hie Oe ds are 208 
i=1 î=1 


Also, let {V,} be a unitary group (— co<#< co) such that lim Vi U,f and 
lim V*U,f when #t-+ 7 co exist for arbitrary elements f and f in the closed 
domains 2 and 9, respectively. Then if H,4H,, R_LR_ and RESA 
where R_, R_, R,, R, are the ranges of the partial isometries OOO RIO 
respectively, which are defined by the three groups {U}, {U.} and {V.}. 

The proof of this theorem will be carried out with the aid of three lemmas, 
the first of which is 


Lemma 7. Let {U;} and ST} be two unitary groups defined for — co < 
<t< co such that 


(15) U,= exp[— iHoi]; Ul exp Het, 

where 

(16) nh, SOM Pic i= > a. Pr; ai, a; > 0. 
i=1 i=1 


Mieneit J,-H, 
(17) U*U,—0 as bas Logs 


Proof. Just as in the proof of Theorem 1 introduce as a realization of 
the underlying Hilbert space # the set F(k,) of L2-functions defined on the 
entire 3n-dimensional Euclidean space Egn. 3(k;) has properties similar to 
those espressed in (1) and (2), 4.65 


(18) w(P?, P? DICE] P?)g(k,, k., UL) kn) = (ki, ki, wales ki) p(k, ko, erie k,) bl 


2) n 


96 I. I. ZINNES 


where (k,, k,,.... k,)=(k,) is an arbitrary element of f(k,) such that 
|o(P2)p| < 00. 

Consider first two special elements f and g in # which are represented 
in 3(k;) by the same step-function: 


(od oo 


(19) {k)=g(k)= 4 hii 
[| 0, otherwise, 


where A is a rectangular region in H,, whose sides are parallel to the axes of 
#E;, and whose measure is finite. Then 


(20) 100) = (f, 070g) =| exp(—it > (a —@)k]dr, 
È =1 


where dt is the 3n-dimensional volume element in £,,. J(t) can be written 
as a 3-fold iterated integral, each of whose iterates is a one-dimensional integral 
having the form 


o 


(21) J (t) =[exp [— ittqk?| dk , 


O71 
3 


where q=a,—a, for some s and, if ki= > ki 


sa) 


DI 


kh? =k, for some s ‘andile 
oa 


Also, (01, 62) is the interval over which % can vary in A. Since according to 
the assumptions H, ~H,, at least one value of r exists for which a,—a, 40. 
For such values of r 


d0) 


n at 
I) = 57 i ee aa, 


2 
107 


The range of integration may be extended from — co to + co if the integrand 
is replaced by the function (4) defined in the following way: 


dea 
wd) = VA 
0, 


otherwise. 


2 2 
» Goi <A<qo? , 


Since y(4) is absolutely integrable over the whole line, J(t) > 0 ast > + oc 
according to the Riemann-Lebesgue lemma (11). For other values of r, i.e.. 


(!!) E. C. TircHMaRsH: The Theory of Fourier Integrals (Oxford, 1937). 
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those for which a,—a,=0, q=0 and 
= (f, 0; Ug) —>0 as t > + 0. 

Although this result has been obtained for the very special step-function (19), 
it is also true for arbitrary square-integrable step-functions, since the scalar 
product is bilinear in its two arguments. Since the set of square-integrable 
step-functions is dense in # (#2), the lemma has been proved for elements in 
a domain 7, dense in &#, le., 


J(t)|<0.— 01. Therefore Le: 


(22) REMO ESA for {,,9€7. 


To complete the proof of this lemma we must show that (22) is true for 
elements in the closure 7 of 7. This is most easily shown with the aid of the 
following variant of the Banach-Steinhaus Theorem. 

Lemma 8. Let {T,} be a sequence of bounded operators such that 

) |Z.| <M for all n 
2) (7,p,q) converges as n + co for each p and q in a domain 7 which 
is dense in DC. 


Then 7, —T, where 7 is bounded and |7|<M. It will be noticed that 
this lemma differs from Lemma 5 only in the types of convergences consi. 
dered. Neither lemma implies the other, however. Since the theorem is 
trivial if 17 = 0, we shall assume that M > 0. 


Proof. Let f and g be any two non-null elements in Sf. Then, since 7 
is dense in 20, sequences {f,} and {gn} exist in 7 such that fr >f and gn > 9g 
as n —=oc0. Then 
(23) MERE 9) sa (Et, 9) | sai (ET 9) MIT) 9) ise 
(Li Ladies g gs) | owas EEN ORG AN ase 
<= Lol lf —fell Val +1 Lo— Lol tel Lg — gol +1 (LZ 9) 1< 
<2M-|f—fell-|g| + 2M Fell 19 —gel + | (Pn — Tmdfor 99. 


Since f.>f/+0 and g,>g#A0 as 7, § > 00, R and S exist such that f, 40 and 


(24) ff] <e/6 2: lg —9s| <e/6M [fal - 


According to assumption 2, N exists such that for fixed R and S and for all 
n, m>N 


(25) | (Ln — Tm)far 98) |< 88 


(12) E. J. MoSHANE: Integration (Princeton, 1944), p. 230. 


7 - Supplemento al Nuovo Cimento. 
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Setting r= R and s=S in (23) as determined in (24), we obtain 


(26) (Luh, 9) —(Tmhs 9) |< 2M aoe pay lll + 


€ 


Vel Gar. [fal 239 


+2M- 


ee te M SING 


If either or both the relations: f= 0, g= 0 hold, the left-hand side of (26) 
automatically vanishes. Therefore (26) is true for all elements in # and 
asserts that 7 exists such that 7, —T as n+ oo. 

The boundedness of 7 follows from the result just obtained and the re- 
lation 


M?|f 


>| (Lat, Lf) |—> (LF, T)=|Tf 
as n, n'> oo, a relation which is correct for arbitrary elements of dé. 


Corollary. If {T,} has the properties required in Lemma 8, and if in ad- 
dition (7,p,9) >0 asn+ co for all p and gq in J, then (T, f, g) > 0 as n + co 
for all j and g in J= 93€. 


Proof. According to Lemma 8, 7,,—T as noo and T is bounded. 
Therefore, for arbitrary f and g, (7,f,9) >(Tf, g) as n > co. If f and ge7, 
(Tf,g)=0 by assumption. If f and g e 7, they are each limits of sequences 
of elements in 7. And for elements of these sequences (7f,, gm) =0, where 
In >f and gg >9g a8 n-+ co. Since 7 is bounded, (Tf,, 9) is continuous in 
f, and Ym, and therefore [continuity of the scalar product] vanishes also when, 
NM, M—>co. (g.6.4.). 


The corollary just proved completes the proof of Lemma 7. This can be 
seen in the following way. In (22) it was shown that (f, 07U,g) >0 as 
t>x+tco for j and ge 7, where 7 is dense in 30. This result implies that for 
elements of J, (f, Uj,U,g) >0 as n+co if t,-> + c0 as n1+ co. In addi- 
tion, since U, and U; are unitary, ||0,°U,|<1. Therefore, since {TT U} sa- 
tisfies the requirements of the corollary, we conclude that UU, = 0 as nm + co. 
Finally, since this result holds for arbitrary sequences {t,} such that t, > +0 
as n>co, UU, ->0 as t > +o; and this is just the conclusion of Lemma 7. 


We now prove Lemma 9, the final lemma required for the proof of Theo- 
rem 3. 


Lemma 9. Let {U} and {U;} satisfy the requirements of Lemma 7. In 
addition let {V,} be another unitary family defined for — co<t< co with 
the property that for the two non-null elements f and f at least one of the 


cel 
D 
to) 
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following two pairs of limits exists: 


(27), Vil VEE as t>+ co, 
(27) RI Vetere Hypa os 
Then 
(28). Ei 

2 (28) 0 


Proof. Since the proof of (28)_ is similar to that of (28),, we shall prove 
only (28),. According to (27), and because the scalar product is continuous 
in its two arguments, 


(ff) = lim (MEUS, VU.) 


Since the left-hand side is independent of the exact manner in which ¢, #' 0, 
we may take t=?’. Then 


(f-, f) = lim (V5U.j, V70,7) = lim (T7U,j, f) = 0- 


The last equality follows from Lemma 7. (q.¢.d.). 

Theorem 3 is now effectively a trivial corollary of Lemma 9, since according 
to the lemma each element of D and of Min Theorem 3 is mapped by the iso- 
metries into two corresponding elements, which are mutually orthogonal. 

As stated earlier (see paragraph preceding the statement of Theorem 3), 
proof of Theorem 3 is equivalent with the proof of Theorem 2. With the proof 
of Theorem 3 we have therefore completed the proof of the main theorem. 


OK 


It is a great pleasure to acknowledge my indebtedness to Professor J. M. 
Jaucu for introducing me to the problems of this paper, for emphasizing the 
practical power available in the formal theory of HILBERT space, and for many 
stimulating discussions. I also wish to thank him most warmly for the generous 
hospitality he made available to me at the State University of Iowa. 


579 


SUPPLEMENTO AL VOLUME XII, SERIE X N. 1, 1959 
20 i ° 
DEL NUOVO CIMENTO 20 Trimestre 


Gli organismi biologici come sistemi aperti stazionari 
nel modello teorico di L. von Bertalanffy (*). 


R. MorcHIO 


Istituto di Fisica Teorica dell Università - Genova 


(ricevuto il 4 Maggio 1959) 


Sommario. — 1. Introduzione. — 2. Le proposte del von Bertalanffy. — 
3. Aspetti cinetici e termodinamici dei « sistemi aperti stazionari »). — 
4. Problemi e conseguenze del modello proposto. — 5. Formalizzabilità 


e accordo con l’esperienza del modello. — 6. Operatività del modello. — 
7. Conclusione. 


1. — Introduzione. 


Il lavoro del biologo che affronta i tradizionali problemi della biologia 
secondo le tradizionali modalità, appare, a prima vista, metodologicamente 
diverso non soltanto dal tipo di lavoro del fisico, ma anche del biofisico 0 del 
biochimico, che, nei capitoli più recenti delle scienze biologiche, sempre più 
tendono ad un ampio uso della formalizzazione matematica. 

Se però l'oggetto ultimo della ricerca dei biologi, dei biochimici e dei bio- 
fisici è fondamentalmente lo stesso, perchè questa unità di ricerca non sia 
limitata a contenuti più intuitivi che teorici, è necessario che proprio sul piano 
metodologico si arrivi ad un sostanziale chiarimento nel senso del superamento 
di ogni residuo dogmatismo sulla non applicabilità di metodi rigorosi e quan- 
titativi (*) alla ricerca biologica a tutti i suoi livelli, che è poi quanto riaffer- 
mare in altri termini la omogeneità sostanziale della ricerca sulla natura. 

(*) Lavoro fatto con un contributo del C.N.R. 

(1) Nota giustamente il WooDGER (Biology and Language, Cambridge, (1952)) che 
veramente il problema non è tanto quello di rendere la biologia più quantitativa, 


quanto di trasformare i metodi « intensionali » in « estensionali ». Essendo però la ma- 


tematica essenzialmente « estensionale » in prima approssimazione le due affermazioni 
non differiscono di molto. 
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L'affermazione dell’applicabilità dei metodi quantitativi alla ricerca bio- 
logica implica la effettiva applicazione concreta di tali metodi ai tradizionali 
problemi biologici (essa è già avvenuta in qualche singolo caso e per il resto 
sarà frutto della ricerca); ma può essere sostenuta, in linea teorica, anche 
mediante la dimostrata proponibilità di modelli teorici dell'organismo, che, 
suscettibili di formulazione matematica, siano in accordo con l’esperienza. 
In questo senso vale la pena esaminare se esista qualche modello rispondente 
alle condizioni accennate: che non sia cioè di solo interesse filosofico (*?) ma 
che possa costituire piuttosto una ipotesi di lavoro. 


2. — Le proposte del von Bertalanffy. 


LUDWING VON BERTALANFFY, il noto biologo viennese che vive attualmente 
negli Stati Uniti, ha proposto già da molti anni di considerare l'organismo 
biologico come un « sistema », complesso fin che si vuole, rispondente a sue 
particolari leggi di comportamento, ma aggredibile e studiabile con metodi 
di tipo fisicalistico. Naturalmente la parola « sistema » non bastava. Bisognava 
rispondere a problemi e difficoltà ben note, lasciando spazio per le soluzioni 
originali che i problemi originali dell'organismo biologico richiedono. Per questo 
il von BERTALANFFY fin dal 1932 (*) avanzò la concezione dell'organismo come 
un «sistema aperto », un sistema cioò in cui lo scambio di materia con l'esterno 
modificava completamente lo schema del « sistema chiuso » descritto dalla 
fisica ed in particolare dalla termodinamica ordinaria. 

Le difficoltà in ogni modo erano notevoli. Da un lato la biologia ha spesso 
usato espressioni come « sistema », « organismo », «interazione », « il tutto è più 
che la somma delle sue parti », « equilibri biologici », e simili, ma sempre, come 
nota lo stesso VON BERTALANFFY (4), in modo vago e talora mistico; dall'altro 
la stessa descrizione dei « sistemi chiusi» della fisica era assolutamente inap- 
plicabile agli organismi viventi e non poteva dare molti aiuti. Scrive il PRI- 
GOGINE, a proposito della termodinamica: « Una seria limitazione della ter- 
modinamica classica come strumento generale per una descrizione macro- 
scopica dei processi fisico-chimici sta nel fatto che il suo metodo è basato sui 
concetti di « processi reversibili » e sugli stati di vero equilibrio » (°). Non a caso 


(2) Anche il lavoro del MAINX (Foundations of Biology, in: Intern. Encycl. of Unified 
Science, vol. 1, parte 28, (Chicago, 1955), p. 567), pur interessante per molti aspetti, 
non soddisfa questo punto di vista. 

(3) L. von BERTALANFFY: Theoretische Biologie, Band II, I ed. (Berlin, 1932-42) e 
Naturwiss., 28, 521 (1940). 

(4) L. von BERTALANFFY: Brit. Journ. Phil. Sci., 2, 134 (1950). 

(5) I. PRIGOGINE: Introduction to Thermodynamics of Irreversible Processes (Spring- 


field, 1955), p. V. 
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proprio su questo fatto nacque Vequivoco di considerare la termodinamica 
ordinaria, a causa della sua inapplicabilità agli organismi viventi, quale prova 
della differenza sostanziale e incolmabile tra vita e non vita. Il Lecomte du 
Notiy, in particolare, diffuse tra i biologi questo argomento che acquistò in 
molti ambienti un notevole prestigio, anche se immeritato. 

Per tutti questi motivi il von BERTALANFFY da circa vent’anni ha sotto- 
lineato la necessità di generalizzare alcuni dei concetti della fisica e della chimica, 
perchè potessero essere estesi anche ai « sistemi aperti ». Lo stesso VON BER- 
TALANFFY (6) e il BURTON (7) cominciarono a sviluppare qualche principio 
della cinetica dei « sistemi aperti »; mentre i processi irreversibili e i « sistemi 
aperti» stessi sono stati studiati termodinamicamente in particolare dalla 
scuola belga: DONDER, ONSAGER, PRIGOGINE e altri (8). 

Diveniva così possibile al von BERTALANFFY raccogliere tutti questi elementi 
e delineare un primo schema di un modello d’interpretazione dell’organismo 
vivente come «sistema aperto ». 

Per la verità questa visione dell'organismo fu ampiamente generalizzata 
dal vox BERTALANFFY. La sua conclusione fu che ciascuna scienza studiando 
dei « sistemi)» e presentando relazioni matematiche, leggi, principi formal. 
mente equivalenti, sarebbe possibile formulare alcune relazioni generalissime, 
le quali, riferendosi ad un sistema in generale, potrebbero valere per tutte 
le scienze, rendendo finalmente possibile, su un piano formale, un’effettiva 
unità delle scienze stesse. Tali concetti costitutiscono il fondamento di quella 
che il von BERTALANFFY chiamò « General System Theory », strumento, nel 
pensiero dell’autore, di controllo e di spinta a trasferire principi da un campo 
all’altro. La estrema generalità della « General System Theory » (?) la pone, 
in questa sua forma, al di fuori dei limiti e dell’interesse di questo lavoro. 

Qui non si vuole, del resto nè esporre dettagliatamente il modello dell’or- 
ganismo come «sistema aperto », nè analizzarne la validità in modo esauriente. 
Si vuole solo esaminare se esso è proponibile. Per questo è sufficiente un breve 
« survey » della teoria e delle eventuali conferme. Sarà l’affermazione di pro- 
ponibilità del modello che potrà aprire, in un secondo tempo, un discorso critico 
più impegnato. 


) L. von BERTALANFFY: Naturwiss., 28, 521 (1940), già citato. 
) A. C. Burton: Journ, Cell. Comp. Physiol., 14, 327 (1939). 

($) Vedi in particolare: I. PRIGOGINE: Etude thermodynamique des phenomènes irré- 
versibles (Parigi, 1947) e Introduction to Thermodynamics of Irreversible Processes (già 
citato); S. R. DE Groor: Thermodynamics of Irreversible Processes (New York, 1951) 
I. PRIGOGINE e J. M. WIAME: Experientia, 2, 450 (1946). 

(*) Il von BERTALANFFY, citando il filosofo N. HARTMANN, parla della G.S.T. addi- 
rittura come di uno strumento per rimpiazzare la teoria delle categorie con un sistema 
esatto di leggi logico-matematiche. 
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Ma che significa chiedersi se questo modello è proponibile? Nello spirito 
di quanto è stato detto si possono delineare i seguenti criteri: a) che il « modello » 
sia formalizzabile, sia suscettibile cioè di una formulazione matematica, d) che 
non sia, ovviamente, in disaccordo con l’esperienza, e €) che sia operativo, 
che permetta cioè di trarre dalla sua formulazione generale conseguenze pas- 
sibili di verifica sperimentale. 

I punti essenziali del modello, che verranno qui brevemente proposti, 
saranno valutati appunto in base a questi criteri. 


3. — Aspetti cinetici e termodinamici dei « sistemi aperti stazionari ». 


Il vox BERTALANFFY intende in generale per «sistema»: « Un complesso 
di elementi p,, P», -.- Pn interagenti. Interazione significa che gli elementi 
stanno in una certa relazione R, cosicchè il loro comportamento in R è diverso 
dal loro comportamento in un’altra relazione R'. Se il comportamento in Io, 
ed in R' non è differente, non vi è interazione e gli elementi si comportano 
indipendentemente rispetto alle relazioni & e Re ( 20s 

Un tale sistema pud biologicamente riferirsi, a livelli diversi, ad una popo- 
lazione (gli individui ne costituiscono gli elementi), ad un organismo solo (gli 
organi, o le cellule ne costituiranno gli elementi), ad una cellula, ece. Tali 
sistemi possono essere descritti da. sistemi di equazioni ben note. Indicando 
la misura di qualche aspetto quantitativo degli elementi p di cui sopra con 
Q,, 92, Qa, ---1 Qn, la sua variazione può essere definita appunto da un sistema 
di equazioni differenziali simultanee che per un numero finito di elementi, 
nel caso più semplice, avrà la forma: 


dQ, 

de CORO), 

dQ» 

a aie Ore ( 23 teh ’ 
È, Ap CITAZIONI 

dQn — i 

dg sh fa(01, Gay oo On). 


| Tale sistema esprime, secondo VON BERTALANFFY, un principio generale 
di cinetica valido per qualunque sistema; ovvio che, a livelli diversi, sarà diversa 
l'applicazione, come diverso sarà il significato dei parametri e delle variabili. 
Chiaro è il significato di queste equazioni. È pertanto sufficiente qui ricor- 


(1) L. von BERTALANFFY: Brit. Journ. Phil. Sci., 2, 134 (1950), già citato. 
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darne le conseguenze solo per alcuni casi particolarmente semplici ma non 
privi di importanza per la biologia, come si vedrà successivamente. 
Se il sistema studiato è composto di elementi di una sola specie. 1°(1) diviene 
semplicemente: 
nl f(Q) 
in 
che, supposta la funzione analitica, sviluppata in serie di Taylor nell'intorno 
di Q=0, diviene: 
do 
(2) ae DAMIANI LES 
dt 
Se della (2) si considera soltanto il primo termine del membro di destra, 
se cioè l'accrescimento (0 il decremento) numerico di Q dipende esclusivamente 
in qualche maniera dallo stesso valore che assume @ in ogni singolo momento, 
si ha il caso più semplice: 
dQ 


ap ae 


cioè integrando: 


(= o exp [at] ’ 


indicante un incremento o un decremento (a, positivo o negativo) esponenziale. 
Considerando invece i primi due termini del membro di destra della (2), 


nel caso che l'accrescimento numerico (o il decremento) di Q sia funzione 
non soltanto del valore dello stesso Q, ma anche di un altro fattore, si ha: 


dQ 


Fes = 4,9 + dn? , 


che, integrata, da la soluzione: 


rere a, C exp[a,t] 
1—a,,C exp[a,t]’ 


descrivente un accrescimento indicato da una curva logistica. 

Nel caso di un sistema consistente di due specie di elementi, il sistema si 
svolgera secondo linee diverse a seconda dei casi. Sara rappresentato dalle 
due equazioni: 


dq 
di 


dQ» ) 
o fa(Q1 302), 


>= f(0,,,03), 


x 
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che, svolte in serie di Taylor e risolte, in base al valore assunto caso per caso 
dai parametri, daranno tre possibilità: 1) descrizione di un sistema che tende 
verso uno stato stabilmente stazionario (si vedrà poi quale significato esatto 
assuma la parola « stazionario »), 2) descrizione, mediante termini periodici, 
di un sistema che oscilli attorno a valori stazionari, 3) descrizione di un sistema 
non in stato stazionario, nè oscillatorio. 

Gli aspetti termodinamici dei «sistemi aperti» (ne abbiamo visto quelli 
cinetici) devono essere esaminati tenendo conto, come s'è detto della recente 
termodinamica dei processi irreversibili. Di essa basta ricordare qui alcuni 
punti, rimandando ai testi indicati nella nota (*) per una maggior docu- 
mentazione. 

La termodinamica ordinaria ha a che fare con sistemi chiusi, sistemi per 
i quali non si ha scambio di materia con l’ambiente esterno, ma solo di energia. 
Per essi vale la relazione di Carnot-Clausius: 


7 È , d 
dS =d.8S+d4,8 = vl + d,S 


in cui S indica l'entropia, gli indici e ed è si riferiscono rispettivamente all’en- 
tropia scambiata dal sistema con l'esterno (e) e a quella prodotta all’interno (2) 
del sistema dai processi irreversibili, dQ è il calore fornito al sistema dall’am- 
biente e 7 la temperatura (assoluta). d;S può essere soltanto zero o positivo. 
La termodinamica ordinaria considera il caso in cui d,S=0 che dà la seconda 
legge per i processi reversibili (aS =dQ/T). 

Se abbiamo d;S> 0, cioè se dS >dQ/T, che si riferisce evidentemente 
ad un processo irreversibile, la termodinamica ordinaria descrive quest’ultimo 
come un cambiamento infinitamente lento, da uno stato di equilibrio ad un 
altro. Questo schema, per essere utilizzabile biologicamente, deve evidente- 
mente esser reso più complesso. | 

Già qualche sistema chiuso può ampliarlo: il secondo principio richiede 
che in un sistema chiuso la produzione di entropia risultante da eventuali rea- 
zioni simultanee sia positiva. Può però accadere che, indicando con A l'affinità 
chimica e con v la velocità di reazione, un sistema subisca due reazioni simul- 
tanee tali che A,v,< 0 e A30,>0, purchè 4,0; +A2% > 0. Tali reazioni sono 
dette reazioni accoppiate, ed è così possibile che una delle reazioni progredisca 
in direzione contraria a quella prescritta dalla sua affinità. Così nella termo- 
diffusione, ad esempio, si ha una diffusione di materia che avviene contro il 
suo gradiente di concentrazione, producendo entropia negativa, compensata 
peraltro dalla produzione di entropia positiva dovuta al flusso di SUO) Tale 
descrizione ha già un interesse biologico perchè negli organismi viventi sono 
estremamente importanti le reazioni accoppiate. 

Ma ancor più interessante sa ‘à questa descrizione Sse, prendendo in consi- 
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derazione i sistemi aperti, la relazione di Carnot-Clausius, in d,S, verra comple- 
tata con un termine ulteriore rappresentante la fornitura di entropia da parte 
dell’ambiente attraverso il trasporto di materia. Ora se l’entropia assunta 
dal sistema aperto in relazione ad una assunzione di materia è negativa e 
pareggia la produzione all’interno del sistema, si ha una situazione che non 
è di equilibrio nel senso della termodinamica dei processi reversibili perchè 
il sistema subisce delle trasformazioni, ma queste trasformazioni sono tali 
per cui il loro risultato complessivo è di mantenere il sistema in una situazione 
costante. La termodinamica dei processi irreversibili chiama questo stato 
stazionario. 

Un esempio di stato stazionario può esser dato da un sistema che riceve 
in continuazione un componente M dell’ambiente esterno e lo trasforma con- 
tinuamente, attraverso un certo numero di composti intermedi, in un pro- 
dotto finale # che è restituito in continuazione all'ambiente esterno, mentre 
le concentrazioni dei componenti intermedi non variano col tempo. La con- 
dizione perchè si abbia appunto lo stato stazionario stà nel fatto che la pro- 
duzione di entropia abbia il suo valore minimo compatibile con le condizioni 
stabilite caso per caso dalla concreta situazione e che, nell'esempio testè pre- 
sentato, potrebbero essere le concentrazioni dei prodotti iniziali e finali M e F 
nell'ambiente esterno cui il sistema deve adeguarsi. 

I problemi relativi agli stati stazionari, per il loro interesse, sono stati 
discussi non solo da fisici, ma anche da biologi (!!). La termodinamica dei pro- 
cessi irreversibili ha portato infatti a conseguenze di non poco interesse bio- 
logico, dimostrando che, operando sempre i processi irreversibili interni in 
modo da abbassare il valore della produzione di entropia per unità di tempo, 
quando un sistema è in uno stato di minima produzione di entropia, allora 
non può lasciare questo stato mediante un cambiamento irreversibile spon- 
taneo. Se per qualche fluttuazione devia leggermente da questo stato, avranno 
luogo cambiamenti interni che riporteranno indietro il sistema al suo stato 
iniziale, il quale può considerarsi perciò uno stato stabile come pure stabili 
possono considerarsi le trasformazioni che appunto in tale stato intervengono. 

_In generale si può dire dunque, più esattamente, che in un sistema carat- 
terizzato da n forze indipendenti X,, X,,..., X,, lo stato stazionario consiste 
nella presenza di un certo numero di costrizioni, fissate da un certo numero 
delle n forze X,, X,,..., X, (k è uno dei numeri 0, 1, 2, ..., n) tenute a valori 
costanti (stato stazionario di ordine k). In un tale sistema, con una produ- 
zione minima di entropia o, si vanificano i flussi J; con gli indici j=%+1, 
Bae vie Ny COS ag I ey, COVE, Clascuned 6 1 asso legato 
alla relativa forza X;. Questo teorema, che fu formulato per primo da Pri- 


(!!) Trai biologi sono da ricordare: A. J. LorKa: Proc. Nat. Acad. Sci., 8, 146 (1922); 
Bull. Math. Biophys., 10, 103 (1948) e J. Z. HEARON: Bull. Math. Biophys., 12, 57 (1950). 
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gogine per il caso di un singolo parametro fissato, è una generalizzazione del 
principio di moderazione di Le Chatelier. 

È ovvio che in uno stato stazionario una produzione di entropia dev'essere 
compensata da un flusso negativo di entropia in modo che la variazione totale 
di entropia nel tempo sia zero: 


as GUIS ds 


—— = — 0 
dt ieee 3 
Ma: 
d,S = 0 
di Re; 
per cui: 
das 
0. 
di 


L’entropia della materia che entra in un sistema in stato stazionario deve 
perciò essere più piccola dell’entropia della materia che il sistema restituisce 
al mondo esterno. Si ha perciò una degradazione della materia che entra nel 
sistema, degradazione che mantiene lo stato stazionario, contribuendo all’au- 
mento dell’entropia esterna. Durante l'evoluzione verso lo stato stazionario 
la produzione di entropia e spesso la stessa entropia decrescono prendendo 
jl valore più basso compatibile con le costrizioni esterne quando lo stato sta- 
zionario è raggiunto. 

Questo tipo di descrizione è particolarmente efficace per il comportamento 
globale degli organismi viventi a differenza della descrizione dei sistemi della 
termodinamica ordinaria (!*). L'evoluzione degli organismi fino allo stato 
adulto può vedersi come l'evoluzione del sistema verso lo stato stazionario, 
avente luogo sotto un certo numero di costrizioni determinate dal mondo 
esterno e questo comporta appunto una diminuzione nella produzione di entropia. 
Il fatto che l’organizzazione di un simile sistema durante l'accrescimento 
aumenti, spiega la diminuzione di entropia, mentre la sua stabilità è chiarita 
dalla stabilità dei processi stazionari e la sua capacità a riequilibrarsi, dopo 
modificazioni non eccessive subite, è spiegata in base al principio di Le Chatelier 
generalizzato. Il continuo dinamismo dei processi metabolici in cui quelli 
anabolici pareggiano quelli catabolici ed il continuo rinnovarsi della materia 
componente gli organismi, coesistenti con la permanenza della forma e delle 
strutture così autoregolantisi, trovano dunque una chiarificazione nella dina- 
mica degli stati stazionari. 


(12) Si veda in particolare lo scritto di PRIGOGINE e WIAME in Experientia (1946) 
già citato. Anche: L. von BERTALANFFY: Science, 3, 23 (1950). 
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4. — Problemi e conseguenze del modello proposto. 


Il modello proposto consiste perciò nel considerare un organismo vivente 
come un sistema aperto in uno stato stazionario e nel descrivere tale sistema 
cineticamente con le equazioni differenziali indicate all’inizio della Sezione 3 
e termodinamicamente con le equazioni della termodinamica dei processi irre- 
versibili. 

Tale descrizione è sostanzialmente omogenea non sussistendo contraddizioni 
tra i due tipi di descrizione. Non è qui il caso, naturalmente, di esaminare se 
tale descrizione implichi o meno l'assioma di causalità come qualcuno ha fatto (7*), 
ma è utile sottolineare come questa sia una descrizione globale del comporta- 
mento di un organismo; o anche di sotto-sistemi di un organismo, ma sempre 
nei loro comportamenti globali. Questo modo di vedere un organismo, a nostro 
avviso, apre proprio il problema di come questo modello possa eventualmente 
accordarsi con la considerazione di singoli eventi microscopici, notevolmente 
importanti in un vivente. 

Si può anche porre il quesito di quanto grande sia il contributo del von 
BERTALANFFY nell’elaborazione di questo modello. La termodinamica dei pro- 
cessi irreversibili è stata sviluppata da un gruppo di fisici. Le equazioni diffe- 
renziali cinetiche non si devono completamente a lui; del resto la legge del- 
l’azione di massa poteva già dare una indicazione. Il von BERTALANFFY però, 
come è stato detto, già da tempo aveva intuito che un organismo poteva con- 
siderarsi come sistema aperto, prima che la termodinamica dei processi irre- 
versibili fornisse precise indicazioni in proposito. Inoltre, e non è poco, il Nostro 
raccolse tutti questi elementi e contribuì notevolmente a chiarirne l'interesse 
biologico ed a farli conoscere ai biologi. Infine utilizzò queste idee proprio 
all’interno delle classiche concezioni della biologia tradizionale. Su quest’ultimo 
punto, sia pure brevemente, vale la pena spender qualche parola (14). 

L’unitarietà del comportamento di un organismo vivente assume, ad esempio, 
per VOn BERTALANFFY un significato ben preciso cineticamente e termodina- 
micamente. Certe espressioni generiche come: «una modificazione in una parte 
influenza il tutto », ecc., si traducono da un lato nell’affermazione termodina- 
mica secondo cui una qualsiasi modificazione nel sistema determina comples- 
sivamente nel sistema stesso tali modificazioni da ristabilire l'equilibrio, dal- 
l’altro lato nella constatazione che dall’equazione, sviluppata in serie di Taylor: 


dq ; 


di = Wi1Qi + 49004... + indu + di1107 Soto) 


(13) W. YourGRAv: Scientia, 11, 307 (1952). 

(14) Indicazioni maggiori possono vedersi in: L. von BERTLANEFY: Brit. Journ. 
Phil. Sci., 2, 134 (1950), già citato; Human Biology, 23, 302 (1951); Journ. Personality, 
20, 24 (1951); Main Currents in Modern Thought, 4, 75 (1955). 
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risulta che ogni mutamento di una quantità Y,, è funzione polinomiale di tutti 
gli elementi @,, Q...., Qn. A sua volta un mutamento di uno dei Q; causa 
una variazione in tutti gli altri Y; e quindi nel sistema nel suo complesso. Di 
particolare interesse è l'applicazione che il von BERTALANFFY fa di questi 
concetti a quel tipico comportamento biologico che il Nostro ha definito 
« equifinalita ». In un sistema fisico ordinario lo stato finale è determinato 
dalle condizioni iniziali; in un organismo vivente invece lo stato finale può 
essere raggiunto partendo da condizioni iniziali diverse e attraverso vie diverse. 
Il famoso esperimento su un germe di riccio di mare in base al quale lo stesso 
stato adulto è raggiunto partendo dal germe normale completo, quanto da 
una metà, dopo una divisione sperimentale, fu, come è noto, la base per la 
formulazione delle teorie vitalistiche di Driesch. Il VON BERTALANFFY vede 
nell’autoregolazione globale, risultante dalle caratteristiche sopra descritte del 
sistema aperto in stato stazionario, la spiegazione di questo fondamentale 
aspetto degli organismi. L’accrescimento è infatti il risultato del metabolismo, 
inteso come una contro-azione tra anabolismo e catabolismo dei materiali di 
costruzione, e il rapporto relativo è regolabile in termini di autoregolazione 
dei sistemi considerati. 

Tl modello del sistema aperto in stato stazionario conduce, come s’é visto, 
all’ammissione di meccanismi autoregolativi impliciti nel sistema. Ciò si accorda 
con le caratteristiche presentate dagli organismi. Ma un modello di autocon- 
trollo valido per gli organismi è presentato anche dalla cibernetica. Si tratta 
dei meccanismi a « feed-back ». Il von BERTALANFFY non vede alcuna contrad- 
dizione tra i due modelli. Il problema è trattato ripetutamente (!5) e il Nostro 
conclude che certamente vi sono molti fenomeni biologici inquadrabili nello 
schema del « feed-back », dai fenomeni di omeostasi (la termoregolazione, ad 
esempio, negli animali a sangue caldo), ai meccanismi del sistema nervoso. 
Questi, per potersi verificare, presuppongono però arrangiamenti strutturali. 
Negli organismi esistono invece molte regolazioni di natura diversa. Esse non 
si basano su strutture preesistenti in qualche modo, ma si fondano sull’in- 
terazione dinamica dei processi. Sono esse, come dice il von BERTALANFFY, 
regolazioni primarie e sono spiegabili in termini di stati stazionari 0 di tendenza 
allo stato stazionario in sistemi aperti. Le regolazioni secondarie, del tipo « feed- 
back », sovrimporrebbero successivamente, sia nel senso dello sviluppo embrio- 
nale, quanto in quello evolutivo, alla regolazione dinamica dei componenti, 
arrangiamenti fissi e condizioni di costrizione, rendendo il sistema e le sue 
parti più efficienti. 

In Problems of Life (19) von BERTALANFYY scrive: «A living organism is a 


(5) Vedi in particolare: L. VON BERTALANNFY: Human Biology, 23, 346 (1951). 
(16) L. von BERTALANFFY: Problems of Life. An Evaluation of Modern Biological 
Thought (Londra, 1952). Ed. originale: Das biologische Weltbild (Berna, 1949). 
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hierarchy of open systems maintaining itself in a steady state due to its inherent 
system conditions ». Anche altrove ritorna su quello che chiama: « ordine ge- 
rarchico » secondo cui, in molti sistemi, i componenti sono essi stessi sistemi 
di un ulteriore ordine più basso, significando con cid che ogni Q; è esso stesso 
un sistema di elementi R, definibile pure da un sistema di equazioni cinetiche 
quali quelle già usate. A parte l’espressione « ordine gerarchico » che, se non 
definita rigorosamente, rischia di riportare il discorso su di un piano quali 
tativo, questo concetto di « sistemi di sistemi » non è stato approfondito dal 
VON BERTALANFFY quanto meriterebbe, e non consente pertanto qui un lungo 
esame. 


5. — Formalizzabilità e accordo con l’esperienza del modello. 


Occorre adesso tentare di rispondere alla domanda posta all’inizio. È questo 
modello proponibile? Si può confrontare con le tre condizioni poste. Si risponde 
subito facilmente alla prima. Il modello è certamente suscettibile di una for- 
mulazione matematica. Sia sotto il profilo cinetico, sia sotto quello termodina- 
mico è utilizzabile il formalismo matematico per corrispondere al quale, in 
fondo, il modello è stato formulato. 

Anche la seconda condizione non richiede un lunghissimo discorso. Il mo- 
dello non è contraddetto da alcun fatto sperimentale specifico. Le stesse idee 
fondamentali su cui si basa oggi la biologia non sono in contrasto con esso. 
La teoria cellulare vi si inquadra perfettamente; anzi, si potrebbe dire, è espli- 
citamente richiesta dalla concezione di sistema composto di elementi, descritto 
dalle equazioni cinetiche scritte sopra. L’accrescimento individuale e l’evolu- 
zione, cioè l'aumento dell’ordine è ammesso dalla termodinamica dei processi 
irreversibili, purchè siano rispettate le condizioni cui si è accennato sopra. 
Un discorso un po’ diverso si potrebbe fare per la genetica e i problemi rela- 
tivi. Il von BERTALANFFY pur avendo studiato geni e cromosomi alla luce del 
modello considerato (vedi oltre), non ci pare abbia affrontato sistematicamente 
il problema di come dei sistemi aperti (geni e cromosomi) possono controllare 
il sistema più ampio (organismo) di cui sono elementi. È vero che (vedi sopra) 
secondo VON BERTALANFFY ogni mutamento di qualche quantità Q, causa un 
cambiamento in tutti gli altri elementi e quindi nel sistema nel suo complesso. 
Si può anche ammettere che il Q,, rappresentante un’unità genetica, appaia 
nell’equazione con parametri tali, rispetto agli altri elementi, che sue even- 
tuali modificazioni abbiano un'influenza particolarmente ampia su altri ele- 
menti, o su alcuni tra questi. Però una simile formulazione è eccessivamente 
generica e insufficiente di fronte all’importanza ed alla regolarità che gli eventi 
genetici hanno negli organismi viventi. In ogni modo si tratta qui di una insuf- 
ficienza piuttosto che di una smentita. 
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Non si può d’altra parte non sottolineare le non trascurabili conferme e 
concordanze che su singoli punti il modello presentato ha trovato. Ricorde- 
remo qui alcuni dati. 

Per quanto riguarda gli aspetti termodinamici già la termodinamica ordi- 
naria permetteva di sviluppare, in casi particolari, concetti che saranno piena- 
mente compresi alla luce della termodinamica nuova e del modello del « sistema 
aperto ». La termodinamica classica mostrava che talora è possibile trascurare, 
nelle equazioni relative, il termine della variazione dell’entropia. In questi casi 
la produzione di entropia dovuta ad una reazione chimica è semplicemente 
proporzionale al calore di reazione. In tale approssimazione la produzione di 
entropia di un organismo vivente può essere misurata dal suo metabolismo 
in base alle registrazioni calorimetriche. Una reazione accoppiata del tipo 
Av<0 sarà perciò equivalente ad una reazione endotermica e la sua pre- 
senza si manifesterà mediante la diminuzione del calore totale che il sistema 
dà al mondo esterno. Molti autori hanno lavorato a questo problema, parti- 
colarmente nella fase embrionale, quando, per i processi di organizzazione, 
reazioni accoppiate devono giocare un ruolo particolarmente importante (!). 

Gli aspetti biologici di un altro punto della termodinamica sono stati svi- 
luppati dal Lrpman (18). Per reazioni vicine all’equilibrio si ha [Al Rae 
(R è la costante dei gas). Se laffinita di una data reazione chimica è grande 
e questa relazione non è soddisfatta, la reazione può spesso esser divisa in un 
certo numero di reazioni elementari, aventi ciascuna affinità sufficientemente 
piccola per giustificare l'applicazione delle leggi della termodinamica (le così 
dette leggi fenomenologiche lineari). Lo sviluppo di questi casi interessa il 
biologo proprio perchè costui ha spesso a che fare con reazioni chimiche che 
usualmente hanno luogo in molti passi, ciascuno dei quali è vicino alla rever- 
sibilità. 

Ricordiamo infine che abbastanza recentemente HEARON (!°) ha messo in 
evidenza l’importanza, per i processi biologici, dell’accoppiamento tra diffu- 
sione e reazioni chimiche, studiato dalla termodinamica e in base a cui si può 
capire, in una cellula vivente, una distribuzione stazionaria (non di equilibrio) 
di materia che differisca da quella ambientale in base ai rapporti metabolici 
della cellula stessa. 

Anche la validità delle equazioni cinetiche usate dal von BERTALANFFY è 
stata verificata in biologia. In capitoli diversi, studiosi diversi, per risolvere 
problemi diversi, hanno usato questo particolare tipo di equazioni differenziali. 

Un caso notissimo e di particolare importanza è costituito dagli studi del 


(7) Vedi in particolare: J. NEEDHAM: Chem. Embryol., 3 vols. (Cambridge, 1931). 
(18) F. Lipman: Adv. Enzym., 1, 99 (1941). 
(19) J. Z. Hearon: Bull. Math. Biophys., 12, 135 (1950). 
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Lotka, ma particolarmente del VOLTERRA (2°), sulla dinamica delle popola- 
zioni. Non è evidentemente qui possibile soffermarsi su questi problemi, ma 
è sufficiente ricordare che l’accrescimento teorico di una specie, supposta iso- 
lata e non impedita da alcuna limitazione, si sviluppa proprio secondo Pequa- 
zione dQ/dt=a,Q che, come s'è visto, integrata da Q = @ exp[at], indicante, 
per Q = numero degli individui al tempo ¢ e @, = numero iniziale degli indi- 
vidui, un aumento secondo una curva esponenziale. Con esponente negativo 
descrive invece il decremento di una popolazione in cui la mortalità sia mag- 
giore della natalità, ma si applica anche all’uccisione di batteri da parte di 
disinfettanti o alla perdita di sostanza corporea in un animale che digiuni. 
In pratica però una popolazione, anche supposta isolata, oltre certi limiti, 
trova sempre ostacoli al suo sviluppo che pertanto va gradatamente dimi- 
nuendo, sì da descrivere piuttosto che una curva esponenziale, una curva logi- 
stica. Anche questa curva è descritta analiticamente da un’altra equazione già 
vista, e cioè: dQ/dt=a,Q+a4,,Q?. La stessa equazione può descrivere anche 
certi casi di crescita organica e di relativo aumento di determinati organi e 
attività fisiologiche rispetto alle dimensioni del corpo. Ricordiamo infine, per 
concludere su questo punto, che anche il caso dell’associazione di due specie, 
di cui una si nutra dell’altra, è stato descritto dal VOLTERRA con un sistema 
di equazioni che fa pure parte della descrizione generale di un sistema: 


dq 

| "i == h(@1, (2) ’ 
dQ. 

| ay (050) 


Trascuriamo di esaminare lo sviluppo del sistema nel caso specifico consi- 
derato. Basta ricordare che ne derivano delle fluttuazioni regolari per cui lau- 
mento numerico della specie predatrice determina una riduzione della specie 
predata; onde, a causa della successiva mancanza di alimenti, si ha una ridu- 
zione anche della specie predatrice. Questo determina una ripresa della specie 
predata che fornisce nuovo alimento, e così via. L'esperienza ha confermato 
queste fluttuazioni. 

Per concludere questo punto, ricordiamo ancora che le stesse equazioni 
sono state usate dallo SPIEGELMAN (2!) per la cinetica dei processi cellulari e 
per la teoria della competizione entro un organismo. Il WERNER (22), infine, 
ha usato, come legge di base per l’azione farmacodinamica, un sistema di equa- 
zioni simili, ma più generali. Egli ha considerato infatti il sistema come con- 

(2°) Per un prospetto generale di questi studi vedi: U. D’Ancona: La lotta per 


l’esistenza (Torino, 1942). Questo volume è stato preparato dall’autore proprio su 
suggerimento di von BERTALANFFY. 


(21) S. SPIEGELMAN: Quart. Rev. Biol., 20, 121 (1945). 
(22) G. WERNER: Site. Akad. Wiss. Wien Math.-Naturw. KI., Abt. Ia, 156, 457 (1947). 
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tinuo, usando perciò equazioni alle derivate parziali rispetto alle variabili , 
y, 2 e t. Da tale legge di base si ricavano le varie leggi dell’azione della droga 
introducendo le condizioni speciali corrispondenti. 


6. — Operatività del modello. 


Non resta che esaminare l’ultima condizione cui era stato accennato: se 
cioè il modello può essere operativo, suscettibile cioè di trarne conseguenze 
verificabili. La risposta può esser data tenendo presente che il von BERTALANFEY, 
dallo schema generale di cui si parla, ha ritenuto di trarre alcune conseguenze 
che hanno poi avuto una conferma sperimentale. Sulla operatività di questo 
modello del resto il von BERTALANFFY non ha dubbi. « Ciò che è fondamentale 
per una vera biofisica, fisica dell'organismo vivente, è una teoria dei sistemi 
aperti e degli stati stazionari » (?*), scrive. 

Nel quadro della visione dell'organismo vivente come sistema aperto in 
stato stazionario, si collocano gli interessantissimi studi del von BERTALANFFY 
tendenti a stabilire connessioni tra i due aspetti fondamentali degli organismi 
viventi: il metabolismo e l'accrescimento (24). Il metabolismo varia ovvia- 
mente con le dimensioni del corpo. Misurando il metabolismo in base al con- 
sumo di O, per chilogrammo di peso del corpo si era già notato che questo 
rapporto diminuiva coll’aumentare del peso del corpo e che il metabolismo 
aumentava circa in proporzione all’aumento della superficie, anzichè del peso. 
La superficie di un organismo (S) si può ricavare, come è noto, mediante la 
formula di Mech, dal peso (P) mediante la relazione approssimativa: S =bP? (*), 
dove b è un parametro legato alla particolare forma della specie considerata. 
Si hanno naturalmente anche formule rese più esatte mediante alcune corre- 
zioni, per esempio quella di Dubois valida nel caso dell’uomo, ma la (*) è suf- 
ficiente per il ragionamento che segue. Se il metabolismo (M) varia col variare 
della superficie, si avrà: M — pP?, È facile vedere come questa formula si 
inquadri perfettamente nelle equazioni cinetiche dei sistemi. Consideriamo due 
elementi di un sistema, @, e 0, nel caso più semplice in cui tutti i coefficienti 
sono uguali a zero tranne il primo, si avrà: 


dQ: G 

tip oe dio Le [Q1 = € exp Pad 
i da cui integrando: c bs 

dQ. | Qs = e exp [dot] . 

ath = Ag); 


(23) L. von BERTALANFFY: Scientia, 11, 361 (1954). 

(24) L. von BERTALANFFY: Quantitative Laws in Metabolism and Growth. Ciclostilato 
senza indicazioni di luogo e di data, ma redatto nel 1956. Vedi anche, sempre di 
von BERTALANFFY: Nature, 163 (1946) e Amer. Nat. , €5, 111 (1951). 
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Risolvendo rispetto a È: 


i log, Q,—logse, log, d. — 108.6, 
==, SS = = dna oaeoe af , 


Gli Ay 
da cui 
di x di 
log, (0A gra log, Qt log, (Ge log. Co 3 
= (by As 
cioè 
4 log Oe la.) | 1 o Cy 
log, da = og, ds ae Og. cl1%) , 
2 
e ponendo 
a Cy 
A = ei (©) = o , 
As C; 


e passando dai logaritmi ai numeri: 
i yA 
Oe 


che corrisponde alla M—=bP?. In questo particolare caso A =3, ma non era 
evidentemente da escludere che questo valore non fosse universalmente valido, 
tanto più che la relazione metabolismo-superficie era in origine una relazione 
empirica, formulata in base a un numero relativamente ristretto di casi con- 
siderati. Il primo passo del von BERTALANFFY fu di pensare ad una genera- 
lizzazione. Dal 1941 von BERTALANFFY ed altri cercarono sperimentalmente 
come si comportasse la relazione considerata in molteplici gruppi animali e 
risultò che effettivamente l’esponente 3 non valeva in tuttii casi. Valeva, salvo 
qualche correzione, per i mammiferi, i pesci, in genere per i vertebrati pecilo- 
termi (per i quali, come è noto, non v’è termoregolazione) e per qualche inver- 
tebrato, ma non per gli altri animali. Si doveva perciò correggere la relazione: 
M=bP". x assumeva valori tra 2 e 1. Si potevano quindi delineare tre gruppi 
metabolici: il primo composto dagli animali con metabolismo in relazione alla 
superficie (x= 3), il secondo (larve di insetti, lumache, vermi) composto dagli 
animali in cui il metabolismo è invece in relazione direttamente al peso 
(x=1), il terzo (platelminti e alcune lumache) con valori del metabolismo 
intermedi a quelli dei due primi gruppi. L’interpretazione di questi tre tipi 
metabolici resta ancora aperta, comunque il von BERTALANFFY utilizzò questi 
risultati per studiare l'accrescimento. La concezione dell'organismo come si- 
stema aperto in stato stazionario dette l’avvio. Lo stato stazionario implica 
che il risultato globale di crescita o di equilibrio sia il risultato di processi 
dinamici continui, costruttivi e distruttivi. L’accrescimento deve perciò essere 


boi) 
lord 
wD 


GLI ORGANISMI BIOLOGICI COME « SISTEMI APERTI STAZIONARI ) ECC. is) 


il risultato della differenza tra i processi anabolici e quelli catabolici: dP/dt = 
— yy P™—kP", dove n indica la costante dell’anabolismo, k la costante del 
catabolismo e m ed » indicano una qualche dipendenza col peso del corpo P. 
Questa concezione della crescita come contro-azione di anabolismo e di cata- 
bolismo, per cui, quando l’anabolismo non prevale più, si raggiunge lo stato 
stazionario, ha permesso, come si vedrà, di determinare i vari tipi di accresci- 
mento in rapporto ai tipi metabolici sopra considerati. La relazione sopra 
scritta è evidentemente approssimata in notevole misura se si pensa alla sem- 
plificazione, che essa implica, di fenomeni che sono invece di una grande com- 
plessità. Il fatto che se ne traggano risultati soddisfacenti depone perciò a 
favore del modello di ragionamento. 

Il von BERTALANFFY osserva che i processi catabolici significano: biochimi- 
camente, la degradazione dei materiali costruenti, particolarmente delle pro- 
teine, come è dimostrato dalle tecniche isotopiche; citologicamente, il rinnovo 
delle cellule. Il catabolismo si può pertanto assumere, in prima approssimazione, 
come direttamente proporzionale al peso. Perciò, anche considerando il fatto 
che l'equazione di cui sopra è poco sensibile a piccoli spostamenti di n» dal- 
l’unità, si può porre n=1. 

La situazione per l’anabolismo è diversa. La sintesi dei componenti cellu- 
Jari macromolecolari ha bisogno da un lato di costruire aminoacidi, zuccheri, 
fosfati, ecc. e dall'altro lato dell’energia fornita dai processi ossidativi. Del 
resto i risultati sperimentali indicano che, almeno negli animali superiori, vi 
è una regolare connessione tra respirazione, anabolismo e accrescimento. Perciò 
VON BERTALANFFY, con un’altra approssimazione, attribuisce all’esponente m 
i valori trovati per « nel caso dei tipi metabolici. Si avranno così, anche per 
l'accrescimento, tre casi,: m=3, m=1, 2<m<1. Le predizioni dei tre tipi 
di accrescimento in funzione dei tipi metabolici sono state sperimentate e trovate 
corrette in gran numero di casi, spesso in modo insospettato. Nel primo caso 
la proporzione dell’accrescimento diminuisce progressivamente secondo una 
curva sigmoide, finchè si raggiunge eventualmente’ lo stato stazionario (mam- 
miferi, con alcune restrizioni, pesci, alcuni invertebrati). Nel secondo caso 
(larve di insetti, ecc.) l’anabolismo ed il catabolismo procedono secondo lo 
stesso andamento e la proporzione dell’accrescimento tende ad aumentare 
indefinitamente (accrescimento esponenziale) senza raggiungere alcun stato 
stazionario. Il paradosso di un accrescimento indefinito è naturalmente solo 
apparente: l'accrescimento di una larva è infatti ad un certo punto tagliato 
dalla metamorfosi, quello di una lumaca di questo tipo dal ciclo stagionale, ecc. 
Nel terzo caso infine si hanno valori intermedi. 

Dalle curve di accrescimento ottenute in relazione a questo modello possono 
ricavarsi i valori della costante catabolica hk. Tali valori possono confrontarsi 
con quelli ottenuti sperimentalmente dalla distruzione delle proteine. La coinci- 
denza tra i valori ottenuti con i due metodi, secondo von BERTALANFFY, È 
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stata talora sorprendente, se si considera da un lato che il modello teorico 
di partenza è estremamente semplificato e dall’altro che l’errore nelle deter- 
minazioni fisiologiche della degradazione delle proteine è notevole. Ad esempio, 
nel 1938, il von BERTALANFFY calcolò il ritmo di degradazione delle proteine 
dalle curve di accrescimento dell’uomo, ottenendo il valore di 1.165 g/kg di 
peso del corpo al giorno. Nel 1949 SPRINSON e RITTENBERG fecero lo stesso 
calcolo partendo dai loro esperimenti con 15N e trovarono il valore di 3 g/kg 
di peso del corpo al giorno. L’accordo tra le due cifre non è molto, per la 
verità, ma tenendo conto della situazione è da considerare favorevolmente e 
quanto meno indicativo dell’ordine di grandezza del fenomeno. Non è qui il 
caso di discutere dettagliatamente questi problemi, nè la possibilità di appor- 
tare al modello correzioni, nè quelle già apportate, ma piuttosto di vedere 
altre conseguenze che il von BERTALANFFY ha tratto dalla concezione del- 
Vorganismo come sistema aperto in stato stazionario e che sono state poi con- 
fermate. 

Le considerazioni che precedono, ed altre, tutte legate alla teoria della 
crescita animale intesa come il risultato d’una azione antagonista di processi 
anabolici e catabolici (25) ha condotto il nostro autore ad interessanti conclu- 
sioni (26): 1) il catabolismo e la rigenerazione delle proteine hanno luogo a 
una velocità più grande di quanto non immagini la fisiologia classica (27); 
2) si ha sintesi e risintesi di aminoacidi a partire da ammoniaca e da catene 
non azotate; 3) nel caso di una alimentazione povera in proteine, l’azoto pro- 
veniente dai costituenti dei corpi catabolizzati è trattenuto e riutilizzato per 
la rigenerazione proteica. Questo spiegherebbe l’azione di risparmio proteico 
giocata dagli idrati di carbonio, per cui una dieta ricca appunto in idrati di 
carbonio necessita di un minore apporto in proteine. 

Scrive il von BERTALANFFY: « Queste conclusioni sono state verificate da 
successivi lavori con isotopi, particolarmente con !N. La tecnica degli iso- 
topi ha dimostrato il rapido rivolgimento delle proteine nel corpo animale, 
la rigenerazione delle protéine da un ,,metabolic pool” contenente azoto deri- 
vato dalle proteine degradate e l’utilizzazione dell’azoto ammoniacale del re- 
gime alimentare per la sintesi di proteine nelle diete povere di proteine » (28). 


(°°) Il von BERTALANFFY si sofferma su questa teoria particolarmente in: Expe- 
rientia, 4, 225 (1948) e nel saggio già citato in: Nature (Londra, 1949), p. 163. 

(26) L. von BERTALANFFY: Human Biology, 10, 181 (1938). 

(2°) Il von BERTALANFFY cita a conferma particolarmente: LEBLOND e STEVENS: 
Anat. Rec., 100, 357 (1948); SroREy e LEBLOND: Ann. N.Y. Acad. Sci., 58, 537 (1951); 
F. D. BERTALANFFY e LEBLOND: Anat. Rec., 115, 515 (1953); LEBLOND e WALKER:. 
Physiol. Rev., 36, 225 (1956) e L. von BERTALANFFY: Biophysik des FlieRgleichgewichts 
(Braunschweig, 1953) per mostrare che anche il rinnovo delle cellule procede spesso 
ad un ritmo inaspettatamente alto. 

(28) L. von BERTALANFFY: Scientia, 11, 361 (1954). 
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Il von BERTALANFFY Cita a sostegno della sua test i lavori con !5N di SPRINSON 
e RITTENBERG (°°). Effettivamente questo continuo dinamismo negli organismi 
viventi, interessante anche sostanze in apparenza inerti come la dentina o 
l’osso, così lontano dall’equilibrio inteso in senso classico, ma pure tendente 
o mantenente uno stato stazionario, ci pare di non trascurabile peso in favore 
della teoria. Tanto più che questo ordine di idee ha permesso di dedurre altre 
conseguenze tendenti ad illustrare ulteriormente il lavoro d’insieme della cel- 
lula e dell’organismo. Partendo da queste premesse, infatti, KUHN (°°) e VON 
BERTALANFFY (3!) indicarono il lavoro minimo richiesto per produrre una 
molecola d’una sostanza polimera a concentrazione e diversa dalla concentra- 
zione d’equilibrio c* con la relazione: W= RT log e/e*. I calcoli di SCHULZ (8°) 
esprimono approssimativamente il lavoro necessario per grammo di proteina 
con: dw/dt= r(RT/Mw)log(nw/K), ove r indica il rivolgimento di proteine per 
giorno, Mw il peso molecolare degli aminoacidi condensati in polimeri, mw il 
numero di molecole d’acqua per molecola di sostanza mono e polimera e K 
una costante d’equilibrio. Con opportuni valori ne deriva: dw/dt= 7-93 cal/g 
proteina. Il von BERTALANFFY attribuisce ad r valori medi ricavati dai lavori 
già citati di SPRINSON e RITTENBERG con l’ ¥N, e ottiene che è necessario 
da uno a due per cento del metabolismo basale per mantenere le proteine dei 
corpi in stato stazionario. La sintesi proteica richiede dunque una ben mo- 
desta parte del metabolismo totale. Stime fatte da Borsook (°°) e da LANG (34) 
con metodi diversi (esperienze con !5N ma basate su ipotesi diverse e sulla 
duplicazione del materiale nucleare durante la mitosi) corrispondono ai dati 
ottenuti con il metodo sopra descritto sia pure a grandi linee. Commenta il 
VON BERTALANFFY: « Una conseguenza è che sarà futile cercare differenze nella 
respirazione cellulare tra organi con alta o bassa attività mitotica. Rispetto 
alla crescita maligna segue che un approccio dovrebbe esservi nella linea della 
sintesi proteica piuttosto che in quella del metabolismo dell’energia, come ha 
proposto RONDONI » (°°). 

Infine, prima di concludere questa serie di considerazioni sull’operatività 
del modello dell’organismo vivente come sistema aperto stazionario, ricorde- 
remo un'ultima deduzione del von BERTALANFFY non priva di interesse. Come 
è noto il DRIESCH considerò l’automoltiplicazione delle unità biologiche elemen- 
tari (geni, cromosoni) come una seconda prova del vitalismo (la prima era 


(29) D. B. SPRINSON e D. RITTENBERG: Journ. Biol. Chem., 180, 707, 715 (1949). 
(89) W. Kunn: Ergebn. Enzymforsch., iy di (SS) 
(31) L. VON BERTALANFFY: Theoretische Biologie, vol. 2, già citato. 
(32) G. V. ScHuLz: Naturwiss., 37, 196 (1950). 
(33) H. BorsooK: Physiol. Rev., 30, 206 (1950). 
(34) K. Lana: Mikroskopische und chemische Organisation der Zelle. 2. Colloquium 
d. deutschen Ges. f. physiol. Chemie (Berlin, 1952). 
(89) L. von BERTALANFFY: Scientia, 11, 361 (1954), già citato. 
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quella dell’« equifinalità » di cui si è già parlato). Ora il VON BERTALANFFY, 
partendo dalle sue idee, propose uno schema di autoduplicazione di tali unità 
biologiche elementari basato sull’intuizione che queste unità non fossero corpi 
stabili ed inerti, ma sistemi metabolizzati i quali, proprio mediante il meta- 
bolismo, si mantengono in uno stato stazionario crescendo per fissazione di 
gruppi molecolari appropriati su una matrice e subendo processi catabolici (88). 

Lavori di LEBLOND e collaboratori con il radiofosforo (**) hanno mostrato 
che non si ha scambio di acido desossiribonucleico (DNA) nell’interfase, ma 
bensì un rimpiazzamento completo del DNA nello stadio premitotico, onde 
soltanto il DNA nuovo passa ai nuclei figli. Ciò sarebbe dimostrato dal fatto 
che nei tessuti con attività mitotica bassa, media, elevata (fegato, polmoni, 
mucosa intestinale) la quantità di DNA P formata è circa doppia di quella 
che corrisponderebbe al numero delle cellule neoformate. Ciò mostrerebbe 
appunto, secondo il von BERTALANFFY, che il DNA è catabolizzato e rigene- 
rato durante la mitosi. Ma non possiamo purtroppo entrare nel merito di questo 
problema peraltro quanto mai importante. 

Questi cenni sono sufficienti, ci pare, per concludere sulla operatività del 
modello proposto. Ovviamente qui non interessava tanto entrare nel merito 
delle singole previsioni e delle relative conferme. Qui si voleva soltanto mo- 
strare che dal modello si possono ricavare previsioni non solo qualitative, ma 
anche quantitative su singoli eventi, tali da poter essere sottoposte al con- 
fronto con l’esperienza che le verificherà o meno. Il modello generale proposto 
è pertanto operativo. 


7. — Conclusione. 


Il modello dell’organismo vivente come «sistema aperto stazionario », in 
base a quanto è stato detto nella Sezione 2 è dunque seriamente proponibile. 
Come si è detto la proponibilità soddisfa alcune condizioni fondamentali, ma 
non è di per sè un giudizio di rispondenza, di completezza, di validità. Essa 
pertanto non chiude, ma apre un discorso critico impegnato. Restano infatti 
da affrontare numerosi problemi. Innanzi tutto bisogna dare una risposta al 
quesito se il modello sia completo o no. Rispondere a questa domanda signi- 
fica in fondo confrontare il modello con tutta la problematica biologica, oltre 
ai punti qui accennati. Si è fatto riferimento, per esempio, alla necessità di 


(86) L. von BERTALANFFY: Naturwiss., 32, 46 (1944). 
(*?) C. E. STEVENS, R. Daouste C. P. LEBLOND: Journ. Biol. Chem., 202, 177 (1953); 


C. E. STEVENS, R. Daoust e C. P. LeBLOND: Can. Journ. Med. Sci., 31, 263 (1953); 
R. DaAousr, F. D. BERTALANFFY e C. P. LEBLOND: Journ. Biol. Chem.. 207, 405 (1954). 
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esaminare la corrispondenza tra il modello proposto e la genetica, e, soprat- 
tutto, i problemi legati alla modalità con cui il gene agisce nel determinare 
i caratteri che controlla, ma si potrebbero aggiungere altri punti su cui certa- 
mente è fondamentale una verifica dello schema: dai problemi della contra- 
zione del muscolo ai modelli matematici della cellula di Rashevsky; dai più 
recenti studi di embriologia a tutta la vasta problematica della biochimica, ecc. 

E non si è parlato di tutti gli aspetti più strettamente microfisici, di così 
grande interesse per la biologia anche se ci si volesse limitare, nella zona che 
sta tra il livello molecolare e quello atomico, ai loro aspetti teorici più generali 
quali quelli suggeriti panoramicamente dal MONDELLA nel suo bel lavoro re- 
cente (35). 

In ultima analisi, del resto, solo quando i Q, gli X e gli J delle equazioni 
esaminate potranno acquistare dei significati e dei valori ben precisi la teoria 
avrà assolto il suo compito. 

Qui intanto volevamo solo segnalare l’interesse di questo ordine di idee e 
mostrare come oggi siano possibili in biologia, in linea di massima, modelli 
teorici tali da promuovere la ricerca e da rendere utilizzabile il formalismo 
matematico. 


(88) F. MonpELLA: Il Pensiero, 3 (1957). 
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1. — Sostituire la Tabella I a pag. 550 con la seguente. 


TABLE I. — Data on cosmic abundances, including the primary cosmic radiation (beyond 
Z=10, only the most abundant species are listed). 


N 
[©] 


b | e d 
1 H 4 210 12.00 | 620 180 
2 He 3 = 10? 11.06 | 90 20 
3 Li 100 0.86 |) 
4 Be | 20 2.43 | 3 0.4 
5 B 24 [5.0] | 
6 © 3.5 -108 8.56 | | 
if N 6.6 -106 7.98 6 0.8 
8 O 2.15.10? 9.00 | | 
9 F 11600 | 
10 Ne 8.6 «105 8.81 1) | 
14 Si Pah awe 7.60 | | | 
18 A PO Ay okey TA | 2:5 0.3 
20 Ca 4.9 -104 6.38 | | 
26 Te ose a6 e108 6.76 | 
(a) H. E. Suess and H. C. UREyY, Rev. Mod. Phys. 28, 53 (1956). (These values are re. | 
ferred to silicon=1-10°). | 
(6) log N for the sun, except that for He and Ne we give their ratios to H for stars | 
(after ALLER). 
(c) Vertical cosmic ray flux (in particles per m? ster s (peters) at A= 41° [8]. | 


| (d) Vertical cosmic ray flux (in peters) at 2= 0° [S]. 
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2° — Sostituire le References di pag. 555 con le seguenti. 
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